Changeset 3825 for dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/regr_lat_time_climoz_m.F90

Timestamp:

May 6, 2015, 12:14:12 PM (10 years ago)

Author:

ymipsl

Message:

Reorganize geometry and grid modules. Prepare physics for unstructutured grid support. Simplify initialization of physics from dynamic.
Compiled only with dynd3dmem, but not tested for moment.

YM

File:

• dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/regr_lat_time_climoz_m.F90 (modified) (5 diffs)

dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/regr_lat_time_climoz_m.F90

@@ -73 +73 @@
     use netcdf, only: nf90_get_att, nf90_get_var, nf90_noerr, nf90_nowrite
     use assert_m, only: assert
-    !use comconst_phy_mod, only : pi
-    !use comgeom2_phy_mod, only : rlatv
-    use comgeomphy, only: rlatv
+    use regular_lonlat_mod, only : boundslat_reg, south
     use nrtype, only: pi
     integer, intent(in):: read_climoz ! read ozone climatology
@@ -293 +291 @@
             "Found 12 months in ozone climatologies, assuming periodicity..."
        o3_regr_lat(nbp_lat:1:-1, :, 1:12, :) = regr1_step_av(o3_in, &
-            xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(nbp_lat-1:1:-1), pi / 2/)))
+            xs=sin(lat_in_edg), xt=sin((/- pi / 2, boundslat_reg(nbp_lat-1:1:-1,south), pi / 2/)))
        ! (invert order of indices in "o3_regr_lat" because "rlatu" is
        ! in descending order)
@@ -304 +302 @@
        print *, "Using 14 months in ozone climatologies..."
        o3_regr_lat(nbp_lat:1:-1, :, :, :) = regr1_step_av(o3_in, &
-            xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(nbp_lat-1:1:-1), pi / 2/)))
+            xs=sin(lat_in_edg), xt=sin((/- pi / 2, boundslat_reg(nbp_lat-1:1:-1,south), pi / 2/)))
        ! (invert order of indices in "o3_regr_lat" because "rlatu" is
        ! in descending order)
@@ -334 +332 @@
          nf95_put_att, nf95_enddef, nf95_copy_att, nf95_put_var
     use netcdf, only: nf90_clobber, nf90_float, nf90_global
-    !use comconst_phy_mod, only : pi
-    !use comgeom2_phy_mod, only : rlatu
     use nrtype, only: pi
-    use comgeomphy, only: rlatu
+    use regular_lonlat_mod, only : lat_reg
     integer, intent(in):: ncid_in, n_plev
     integer, intent(out):: ncid_out, varid_plev, varid_time
@@ -421 +417 @@
 
     ! Write one of the coordinate variables:
-    call nf95_put_var(ncid_out, varid_rlatu, rlatu(nbp_lat:1:-1) / pi * 180.)
+    call nf95_put_var(ncid_out, varid_rlatu, lat_reg(nbp_lat:1:-1) / pi * 180.)
     ! (convert from rad to degrees and sort in ascending order)