Changeset 1578 for trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new

Timestamp:

Aug 22, 2016, 11:56:58 AM (8 years ago)

Author:

aslmd

Message:

Goal: make the mesoscale model compliant
with new interface philosophy and adapted
to easy interfacing with any physics

Commit 2:

• cleaned makemeso (and removed meso_dimphys.h_ref)
• added linking to phy_common in copy_model
• changed module_lmd_driver.F to reflect that
  • physiq and phys_state_var_init are now within a module
  • their list of arguments changed

NB:

mesoscale model still not working

with HEAD version
... but is now compiling!

need now to fix e.g. geometry

... this will be done within a call
... to the equivalent of iniphysiq
... to make lmd_driver planet-agnostic

retrocompatbility with v1520 is lost

File:

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

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

@@ -62 +62 @@
    USE update_outputs_physiq_mod !! to get outputs from physiq
    USE comm_wrf !! to get fields to be written from physiq
+   USE physiq_mod, only : physiq
+   USE phys_state_var_init_mod, only : phys_state_var_init
+   !!! ajouter init_time (voir iniphysiq. ne plus le faire dans la physiq)
 #endif
 
@@ -139 +142 @@
    REAL :: pday,ptime,MY  
    REAL :: phisfi_val
-   LOGICAL :: firstcall,lastcall,tracerdyn
+   LOGICAL :: firstcall,lastcall
    ! ---------- 
    REAL,DIMENSION(:,:),ALLOCATABLE :: pplev,pplay,pphi,pu,pv,pt,pw
@@ -375 +378 @@
 !    pdq(ngrid,nlayermx)        /                                               !
 !    pdpsrf(ngrid)             /                                                !
-!    tracerdyn                 call tracer in dynamical part of GCM ?           !
 !-------------------------------------------------------------------------------!
 ALLOCATE(pdpsrf(ngrid))
@@ -446 +448 @@
   !! parameters are defined in the module_model_constants.F WRF routine
   PRINT *,'** Mars ** INITIALIZE ARRAYS FOR PHYSICS'
+  !! need to get initial time first
+  CALL update_inputs_physiq_time(&
+            JULYR,JULDAY,GMT,&
+            elaps,&
+            lct_input,lon_input,ls_input,&
+            ptime,pday,MY)
+  !!
   CALL phys_state_var_init(ngrid,nlayer,nq,&
+                           floor(pday),ptime,&
                            wdaysec,ptimestep,1/reradius,g,r_d,cp)
   !! Fill planetary parameters in modules
@@ -653 +663 @@
              firstcall,lastcall,pday,ptime,ptimestep, &
              pplev,pplay,pphi,pu,pv,pt,pq,pw,         &
-             pdu,pdv,pdt,pdq,pdpsrf,tracerdyn)
+             pdu,pdv,pdt,pdq,pdpsrf)
 #endif