Changeset 976 for codes

Timestamp:

09/25/19 16:40:08 (5 years ago)

Author:

jisesh

Message:

devel: Added variable Teq in run.def

File:

• codes/icosagcm/devel/src/dcmip/dcmip_initial_conditions_test_1_2_3_v5.f90 (modified) (6 diffs)

codes/icosagcm/devel/src/dcmip/dcmip_initial_conditions_test_1_2_3_v5.f90

@@ -1026 +1026 @@
 
 SUBROUTINE test2_schaer_mountain (lon,lat,p,z,zcoords,hybrid_eta,hyam,hybm,shear,u,v,w,t,phis,ps,rho,q)
-
+USE getin_mod, ONLY : getin
 IMPLICIT NONE
 !-----------------------------------------------------------------------
@@ -1074 +1074 @@
                                 as      = a/X,                  &       ! New Radius of small Earth     
                                 ueq     = 20.d0,                &       ! Reference Velocity 
-                                Teq     = 300.d0,               &       ! Temperature at Equator    
                                 Peq     = 100000.d0,            &       ! Reference PS at Equator
                                 ztop    = 30000.d0,             &       ! Model Top       
@@ -1090 +1089 @@
       real(rstd) :: c                                                   ! Shear
 
+      REAL(rstd) :: Teq                 = 300.d0   ! default Temperature at Equator, can be modified in run.def
+
+      CALL getin("dcmip2_Teq",Teq)
 !-----------------------------------------------------------------------
 !    PHIS (surface geopotential)
@@ -1197 +1199 @@
 !==========
 SUBROUTINE test3_gravity_wave (X,lon,lat,p,z,zcoords,u,v,w,t,phis,ps,rho,q)
-
+USE getin_mod, ONLY : getin
 IMPLICIT NONE
 !-----------------------------------------------------------------------
@@ -1232 +1234 @@
                                 u0      = 20.d0,                &       ! Reference Velocity 
 !                               u0      = 0.d0,                 &       ! FIXME : no zonal wind for NH tests
-                                Teq     = 300.d0,               &       ! Temperature at Equator    
                                 Peq     = 100000.d0,            &       ! Reference PS at Equator
                                 ztop    = 10000.d0,             &       ! Model Top       
@@ -1253 +1254 @@
       real(rstd) :: theta_pert                                          ! Pot-temp perturbation
 
+      REAL(rstd) :: Teq                 = 300.d0   ! default Temperature at Equator, can be modified in run.def
+
+      CALL getin("dcmip2_Teq",Teq)
+
       as = a/X