Changeset 700 for trunk/MESOSCALE/LMD_MM_MARS/SRC

Timestamp:

Jun 7, 2012, 11:12:36 AM (12 years ago)

Author:

acolaitis

Message:

Mesoscale: minor modif so that you can run a nested simulation with radt=1 without the -DSPECIAL_NEST_SAVE flag. (also replaced the initialisation q2 at start from 0 to 1.E-6)

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F (modified) (6 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F

@@ -460 +460 @@
 ENDIF
 #else
+IF( .NOT. ALLOCATED( dp_save  ) ) THEN
 ALLOCATE(dp_save(ngrid)) !! here are the arrays that would be useful to save physics tendencies
 ALLOCATE(du_save(ngrid,nlayer))
@@ -465 +466 @@
 ALLOCATE(dt_save(ngrid,nlayer))
 ALLOCATE(dq_save(ngrid,nlayer,nq))
+ENDIF
 dp_save(:)=0.    !! initialize these arrays ...
 du_save(:,:)=0.
@@ -473 +475 @@
 #ifndef NORESTART
 ! Restart save arrays
+IF( .NOT. ALLOCATED( save_tsoil_restart  ) ) THEN
 ALLOCATE(save_tsoil_restart(ngrid,nsoil))
 ALLOCATE(save_co2ice_restart(ngrid))
@@ -478 +481 @@
 ALLOCATE(save_qsurf_restart(ngrid,nq))
 ALLOCATE(save_tsurf_restart(ngrid))
+ENDIF
 save_tsoil_restart(:,:)=0.
 save_co2ice_restart(:)=0.
@@ -484 +488 @@
 save_tsurf_restart(:)=0.
 #ifdef NEWPHYS
+IF( .NOT. ALLOCATED( save_wstar_restart  ) ) THEN
 ALLOCATE(save_wstar_restart(ngrid))
 ALLOCATE(save_fluxrad_restart(ngrid))
+ENDIF
 save_wstar_restart(:)=0.
 save_fluxrad_restart(:)=0.
@@ -1138 +1144 @@
 !-------------------------!
 IF (.not. restart) THEN
-   q2_val(:)=0.      !PBL wind variance
+   q2_val(:)=1.E-6      !PBL wind variance
 #ifdef NEWPHYS
    fluxrad_val=0.