Changeset 2160 for LMDZ5/branches/testing/libf/phylmd/limit_netcdf.F90

Timestamp:

Nov 28, 2014, 4:36:29 PM (10 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes -r2070:2158 into testing branch. Compilation problems introduced by revision r2155 have been corrected by hand

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/limit_netcdf.F90 (modified) (9 diffs)

LMDZ5/branches/testing/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
@@ -54 +55 @@
 !--- INPUT NETCDF FILES NAMES --------------------------------------------------
   CHARACTER(LEN=25) :: icefile, sstfile, dumstr
-  CHARACTER(LEN=25), PARAMETER :: famipsst='amipbc_sst_1x1.nc        ',        &
-                                  famipsic='amipbc_sic_1x1.nc        ',        &
-                                  fcpldsst='cpl_atm_sst.nc           ',        &
-                                  fcpldsic='cpl_atm_sic.nc           ',        &
-                                  fhistsst='histmth_sst.nc           ',        &
-                                  fhistsic='histmth_sic.nc           ',        &
-                                  frugo   ='Rugos.nc                 ',        &
-                                  falbe   ='Albedo.nc                '
+  CHARACTER(LEN=25), PARAMETER :: famipsst ='amipbc_sst_1x1.nc        ',        &
+                                  famipsic ='amipbc_sic_1x1.nc        ',        &
+                                  fcpldsst ='cpl_atm_sst.nc           ',        &
+                                  fcpldsic ='cpl_atm_sic.nc           ',        &
+                                  fhistsst ='histmth_sst.nc           ',        &
+                                  fhistsic ='histmth_sic.nc           ',        &
+                                  frugo    ='Rugos.nc                 ',        &
+                                  falbe    ='Albedo.nc                ',        &
+                                  feraisstk='sstk.nc                  ',        &
+                                  feraici  ='ci.nc                    '
   CHARACTER(LEN=10) :: varname
 !--- OUTPUT VARIABLES FOR NETCDF FILE ------------------------------------------
@@ -75 +78 @@
   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
@@ -118 +121 @@
      icefile=TRIM(fhistsic)
      varname='pourc_sic'
+  ELSE IF ( NF90_OPEN(TRIM(feraici),NF90_NOWRITE,nid)==NF90_NOERR ) THEN
+     icefile=TRIM(feraici)
+     varname='ci'
   ELSE
      WRITE(lunout,*) 'ERROR! No sea-ice input file was found.'
-     WRITE(lunout,*) 'One of following files must be availible : ',trim(famipsic),', ',trim(fcpldsic),', ',trim(fhistsic)
+     WRITE(lunout,*) 'One of following files must be availible : ',trim(famipsic),', ',trim(fcpldsic),', ',trim(fhistsic), trim(feraici)
      CALL abort_gcm('limit_netcdf','No sea-ice file was found',1)
   END IF
@@ -179 +185 @@
      sstfile=TRIM(fhistsst)
      varname='tsol_oce'
+  ELSE IF ( NF90_OPEN(TRIM(feraisstk),NF90_NOWRITE,nid)==NF90_NOERR ) THEN
+     sstfile=TRIM(feraisstk)
+     varname='sstk'
   ELSE
      WRITE(lunout,*) 'ERROR! No sst input file was found.'
-     WRITE(lunout,*) 'One of following files must be availible : ',trim(famipsst),trim(fcpldsst),trim(fhistsst)
+     WRITE(lunout,*) 'One of following files must be availible : ',trim(famipsst),trim(fcpldsst),trim(fhistsst),trim(feraisstk)
      CALL abort_gcm('limit_netcdf','No sst file was found',1)
   END IF
@@ -220 +229 @@
   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 +243 @@
   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 +261 @@
   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)