Changeset 2223

Timestamp:

Jan 23, 2020, 10:46:58 AM (5 years ago)

Author:

emillour

Message:

Mars GCM:
First implementation of XIOS in the physics
EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• deftank/README (modified) (1 diff)
• deftank/context_lmdz_physics.xml (added)
• deftank/iodef.xml (added)
• libf/phymars/phys_state_var_init_mod.F90 (modified) (2 diffs)
• libf/phymars/physiq_mod.F (modified) (5 diffs)
• libf/phymars/xios_output_mod.F90 (added)

trunk/LMDZ.MARS/README

@@ -2833 +2833 @@
 Update the nonorographic gravity waves drag parametrization (from the r3599 of Earth s model LMDZ6): 1) add east_ and west_gwdstress variables, 2) delete aleas function: random waves are producted by the MOD function, 3) reproductibility concerning the number of procs is now validated, 4) changing name of some variables like RUW, RVW, ZOP, ZOM,..., 5) tendency of winds due to GW drag are now module variables and written in restart files. 
 
+== 23/01/2020 == EM
+First implementation of XIOS in the physics

trunk/LMDZ.MARS/deftank/README

@@ -20 +20 @@
 - *.MCD5 files correspond to what was used to generate MCDv5.x simulations 
 
+>> XML FILES
+These are required when running with the XIOS library
+- iodef.xml : general settings
+- context_lmdz_physics.xml : definition of variables, grids, output files, etc.
+
 >> SCRIPTS
 - "run0" is a script to run (once or nummax times) chained simulations

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

@@ -34 +34 @@
 !=======================================================================
 
+      use init_print_control_mod, only: init_print_control
       use slope_mod, only: ini_slope_mod,end_slope_mod
       use comsaison_h, only: ini_comsaison_h,end_comsaison_h
@@ -69 +70 @@
       call ini_tracer_mod(nq,tname)
 
+
+      ! initialize "print_control" constants/flags ("prt_level","lunout", etc.)
+      call init_print_control
+      
       ! set parameters in comcstfi_h
       pi=2.*asin(1.) 

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

@@ -71 +71 @@
       use comm_wrf
 #else
-      use planetwide_mod
+      USE planetwide_mod, ONLY: planetwide_maxval, planetwide_minval,
+     &                          planetwide_sumval
       use phyredem, only: physdem0, physdem1
       use phyetat0_mod, only: phyetat0
       use eofdump_mod, only: eofdump
-      USE vertical_layers_mod, ONLY: ap,bp,aps,bps
-     
+      USE vertical_layers_mod, ONLY: ap,bp,aps,bps,presnivs,pseudoalt
+      USE mod_phys_lmdz_omp_data, ONLY: is_omp_master
 #endif
 
-
+#ifdef CPP_XIOS      
+      use xios_output_mod, only: initialize_xios_output,
+     &                           update_xios_timestep,
+     &                           send_xios_field
+      use wxios, only: wxios_context_init, xios_context_finalize
+#endif
 
       IMPLICIT NONE
@@ -491 +497 @@
 #endif
 
+#ifdef CPP_XIOS
+         ! Initialize XIOS context
+         write(*,*) "physiq: call wxios_context_init"
+         CALL wxios_context_init
+#endif
+
 c        read startfi 
 c        ~~~~~~~~~~~~
@@ -630 +642 @@
 #endif
                  
+#ifdef CPP_XIOS
+        ! XIOS outputs
+        write(*,*) "physiq firstcall: call initialize_xios_output"
+        call initialize_xios_output(pday,ptime,ptimestep,daysec,
+     &                              presnivs,pseudoalt)
+#endif
       ENDIF        !  (end of "if firstcall")
 
@@ -637 +655 @@
 c ---------------------------------------------------
 c
+
+#ifdef CPP_XIOS      
+      ! update XIOS time/calendar
+      call update_xios_timestep
+#endif      
 
 c     Initialize various variables
@@ -3444 +3467 @@
       END IF       ! if(ngrid.ne.1)
 
+! XIOS outputs
+#ifdef CPP_XIOS      
+      ! Send fields to XIOS: (NB these fields must also be defined as
+      ! <field id="..." /> in context_lmdz_physics.xml to be correctly used)
+      CALL send_xios_field("ls",zls*180./pi)
+      
+      CALL send_xios_field("ps",ps)
+      CALL send_xios_field("area",cell_area)
+      
+      CALL send_xios_field("temperature",zt)
+      CALL send_xios_field("u",zu)
+      CALL send_xios_field("v",zv)
+
+      if (lastcall.and.is_omp_master) then
+        write(*,*) "physiq lastcall: call xios_context_finalize"
+        call xios_context_finalize
+      endif
+#endif
+
+
       icount=icount+1