Changeset 1590 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Aug 30, 2016, 3:04:29 PM (8 years ago)

Author:

aslmd

Message:

moved interface for outputs from module_lmd_driver to update_outputs_physiq_mod. updated the nophysics case by adding a case dynphy_wrf_void_lmd_new which will serve as a template for a physiq interface

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• comm_wrf.F90 (modified) (2 diffs)
• physiq_mod.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/comm_wrf.F90

@@ -3 +3 @@
 !!
 !! This module is useful to output fields in WRF style
+!! -- useful for diagnostics that are not already in modules
 !!
 
@@ -23 +24 @@
   integer,intent(in) :: ngrid ! number of atmospheric columns
   integer,intent(in) :: nlayer ! number of atmospheric layers
-  allocate(comm_HR_SW(ngrid,nlayer))
-  allocate(comm_HR_LW(ngrid,nlayer))
-  allocate(comm_SWDOWNZ(ngrid))
-  allocate(comm_TAU_DUST(ngrid))
-  allocate(comm_RDUST(ngrid,nlayer))
-  allocate(comm_QSURFDUST(ngrid))
-  allocate(comm_MTOT(ngrid))
-  allocate(comm_ICETOT(ngrid))
-  allocate(comm_VMR_ICE(ngrid,nlayer))
-  allocate(comm_TAU_ICE(ngrid))
-  allocate(comm_RICE(ngrid,nlayer))
+  if (.not.allocated(comm_HR_SW)) allocate(comm_HR_SW(ngrid,nlayer))
+  if (.not.allocated(comm_HR_LW)) allocate(comm_HR_LW(ngrid,nlayer))
+  if (.not.allocated(comm_SWDOWNZ)) allocate(comm_SWDOWNZ(ngrid))
+  if (.not.allocated(comm_TAU_DUST)) allocate(comm_TAU_DUST(ngrid))
+  if (.not.allocated(comm_RDUST)) allocate(comm_RDUST(ngrid,nlayer))
+  if (.not.allocated(comm_QSURFDUST)) allocate(comm_QSURFDUST(ngrid))
+  if (.not.allocated(comm_MTOT)) allocate(comm_MTOT(ngrid))
+  if (.not.allocated(comm_ICETOT)) allocate(comm_ICETOT(ngrid))
+  if (.not.allocated(comm_VMR_ICE)) allocate(comm_VMR_ICE(ngrid,nlayer))
+  if (.not.allocated(comm_TAU_ICE)) allocate(comm_TAU_ICE(ngrid))
+  if (.not.allocated(comm_RICE)) allocate(comm_RICE(ngrid,nlayer))
   end subroutine allocate_comm_wrf
 

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -1993 +1993 @@
       !! --> example : hfmax_th, zmax_th
 
+      CALL allocate_comm_wrf(ngrid,nlayer)
+
       !state  real  HR_SW     ikj   misc  1  -  h  "HR_SW"     "HEATING RATE SW"                 "K/s"
       comm_HR_SW(1:ngrid,1:nlayer) = zdtsw(1:ngrid,1:nlayer)