Ignore:
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/libf/phylmd
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

    r3048 r3059  
    1616!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    1717!  CHARACTER(len=20), dimension(nfiles) :: TEF = type_ecri_files
     18
     19!!! saving lon and lat as variables for CMIP6 DataRequest
     20  TYPE(ctrl_out), SAVE :: o_longitude = ctrl_out((/ 1, 1, 1, 1, 1, 1, 11, 11, 11, 11/), &
     21    'io_lon', '', '', (/ ('once', i=1, 10) /))
     22  TYPE(ctrl_out), SAVE :: o_latitude = ctrl_out((/ 1, 1, 1, 1, 1, 1, 11, 11, 11, 11/), &
     23    'io_lat', '', '', (/ ('once', i=1, 10) /))
    1824
    1925!!! Comosantes de la coordonnee sigma-hybride
  • LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

    r3048 r3059  
    3131    USE vertical_layers_mod, ONLY : ap, bp, aps, bps
    3232    USE phys_output_ctrlout_mod, ONLY: o_phis, o_aire, is_ter, is_lic, is_oce, &
     33         o_longitude, o_latitude, &
    3334         o_Ahyb, o_Bhyb,o_Ahyb_inter, o_Bhyb_inter, &
    3435         is_ave, is_sic, o_contfracATM, o_contfracOR, &
     
    336337    USE indice_sol_mod, ONLY: nbsrf
    337338    USE infotrac_phy, ONLY: nqtot, nqo, type_trac, tname, niadv
    338     USE geometry_mod, ONLY: cell_area
     339    USE geometry_mod, ONLY: cell_area, latitude_deg, longitude_deg
    339340    USE surface_data, ONLY: type_ocean, version_ocean, ok_veget, ok_snow
    340341    USE aero_mod, ONLY: naero_tot, id_STRAT_phy
     
    547548       CALL histwrite_phy(o_Ahyb_inter, Ahyb_bounds)
    548549       CALL histwrite_phy(o_Bhyb_inter, Bhyb_bounds)
     550       CALL histwrite_phy(o_longitude, longitude_deg)
     551       CALL histwrite_phy(o_latitude, latitude_deg)
    549552!
    550553#ifdef CPP_RRTM
Note: See TracChangeset for help on using the changeset viewer.