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

Timestamp:

Aug 23, 2016, 10:42:41 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 4: call to call_physiq in module_lmd_driver
(callphysiq_mod is in WRFV2/phys)

module_lmd_driver.F.new is now supposedly independent
of the considered planet

to build an interface with new planetary physics
consider adapting

• update_inputs_physiq_mod.F
• update_outputs_physiq_mod.F
• iniphysiq_mod.F
• callphysiq_mod.F

for the last two, see examples for a given planet
in LMDZ._yourplanet_/libf/dynphy_lonlat/phy_yourplanet_

NB: checked compatibility of results with previous commit

File:

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

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

@@ -13 +13 @@
 !           - unified module_lmd_driver: old, new phys and LES - February 2011
 !           - a new paradigm to prepare nested simulations - April 2014
+!           - adapted to new interface philosophy & other planets - August 2016
 !*******************************************************************************
 MODULE module_lmd_driver
@@ -63 +64 @@
    USE comm_wrf !! to get fields to be written from physiq
    USE iniphysiq_mod !! to get iniphysiq subroutine
-   USE physiq_mod, only : physiq
+   USE callphysiq_mod !! to call the LMD physics
 #endif
 
@@ -132 +133 @@
    
 
-   ! *** for Mars Mode 20 ***
+   ! *** for tracer Mode 20 ***
    REAL ::    tau_decay
    ! *** ----------------------- ***
@@ -143 +144 @@
    LOGICAL :: firstcall,lastcall
    ! ---------- 
-   REAL,DIMENSION(:,:),ALLOCATABLE :: pplev,pplay,pphi,pu,pv,pt,pw
+   REAL,DIMENSION(:,:),ALLOCATABLE :: pplev,pplay,pphi,pu,pv,pt,flxw
    REAL,DIMENSION(:,:,:),ALLOCATABLE :: pq  
 
@@ -191 +192 @@
 !!!IDEALIZED IDEALIZED
 
-
-!REAL, DIMENSION(:,:,:), ALLOCATABLE, SAVE :: &
-!     UMAX, UMIN, VMAX, VMIN, WMAX, WMIN, TMAX, TMIN
-
-
+   !! arguments to physiq
+   CHARACTER(len=20),ALLOCATABLE :: tname(:) ! tracer names
+   REAL,ALLOCATABLE :: zpk_omp(:,:)
+   REAL,ALLOCATABLE :: zphis_omp(:) ! surface geopotential
+   REAL,ALLOCATABLE :: presnivs_omp(:) ! approximate pressure of atm. layers  
+   REAL,ALLOCATABLE :: zrfi_omp(:,:) ! relative wind vorticity, in s-1
+
+   !! this is temporary and supposed to be set by namelist
+   character(len=10),save :: planet_type ! planet type ('earth','mars',...)
+   planet_type = "mars"
 
 !==================================================================
@@ -225 +231 @@
 
 
-print *,'** Mars ** DOMAIN',id
+print *,'** ',planet_type,' ** DOMAIN',id
 
 !-------------------------!
@@ -356 +362 @@
 !!!***********!!
 IF (JULYR .eq. 9999) THEN
-  PRINT *,'** Mars ** IDEALIZED SIMULATION'
-  PRINT *,'** Mars ** BEWARE: input_coord must be here'
+  PRINT *,'** ',planet_type,'** IDEALIZED SIMULATION'
+  PRINT *,'** ',planet_type,'** BEWARE: input_coord must be here'
   open(unit=14,file='input_coord',form='formatted',status='old')
   rewind(14)
@@ -370 +376 @@
 ! ALLOCATE !
 !----------!
+IF (.not.ALLOCATED(tname)) ALLOCATE(tname(nq)) 
 !-------------------------------------------------------------------------------!
 ! outputs:                                                                      !       
@@ -387 +394 @@
 !!!
 IF (test.NE.0) THEN
-print *,'** Mars ** NO CALL FOR PHYSICS, go to next step...',test
+print *,'** ',planet_type,'** NO CALL FOR PHYSICS, go to next step...',test
 #ifdef SPECIAL_NEST_SAVE
 pdpsrf(:)=dp_save(:,id)
@@ -418 +425 @@
 ALLOCATE(pv(ngrid,nlayer))       !!!!!
 ALLOCATE(pt(ngrid,nlayer))       !!!!!
-ALLOCATE(pw(ngrid,nlayer))       !!!!!
+ALLOCATE(flxw(ngrid,nlayer))       !!!!!
 ALLOCATE(pq(ngrid,nlayer,nq))    !!!!!
