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

Timestamp:

May 24, 2017, 4:02:35 PM (7 years ago)

Author:

dcugnet

Message:

Change for limit files (ce0l:limit_netcdf.F90): daily records are interpolated
at 12h (no longer 0h), so that the 0h-24h mean value is used the whole day.

File:

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

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

@@ -565 +565 @@
     l=1
     DO k=1, ndays
-      timeday = REAL((k-1)*ndays_in)/ndays
+      timeday = (REAL(k)-0.5)*REAL(ndays_in)/ndays
       IF(timeyear(l+1)<timeday) l=l+1
       al=(timeday-timeyear(l))/(timeyear(l+1)-timeyear(l))
@@ -583 +583 @@
               vc_beg=0., vc_end=0.)
          CALL pchfe_95(timeyear, champtime(i, j, :), yder, skip, &
-              arth(0., real(ndays_in) / ndays, ndays), champan(i, j, :), ierr)
+              arth(0.5, real(ndays_in) / ndays, ndays), champan(i, j, :), ierr)
          if (ierr < 0) stop 1
          n_extrap = n_extrap + ierr