Changeset 2128

Timestamp:

Oct 20, 2014, 4:52:34 PM (10 years ago)

Author:

lguez

Message:

Added longitude and latitude to file "limit.nc", so that we can more easily
plot fields from this file.

File:

• LMDZ5/trunk/libf/phylmd/limit_netcdf.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/limit_netcdf.F90

@@ -25 +25 @@
   USE dimphy
   USE ioipsl,             ONLY : ioget_year_len
-  USE phys_state_var_mod, ONLY : pctsrf
+  USE phys_state_var_mod, ONLY : pctsrf, rlon, rlat
   USE netcdf,             ONLY : NF90_OPEN,    NF90_CREATE,  NF90_CLOSE,       &
                    NF90_DEF_DIM, NF90_DEF_VAR, NF90_PUT_VAR, NF90_PUT_ATT,     &
@@ -31 +31 @@
                    NF90_CLOBBER, NF90_ENDDEF,  NF90_UNLIMITED, NF90_FLOAT
   USE inter_barxy_m, only: inter_barxy
+  use netcdf95, only: nf95_def_var, nf95_put_att, nf95_put_var
 #endif
   IMPLICIT NONE
@@ -75 +76 @@
   INTEGER, DIMENSION(2) :: dims
   INTEGER :: id_tim,  id_SST,  id_BILS, id_RUG, id_ALB
-  INTEGER :: id_FOCE, id_FSIC, id_FTER, id_FLIC
+  INTEGER :: id_FOCE, id_FSIC, id_FTER, id_FLIC, varid_longitude, varid_latitude
   INTEGER :: NF90_FORMAT
   INTEGER :: ndays                   !--- Depending on the output calendar
@@ -220 +221 @@
   ierr=NF90_DEF_VAR(nid,"ALB",  NF90_FORMAT,dims,id_ALB)
   ierr=NF90_DEF_VAR(nid,"RUG",  NF90_FORMAT,dims,id_RUG)
+  call nf95_def_var(nid, "longitude", NF90_FLOAT, ndim, varid_longitude)
+  call nf95_def_var(nid, "latitude", NF90_FLOAT, ndim, varid_latitude)
 
   !--- Attributes creation
@@ -232 +235 @@
   ierr=NF90_PUT_ATT(nid,id_RUG, "title","Rugosite")
 
+  call nf95_put_att(nid, varid_longitude, "standard_name", "longitude")
+  call nf95_put_att(nid, varid_longitude, "units", "degrees_east")
+
+  call nf95_put_att(nid, varid_latitude, "standard_name", "latitude")
+  call nf95_put_att(nid, varid_latitude, "units", "degrees_north")
+
   ierr=NF90_ENDDEF(nid)
 
@@ -244 +253 @@
   ierr=NF90_PUT_VAR(nid,id_ALB ,phy_alb(:,:),(/1,1/),(/klon,ndays/))
   ierr=NF90_PUT_VAR(nid,id_RUG ,phy_rug(:,:),(/1,1/),(/klon,ndays/))
+  call nf95_put_var(nid, varid_longitude, rlon)
+  call nf95_put_var(nid, varid_latitude, rlat)
 
   ierr=NF90_CLOSE(nid)