Changeset 957 for codes

Timestamp:

07/16/19 01:45:04 (5 years ago)

Author:

dubos

Message:

devel : fix DCMIP31 when scale factor X=1

File:

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

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

@@ -1230 +1230 @@
         real(rstd), parameter ::        &
                                 Om      = 0.d0,                 &       ! Rotation Rate of Earth
-                                u0      = 20.d0,                &       ! Reference Velocity 
-!                               u0      = 0.d0,                 &       ! FIXME : no zonal wind for NH tests
+!                               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       
                                 lambdac = 2.d0*pi/3.d0,         &       ! Lon of Pert Center
-                                d       = 5000.d0,              &       ! Width for Pert
+                                d       = 625000.d0,            &       ! Width for Pert for X=1
                                 phic    = 0.d0,                 &       ! Lat of Pert Center
                                 delta_theta = 1.d0,             &       ! Max Amplitude of Pert
@@ -1333 +1333 @@
         r  = as * ACOS (sin_tmp + cos_tmp*cos(lon-lambdac)) 
 
-        s = (d**2)/(d**2 + r**2)
+        s = ((d/X)**2)/((d/X)**2 + r**2)
 
         theta_pert = delta_theta*s*sin(2.d0*pi*height/Lz)