Changeset 2665 for LMDZ5/trunk/libf/dynphy_lonlat/phylmd/limit_netcdf.F90

Timestamp:

Oct 12, 2016, 2:53:20 PM (8 years ago)

Author:

dcugnet

Message:

• A (re)startphy.nc file (standard name: "startphy0.nc") can be read by ce0l to get land mask, so mask can be defined (in decreasing priority order) from: 1) "o2a.nc file" if this file is found 2) "startphy0.nc" if this file is found 3) "Relief.nc" otherwise
• Sub-cell scales parameters for orographic gravity waves can be read from file "oro_params.nc" if the configuration key "read_orop" is TRUE. The effect is to bypass the "grid_noro" routine in ce0l, so that any pre-defined mask (from o2a.nc or startphy0.nc) is then overwritten.
• The gcm stops if the "limit.nc" records number differs from the current year number of days. A warning is issued in case the gcm calendar does not match the time axis attribute "calendar" (if available) from the "limit.nc" file. This attribute is now added to the "limit.nc" time axis.
• Few simplifications in grid_noro
• Few parameters changes in acama_gwd and flott_gwd.
• Variable d_u can be saved in the outputs.

File:

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

LMDZ5/trunk/libf/dynphy_lonlat/phylmd/limit_netcdf.F90

@@ -71 +71 @@
   USE netcdf95,           ONLY: nf95_def_var, nf95_put_att, nf95_put_var
   USE comconst_mod, ONLY: pi
+  USE phys_cal_mod, ONLY: calend
   IMPLICIT NONE
 !-------------------------------------------------------------------------------
@@ -244 +245 @@
   !--- Attributes creation
   CALL ncerr(NF90_PUT_ATT(nid,id_tim, "title","Jour dans l annee"),fnam)
+  CALL ncerr(NF90_PUT_ATT(nid,id_tim, "calendar",calend),fnam)
   CALL ncerr(NF90_PUT_ATT(nid,id_FOCE,"title","Fraction ocean"),fnam)
   CALL ncerr(NF90_PUT_ATT(nid,id_FSIC,"title","Fraction glace de mer"),fnam)
@@ -405 +407 @@
       CASE('SIC', 'SST'); cal_in='gregorian'
     END SELECT
-  CALL msg(5,'WARNING: missing "calendar" attribute for "time" in '&
+    CALL msg(5,'WARNING: missing "calendar" attribute for "time" in '&
      &//TRIM(fnam)//'. Choosing default value.')
   END IF
+  CALL strclean(cal_in)                     !--- REMOVE (WEIRD) NULL CHARACTERS
   CALL msg(5,'var, calendar, dim: '//TRIM(dnam)//' '//TRIM(cal_in), lmdep)
   
@@ -477 +480 @@
   fnam_p=fnam(1:idx)//'_p.nc'
   IF(NF90_OPEN(fnam_p,NF90_NOWRITE,ncid)==NF90_NOERR) THEN
-    CALL msg(0,'Reading previous year file ("'//TRIM(fnam_p)//'") first record for '//TRIM(title))
+    CALL msg(0,'Reading next year file ("'//TRIM(fnam_p)//'") first record for '//TRIM(title))
     CALL ncerr(NF90_INQ_VARID(ncid, varname, varid),fnam_p)
     CALL ncerr(NF90_GET_VAR(ncid,varid,champ,[1,1,1],[imdep,jmdep,1]),fnam_p)
@@ -767 +770 @@
 !-------------------------------------------------------------------------------
 
+
+!-------------------------------------------------------------------------------
+!
+SUBROUTINE strclean(s)
+!
+!-------------------------------------------------------------------------------
+  IMPLICIT NONE
+!-------------------------------------------------------------------------------
+! Purpose: Remove tail null characters from the input string.
+!-------------------------------------------------------------------------------
+! Parameters:
+  CHARACTER(LEN=*), INTENT(INOUT) :: s
+!-------------------------------------------------------------------------------
+! Local variable:
+  INTEGER :: k
+!-------------------------------------------------------------------------------
+  k=LEN_TRIM(s); DO WHILE(ICHAR(s(k:k))==0); s(k:k)=' '; k=LEN_TRIM(s); END DO
+
+END SUBROUTINE strclean
+!
+!-------------------------------------------------------------------------------
+
 #endif
 ! of #ifndef CPP_1D