Changeset 1759 for trunk

Timestamp:

Jul 26, 2017, 12:00:08 AM (7 years ago)

Author:

aslmd

Message:

MESOSCALE MARS and VENUS. corrected two bugs left from changes in r1755. 1) phour_ini was misplaced thus not set correctly. 2) pdtheta was not allocated.

File:

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

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

@@ -365 +365 @@
 ! ALLOCATE !
 !----------!
-!!-------------------------------------------------------------------------------!
-!! outputs:                                                                     !       
-!!    pdu(ngrid,nlayermx)        \                                              !
-!!    pdv(ngrid,nlayermx)         \  Temporal derivative of the corresponding   !
-!!    pdt(ngrid,nlayermx)         /  variables due to physical processes.       !
-!!    pdq(ngrid,nlayermx)        /                                              !
-!!    pdpsrf(ngrid)             /                                               !
-!!-------------------------------------------------------------------------------!
-!ALLOCATE(pdpsrf(ngrid))
-!ALLOCATE(pdu(ngrid,nlayer))
-!ALLOCATE(pdv(ngrid,nlayer))
-!ALLOCATE(pdt(ngrid,nlayer))
-!ALLOCATE(pdtheta(ngrid,nlayer))
-!ALLOCATE(pdq(ngrid,nlayer,nq))
 CALL allocate_interface(ngrid,nlayer,nq)
+ALLOCATE(pdtheta(ngrid,nlayer))
 !!!
 !!! BIG LOOP : 1. no call for physics, used saved values
@@ -440 +427 @@
             lct_input,lon_input,ls_input,&
             MY)
+  !! Set up initial time
+  phour_ini = JH_cur_split
   !! Fill planetary parameters in modules
   !! Values defined in the module_model_constants.F WRF routine
@@ -447 +436 @@
                      wdaysec,floor(JD_cur), &
                      1./reradius,g,r_d,cp,1)
-  !! Set up initial time
-  phour_ini = JH_cur_split
 ENDIF allocation_firstcall
 
@@ -775 +762 @@
     RTHBLTEN(i,kps:kpe,j) = pdtheta(subs,kps:kpe)
   ELSE
-   RTHBLTEN(i,kps:kpe,j) = zdtfi_omp(subs,kps:kpe) / exner(i,kps:kpe,j)
+    RTHBLTEN(i,kps:kpe,j) = zdtfi_omp(subs,kps:kpe) / exner(i,kps:kpe,j)
   ENDIF
     ! update surface pressure (cf CO2 cycle in physics)