Changeset 3059 for LMDZ6

Timestamp:

Nov 8, 2017, 9:31:44 AM (7 years ago)

Author:

Laurent Fairhead

Message:

Adding longitude and latitude as output variables as per CMIP6 Data Request
LF

Location:

LMDZ6/trunk

Files:

• DefLists/field_def_lmdz.xml (modified) (1 diff)
• DefLists/file_def_histmth_lmdz.xml (modified) (1 diff)
• libf/phylmd/phys_output_ctrlout_mod.F90 (modified) (1 diff)
• libf/phylmd/phys_output_write_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/DefLists/field_def_lmdz.xml

@@ -6 +6 @@
 <!--        <field id="Alt"  long_name="Pseudo-altitude (for an 8km scale height)"    unit="km" />  -->
 <!--    </field_group>  -->
+
+    <field_group id="coordinates" grid_ref="grid_glo">
+        <field id="io_lon"  />
+        <field id="io_lat"  />
+    </field_group>
 
     <field_group id="coord_hyb">

LMDZ6/trunk/DefLists/file_def_histmth_lmdz.xml

@@ -11 +11 @@
 
             <field_group operation="once">
+              <field field_ref="io_lon" level="1" />
+              <field field_ref="io_lat" level="1" />
               <field field_ref="Ahyb" level="1" />
               <field field_ref="Ahyb_inter" level="1" />

LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -16 +16 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !  CHARACTER(len=20), dimension(nfiles) :: TEF = type_ecri_files
+
+!!! saving lon and lat as variables for CMIP6 DataRequest
+  TYPE(ctrl_out), SAVE :: o_longitude = ctrl_out((/ 1, 1, 1, 1, 1, 1, 11, 11, 11, 11/), &
+    'io_lon', '', '', (/ ('once', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_latitude = ctrl_out((/ 1, 1, 1, 1, 1, 1, 11, 11, 11, 11/), &
+    'io_lat', '', '', (/ ('once', i=1, 10) /))
 
 !!! Comosantes de la coordonnee sigma-hybride

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -31 +31 @@
     USE vertical_layers_mod, ONLY : ap, bp, aps, bps
     USE phys_output_ctrlout_mod, ONLY: o_phis, o_aire, is_ter, is_lic, is_oce, &
+         o_longitude, o_latitude, &
          o_Ahyb, o_Bhyb,o_Ahyb_inter, o_Bhyb_inter, & 
          is_ave, is_sic, o_contfracATM, o_contfracOR, &
@@ -336 +337 @@
     USE indice_sol_mod, ONLY: nbsrf
     USE infotrac_phy, ONLY: nqtot, nqo, type_trac, tname, niadv
-    USE geometry_mod, ONLY: cell_area
+    USE geometry_mod, ONLY: cell_area, latitude_deg, longitude_deg
     USE surface_data, ONLY: type_ocean, version_ocean, ok_veget, ok_snow
     USE aero_mod, ONLY: naero_tot, id_STRAT_phy
@@ -547 +548 @@
        CALL histwrite_phy(o_Ahyb_inter, Ahyb_bounds)
        CALL histwrite_phy(o_Bhyb_inter, Bhyb_bounds)
+       CALL histwrite_phy(o_longitude, longitude_deg)
+       CALL histwrite_phy(o_latitude, latitude_deg)
 !
 #ifdef CPP_RRTM