coords_n_const.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	coords_n_const

Function/Subroutine Documentation

coords_n_const()

subroutine coords_n_const

(

)

Definition at line 41 of file coords_n_const.f90.

  USE mod_input
  USE mod_setup
  IMPLICIT NONE
  INTEGER STATUS, I

!===============================================================================!
! OPEN FILES
! UNLESS OTHERWISE SPECIFED THE ROUTINES CALLED HERE ARE IN mod_input.F
!===============================================================================!
    
  SELECT CASE(startup_type)
     !=================================================
     ! HOTSTART 
  CASE("hotstart")
     !=================================================
     if(dbg_set(dbg_log)) then 
        WRITE(ipt,*)'!               READING FILES FOR HOTSTART                       !'
        WRITE(ipt,*)'!                                                                !'
     end if
               
     ! ALLOCATE SPACE FOR THE LOCAL GRID DATA
     ALLOCATE(y_lcl(0:mt),stat=status)
     IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE Y_LCL")
     y_lcl=0.0_sp
     ALLOCATE(x_lcl(0:mt),stat=status)
     IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE X_LCL")
     x_lcl=0.0_sp
     ALLOCATE(h_lcl(0:mt),stat=status)
     IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE H_LCL")
     h_lcl=0.0_sp

     CALL load_restart_coords(x_lcl,y_lcl)
     CALL coordinate_units(x_lcl,y_lcl)
     CALL setup_center_coords


     DEALLOCATE(x_lcl)
     DEALLOCATE(y_lcl)


     CALL load_restart_depth(h_lcl)
     CALL setup_depth
     DEALLOCATE(h_lcl) ! COULD BE LOADED DIRECTLY - MUST SET MAX/MIN

     CALL load_restart_coriolis(cor)! LOAD DIRECTLY

    !!  ggao for equatoral min (4deg)
!     IF(.NOT. EQUATOR_BETA_PLANE)THEN
!      WHERE(COR <  1.e-5_SP .AND. COR > 0.0_SP) COR =  1.e-5_SP
!      WHERE(COR > -1.e-5_SP .AND. COR < 0.0_SP) COR = -1.e-5_SP
!     END IF 

     CALL setup_gravity
     
     CALL load_restart_sponge(cc_sponge) ! LOAD DIRECTLY
     
     stype = stype_restart
     CALL load_restart_sigma(z,z1) ! LOAD DIRECTLY TO ALL_VARS:Z,Z1
     CALL setup_sigma_derivatives

     CALL setup_obctypes

     CALL setup_horizontal_mixing_coefficient

     CALL setup_bottom_roughness

     !=================================================
     ! CRASHSTART 
  CASE("crashrestart")
     !=================================================
     if(dbg_set(dbg_log)) then 
        WRITE(ipt,*)'!              READING FILES FOR CRASHRESTART                    !'
        WRITE(ipt,*)'!                                                                !'
     end if
     ! ALLOCATE SPACE FOR THE LOCAL GRID DATA
     ALLOCATE(y_lcl(0:mt),stat=status)
     IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE Y_LCL")
     y_lcl=0.0_sp
     ALLOCATE(x_lcl(0:mt),stat=status)
     IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE X_LCL")
     x_lcl=0.0_sp
     ALLOCATE(h_lcl(0:mt),stat=status)
     IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE H_LCL")
     h_lcl=0.0_sp

          
     CALL load_restart_coords(x_lcl,y_lcl)
     CALL coordinate_units(x_lcl,y_lcl)
     CALL setup_center_coords


     DEALLOCATE(x_lcl)
     DEALLOCATE(y_lcl)

     
     CALL load_restart_depth(h_lcl)
     CALL setup_depth
     DEALLOCATE(h_lcl) ! COULD BE LOADED DIRECTLY - MUST SET MAX/MIN

     CALL load_restart_coriolis(cor)

     CALL setup_gravity
     
     CALL load_restart_sponge(cc_sponge) ! LOAD DIRECTLY TO ALL_VARS:CC_SPONGE


     stype = stype_restart
     CALL load_restart_sigma(z,z1) ! LOAD DIRECTLY TO ALL_VARS:Z,Z1
     CALL setup_sigma_derivatives

     CALL setup_obctypes

     CALL setup_horizontal_mixing_coefficient

     CALL setup_bottom_roughness

     !=================================================
     ! COLDSTART 
  CASE("coldstart")
     !=================================================
     if(dbg_set(dbg_log)) then 
        WRITE(ipt,*)'!              READING FILES FOR COLDSTART                       !'
        WRITE(ipt,*)'!                                                                !'
     end if
     
     IF(msr) THEN
        ! ALLOCATE SPACE FOR THE GLOBAL GRID DATA
        ALLOCATE(y_gbl(0:mgl),stat=status)
        IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE Y_GBL")
        y_gbl=0.0_sp
        ALLOCATE(x_gbl(0:mgl),stat=status)
        IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE X_GBL")
        x_gbl=0.0_sp
     END IF

    ! ALLOCATE SPACE FOR THE LOCAL GRID DATA
    ALLOCATE(y_lcl(0:mt),stat=status)
    IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE Y_LCL")
    y_lcl=0.0_sp
    ALLOCATE(x_lcl(0:mt),stat=status)
    IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE X_LCL")
    x_lcl=0.0_sp

    ALLOCATE(h_lcl(0:mt),stat=status)
    IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE H_LCL")
    h_lcl=0.0_sp

    ALLOCATE(c_lcl(0:mt),stat=status)
    IF (status /=0 ) CALL fatal_error("COULD NOT ALLOCATE C_LCL")
    c_lcl=0.0_sp

     
     CALL load_coldstart_coords(x_gbl,y_gbl,x_lcl,y_lcl)
     CALL coordinate_units(x_lcl,y_lcl)
     CALL setup_center_coords


     DEALLOCATE(x_lcl)
     DEALLOCATE(y_lcl)

     CALL load_coldstart_depth(x_gbl,y_gbl,h_lcl)
     CALL setup_depth

     DEALLOCATE(h_lcl)

     CALL load_coldstart_coriolis(x_gbl,y_gbl,c_lcl)
     CALL setup_coriolis

     DEALLOCATE(c_lcl)

     CALL setup_gravity

     ! THESE VARIABLES ARE NOT ALLOCATED YET
     CALL load_coldstart_sponge(x_gbl,y_gbl,nsponge,n_spg,r_spg,c_spg,x_spg,y_spg)
     CALL setup_sponge

 
!     CALL LOAD_COLDSTART_SIGMA
     CALL setup_sigma
     CALL setup_sigma_derivatives

     IF(msr) THEN
        DEALLOCATE(x_gbl)
        DEALLOCATE(y_gbl)
     END IF

     CALL setup_obctypes

     CALL setup_horizontal_mixing_coefficient

     CALL setup_bottom_roughness

  END SELECT
  
  

Here is the call graph for this function:

Here is the caller graph for this function: