Changeset 1220 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/regr_lat_time_climoz_m.F90

Timestamp:

Aug 5, 2009, 4:38:34 PM (15 years ago)

Author:

lguez

Message:

-- Replaced "integer*4" declarations by "integer", "real*8" by

"real(kind=8)" and "real*4" by "real". Note that these are the only
modifications in the files "radiation_AR4.F" and "sw_aeroAR4.F90".

-- Corrected the kind of arguments to "max" and "min".

-- Replaced "nH" edit descriptors, which is a deleted feature in

Fortran 95, by character strings.

-- "regr_lat_time_climoz" now allows the pressure coordinate in the

input file to be in descending order.

-- Replaced call to not standard function "float" by call to intrinsic

function "real".

-- Included file "radepsi.h" in "physiq" was not used. Removed it.

The following set of modifications is related to the management of time.

-- In "gcm", "leapfrog" and "sortvarc0", "day_ini" was defined as 1

plus number of days between the reference date "(annee_ref,
day_ref)" and the first day of the current simulation. Changed
definition: "(annee_ref, day_ini)" is the first day of the current
simulation. There is an accompanying modification for "day_end".

-- Corrected bug in call to "ioconf_startdate" in "gcm".

-- Added call to "ioconf_calendar" in "create_etat0_limit".

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/regr_lat_time_climoz_m.F90 (modified) (9 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/regr_lat_time_climoz_m.F90

@@ -44 +44 @@
 
     ! We assume that in the input file:
-    ! -- the latitude is in degrees and strictly monotonic (as all
-    ! NetCDF coordinate variables should be);
+    ! -- latitude is in degrees;
+    ! -- pressure is in hPa (even though we do not use pressure values
+    ! here, we write the unit of pressure in the NetCDF header, and we
+    ! will use the assumption later, when we regrid in pressure).
+    ! -- latitude and pressure are strictly monotonic (as all NetCDF
+    ! coordinate variables should be);
     ! -- time increases (even though we do not use values of the input
     ! time coordinate);
-    ! -- pressure is in hPa and in strictly ascending order (even
-    ! though we do not use pressure values here, we write the unit of
-    ! pressure in the NetCDF header, and we will use the assumptions later,
-    ! when we regrid in pressure).
 
     use regr1_step_av_m, only: regr1_step_av
@@ -81 +81 @@
     ! ascending order)
 
-    real, pointer:: plev(:) ! pressure level of input data
+    real, pointer:: plev(:)
+    ! pressure levels of input data, sorted in strictly ascending order
+
     logical desc_lat ! latitude in descending order in the input file
+    logical desc_plev ! pressure levels in descending order in the input file
 
     real, pointer:: o3_in(:, :, :) ! (n_lat, n_plev, 12 or 0:13)
@@ -110 +113 @@
 
     integer j
-
-    integer varid_in, varid_out, varid_plev, varid_time
-    integer varid
-    ! (for NetCDF)
+    integer varid_in, varid_out, varid_plev, varid_time, varid ! (for NetCDF)
 
     real, parameter:: tmidmonth(0:13) = (/(-15. + 30. * j, j = 0, 13)/)
@@ -153 +153 @@
     call nf95_gw_var(ncid_in, varid, plev)
     n_plev = size(plev)
-    ! (We only need the pressure coordinate to copy it to the output file.)
+    ! We only need the pressure coordinate to copy it to the output file.
+    ! The program "gcm" will assume that pressure levels are in
+    ! ascending order in the regridded climatology so invert order if
+    ! necessary:
+    desc_plev = plev(1) > plev(n_plev)
+    if (desc_plev) plev = plev(n_plev:1:-1)
 
     ! Create the output file and get the variable IDs:
@@ -168 +173 @@
     call nf95_gw_var(ncid_in, varid_in, o3_in)
     if (desc_lat) o3_in = o3_in(n_lat:1:-1, :, :)
+    if (desc_plev) o3_in = o3_in(:, n_plev:1:-1, :)
 
     call nf95_close(ncid_in)
@@ -219 +225 @@
     use netcdf95, only: nf95_create, nf95_def_dim, nf95_def_var, &
          nf95_put_att, nf95_enddef, nf95_copy_att, nf95_put_var
-    use netcdf, only: nf90_clobber, nf90_float, nf90_copy_att, nf90_global
+    use netcdf, only: nf90_clobber, nf90_float, nf90_global
 
     integer, intent(in):: ncid_in, n_plev
@@ -276 +282 @@
     ! Global attributes:
 
-    call nf95_put_att(ncid_out, nf90_global, "comment", "Regridded for LMDZ")
-
     ! The following commands, copying attributes, may fail.
     ! That is OK.
@@ -290 +294 @@
     call handle_err_copy_att("title")
 
+    call nf95_copy_att(ncid_in, nf90_global, "institution", ncid_out, &
+         nf90_global, ncerr)
+    call handle_err_copy_att("institution")
+
     call nf95_copy_att(ncid_in, nf90_global, "source", ncid_out, nf90_global, &
          ncerr)
     call handle_err_copy_att("source")
 
+    call nf95_put_att(ncid_out, nf90_global, "comment", "Regridded for LMDZ")
+
     call nf95_enddef(ncid_out)
 
@@ -311 +321 @@
 
       if (ncerr /= nf90_noerr) then
-         print *, "regr_lat_time_climoz_m prepare_out nf90_copy_att " &
+         print *, "regr_lat_time_climoz_m prepare_out nf95_copy_att " &
               // att_name // " -- " // trim(nf90_strerror(ncerr))
       end if