+ALLOCATE(zpk_omp(ngrid,nlayer))
+ALLOCATE(zphis_omp(ngrid))
+ALLOCATE(presnivs_omp(nlayer))
+ALLOCATE(zrfi_omp(ngrid,nlayer))
 ! interm
 ALLOCATE(dz8w_prof(nlayer))
@@ -442 +453 @@
 allocation_firstcall: IF (firstcall .EQV. .true.) THEN
   !! tracers' name
-  PRINT *,'** Mars ** TRACERS NAMES'
-  CALL update_inputs_physiq_tracers(nq,MARS_MODE)
+  PRINT *,'** ',planet_type,'** TRACERS NAMES'
+  CALL update_inputs_physiq_tracers(nq,MARS_MODE,tname)
   !! PHYSICS VARIABLES (cf. iniphysiq in LMD GCM)
   !! parameters are defined in the module_model_constants.F WRF routine
-  PRINT *,'** Mars ** INITIALIZE ARRAYS FOR PHYSICS'
+  PRINT *,'** ',planet_type,'** INITIALIZE ARRAYS FOR PHYSICS'
   !! need to get initial time first
   CALL update_inputs_physiq_time(&
@@ -538 +549 @@
 pu(subs,:) = u_prof(:)
 pv(subs,:) = v_prof(:)
-pw(subs,:) = 0   !! NB: not used in the physics, only diagnostic...
+flxw(subs,:) = 0   !! NB: not used in the physics, only diagnostic...
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !! for IDEALIZED CASES ONLY
@@ -584 +595 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 pass_interface: IF ( (firstcall .EQV. .true.) .or. (max_dom .GT. 1) ) THEN
-PRINT *,'** Mars ** PASS INTERFACE. EITHER FIRSTCALL or NESTED SIMULATION.'
+PRINT *,'** ',planet_type,'** PASS INTERFACE. EITHER FIRSTCALL or NESTED SIMULATION.'
 !!*******************************************!!
 !!*******************************************!!
@@ -649 +660 @@
 pdq(:,:,:)=0.
 #ifndef NOPHYS
-print *, '** Mars ** CALL TO LMD PHYSICS'
+print *, '** ',planet_type,'** CALL TO LMD PHYSICS'
 !!!
 CALL update_inputs_physiq_time(&
@@ -657 +668 @@
             ptime,pday,MY)
 !!!
-CALL physiq (ngrid,nlayer,nq,                         &
-             firstcall,lastcall,pday,ptime,ptimestep, &
-             pplev,pplay,pphi,pu,pv,pt,pq,pw,         &
-             pdu,pdv,pdt,pdq,pdpsrf)
+CALL call_physiq(planet_type,ngrid,nlayer,nq,tname, &
+                       firstcall,lastcall,          &
+                       pday,ptime,ptimestep,        &
+                       pplev,pplay,                 &
+                       zpk_omp,pphi,zphis_omp,      &
+                       presnivs_omp,                &
+                       pu,pv,zrfi_omp,pt,pq,        &
+                       flxw,                        &
+                       pdu,pdv,pdt,pdq,pdpsrf)
+!!!
 #endif
 
@@ -670 +687 @@
 #endif
 
-print *, '** Mars ** CALL TO LMD PHYSICS DONE'
+print *, '** ',planet_type,'** CALL TO LMD PHYSICS DONE'
 DEALLOCATE(pplev)
 DEALLOCATE(pplay)
@@ -677 +694 @@
 DEALLOCATE(pv)
 DEALLOCATE(pt)
-DEALLOCATE(pw)
+DEALLOCATE(flxw)
 DEALLOCATE(pq)
+DEALLOCATE(zpk_omp)
+DEALLOCATE(zphis_omp)
+DEALLOCATE(presnivs_omp)
+DEALLOCATE(zrfi_omp)
+
 
 !---------------------------------------------------------------------------------!
@@ -796 +818 @@
       SCALAR(i,kps:kpe,j,:)=0.
     CASE(20)
-      !! Mars mode 20 : add a passive tracer with radioactive-like decay
+      !! tracer mode 20 : add a passive tracer with radioactive-like decay
       IF ( (i == ips) .AND. (j == jps) )   print *, 'RADIOACTIVE-LIKE TRACER WITH SOURCE AT SURFACE LAYER.'
       tau_decay=60.*10. !! why not make it a namelist argument?
@@ -817 +839 @@
 !! END !!
 !!*****!!
-print *,'** Mars ** END LMD PHYSICS'
+print *,'** ',planet_type,'** END LMD PHYSICS'
 !----------------------------------------------------------------!
 ! use debug (see solve_em) if you wanna display some message ... !