Changeset 2154 for LMDZ5/trunk/libf/phylmd/limit_netcdf.F90

Timestamp:

Nov 25, 2014, 10:49:50 AM (10 years ago)

Author:

acozic

Message:

To create limit.nc add the possibility to read sst and sic from erai files
(wait as sst.nc and ci.nc by the model)

File:

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

LMDZ5/trunk/libf/phylmd/limit_netcdf.F90

@@ -55 +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 ------------------------------------------
@@ -119 +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
@@ -180 +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