Changeset 1622 for trunk/LMDZ.GENERIC

Timestamp:

Oct 31, 2016, 4:09:55 PM (8 years ago)

Author:

emillour

Message:

Generic GCM:
Introducing a first prototype of XIOS outputs: requires compiling with
"-io xios" flag and having an XIOS library at hand. Corresponding
companion xml files iodef.xml and context_lmdz_physics.xml in the "deftank"
subdirectory. Still some work to do concerning the XIOS calendar (definition
of an initial date and month number and lenght).
EM

Location:

trunk/LMDZ.GENERIC

Files:

• README (modified) (1 diff)
• deftank/context_lmdz_physics.xml (added)
• deftank/iodef.xml (added)
• libf/phystd/physiq_mod.F90 (modified) (6 diffs)
• libf/phystd/xios_output_mod.F90 (added)

trunk/LMDZ.GENERIC/README

@@ -1275 +1275 @@
 - added scalheight (atmospheric scale height) in comvert_mod.
 
+== 31/10/2016 == EM
+Introducing a first prototype of XIOS outputs: requires compiling with
+"-io xios" flag and having an XIOS library at hand. Corresponding
+companion xml files iodef.xml and context_lmdz_physics.xml in the "deftank"
+subdirectory. Still some work to do concerning the XIOS calendar (definition
+of an initial date and month number and lenght).
+

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -46 +46 @@
       use time_phylmdz_mod, only: daysec
       use callkeys_mod
+      use vertical_layers_mod, only: presnivs, pseudoalt
+      use xios_output_mod, only: initialize_xios_output, &
+                                 update_xios_timestep, &
+                                 send_xios_field
       implicit none
 
@@ -198 +202 @@
 ! ----------------------
 
-      integer day_ini                                              ! Initial date of the run (sol since Ls=0).
-      integer icount                                               ! Counter of calls to physiq during the run.
+      integer,save :: day_ini                                      ! Initial date of the run (sol since Ls=0).
+      integer,save :: icount                                       ! Counter of calls to physiq during the run.
+!$OMP THREADPRIVATE(day_ini,icount)
+
       real, dimension(:),allocatable,save ::  tsurf                ! Surface temperature (K).
       real, dimension(:,:),allocatable,save ::  tsoil              ! Sub-surface temperatures (K).
@@ -225 +231 @@
 !$OMP THREADPRIVATE(emis,dtrad,fluxrad_sky,fluxrad,capcal,fluxgrd,qsurf,q2)
 
-      save day_ini, icount
-      
-!$OMP THREADPRIVATE(day_ini,icount)
 
 ! Local variables :
@@ -706 +709 @@
          endif
          
-               
+         ! XIOS outputs
+#ifdef CPP_XIOS
+
+         write(*,*) "physiq: call initialize_xios_output"
+         call initialize_xios_output(pday,ptime,ptimestep,daysec, &
+                                     presnivs,pseudoalt)
+#endif
       endif ! end of 'firstcall'
 
@@ -712 +721 @@
 ! I.2   Initializations done at every physical timestep:
 ! ------------------------------------------------------
+
+#ifdef CPP_XIOS      
+      ! update XIOS time/calendar
+      call update_xios_timestep
+#endif      
 
       ! Initialize various variables
@@ -2132 +2146 @@
       endif
 
+! 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("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)
+      
+#endif
+
       icount=icount+1