Changeset 2223 for trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

Timestamp:

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

Author:

emillour

Message:

Mars GCM:
First implementation of XIOS in the physics
EM

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (5 diffs)

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