Changeset 2086

Timestamp:

Feb 1, 2019, 12:21:21 PM (6 years ago)

Author:

mlefevre

Message:

Delete redifinition of atmospheric constant for ideal case

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_quarter_ss.F (modified) (3 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_quarter_ss.F

@@ -974 +974 @@
 ! local data
 
-      real p1000mb,cv,cp,r,cvpm,g
+      !real p1000mb,cv,cp,r,cvpm,g
 !****Mars
 !      parameter (p1000mb = 1.e+05, r = 287, cp = 1003., cv = cp-r, cvpm = -cv/cp, g=9.81 )
@@ -1021 +1021 @@
       p_surf = 100.*p_surf  ! convert to pascals
       qvf = 1. + rvovrd*qv_input(1) 
-      rho_surf = 1./((r/p1000mb)*th_surf*qvf*((p_surf/p1000mb)**cvpm))
+      rho_surf = 1./((r_d/p1000mb)*th_surf*qvf*((p_surf/p1000mb)**cvpm))
       pi_surf = (p_surf/p1000mb)**(rcp)
           !!!!!! rcp variable
@@ -1047 +1047 @@
 !                    - dz*(0.25*rho_surf+0.75*rho_input(1))*g*qvf1  !!! BEURK
 !                    - 0.5*dz*(rho_surf+rho_input(1))*g*qvf1  !! parce que couche 1 tres proche
-            rho_input(1) = 1./((r/p1000mb)*th_input(1)*qvf*((pm_input(1)/p1000mb)**cvpm))
+            rho_input(1) = 1./((r_d/p1000mb)*th_input(1)*qvf*((pm_input(1)/p1000mb)**cvpm))
                !!!!!!! rcp variable
                !rho_input(1) = 1./((r_input(1)/p1000mb)*th_input(1)*qvf*((pm_input(1)/p1000mb)**cvpm_input(1)))