Changeset 2333

Timestamp:

May 29, 2020, 9:46:23 AM (5 years ago)

Author:

adelavois

Message:

More outputs with xios from physics
physics is now able to transmit soil fields to xios

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/physiq_mod.F (modified) (2 diffs)
• libf/phymars/xios_output_mod.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/README

@@ -3074 +3074 @@
 - libf/phymars/dyn1d:
   testphys1d.F: idem (the reading interface for traceur.def has been completed to fill the variables nqperes and nqfils).
+
+== 29/05/2020 == AD
+More outputs with xios from physics
+physics is now able to transmit soil fields to xios

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

@@ -697 +697 @@
         write(*,*) "physiq firstcall: call initialize_xios_output"
         call initialize_xios_output(pday,ptime,ptimestep,daysec,
-     &                              presnivs,pseudoalt)
+     &                              presnivs,pseudoalt,mlayer)
 #endif
       ENDIF        !  (end of "if firstcall")
@@ -3745 +3745 @@
       CALL send_xios_field("ps",ps)
       CALL send_xios_field("area",cell_area)
+
+      CALL send_xios_field("ISR",fluxtop_sw_tot)
+      CALL send_xios_field("OLR",fluxtop_lw)
+
+      CALL send_xios_field("tsurf",tsurf)
+      CALL send_xios_field("inertiedat",inertiedat)
+      CALL send_xios_field("tsoil",tsoil)
+      CALL send_xios_field("co2ice",co2ice)
       
       CALL send_xios_field("temperature",zt)
       CALL send_xios_field("u",zu)
       CALL send_xios_field("v",zv)
+      !CREATE IF CO2CYCLE
+      !if (tracer.and.(igcm_co2.ne.0)) then
+      !   CALL send_xios_field("co2",zq(1,1,igcm_co2))
+      !endif
+
 
       if (lastcall.and.is_omp_master) then

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

@@ -18 +18 @@
 
   SUBROUTINE initialize_xios_output(day,timeofday,dtphys,daysec,&
-                                    presnivs,pseudoalt)
+                                    presnivs,pseudoalt,mlayer)
   USE mod_phys_lmdz_para, ONLY: jj_nb, jj_begin, jj_end, ii_begin, ii_end, &
                                 mpi_size, mpi_rank, klon_mpi, &
@@ -39 +39 @@
   REAL,INTENT(IN) :: presnivs(:) ! vertical grid approximate pressure (Pa)
   REAL,INTENT(IN) :: pseudoalt(:) ! vertical grid approximate altitude (km)
+  REAL,INTENT(IN) :: mlayer(:) ! soil layer depth at intermediate level (m)
   
   
@@ -56 +57 @@
     CALL xios_set_axis_attr("altitude", n_glo=size(pseudoalt), value=pseudoalt,&
                             unit="km",positive="up")
+    IF (prt_level>=10) WRITE(lunout,*) "initialize_xios_output: call xios_set_axis_attr for soil"
+    CALL xios_set_axis_attr("soil_layers", n_glo=size(mlayer), value=mlayer,&
+                            unit="m",positive="down")
     
     ! 2. Declare horizontal domain