Changeset 1263 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F

Timestamp:

Nov 17, 2009, 2:00:14 PM (15 years ago)

Author:

lguez

Message:

1) Reactivated ability to read ozone (that was deactivated because of
dependency on version of IOIPSL). Added ability to read a pressure
coordinate in Pa in "regr_lat_time_climoz".

2) Added the ability to read a second ozone climatology, corresponding to
daylight ozone:

-- "read_climoz" is now an integer variable, instead of a logical
variable.

-- Added argument "read_climoz" to "phys_state_var_init",
"phys_output_open" and "regr_lat_time_climoz".

-- Created new variable "ozone_daylight" for "hist*.nc" output files.

-- Added a third dimension to variable "wo" in module
"phys_state_var_mod" and variable "POZON" in "radlwsw": index 1 for
average day-night ozone, index 2 for daylight ozone.

-- Added a fourth dimension to variables "o3_in", "o3_regr_lat" and
"o3_out" in "regr_lat_time_climoz": index 1 for average day-night
ozone, index 2 for daylight ozone.

-- In "physiq", moved call to "conf_phys" before call to
"phys_state_var_init". Thus, "conf_phys" is now inside the block "if
(first)" instead of "IF (debut)". There were definitions of "bl95_b0"
and "bl95_b1" that were useless because the variables were overwritten
by "conf_phys". Removed those definitions.

-- In "radlwsw", we pass the average day-night ozone to "LW_LMDAR4"
and the daylight ozone, if we have it, to "SW_LMDAR4" or
"SW_AEROAR4". If we do not have a specific field for daylight ozone
then "SW_LMDAR4" or "SW_AEROAR4" just get the average day-night ozone.

-- "regr_lat_time_climoz" now manages latitudes where the input ozone
field is missing at all levels (polar night).

-- Encapsulated "radlwsw" in a module.

3) Modifications to make sequential and parallel versions of
"create_etat0_limit" almost identical:

-- In "dyn3dpar/create_etat0_limit.F". No need to call
"phys_state_var_init", removed "use phys_state_var_mod" statement. No
need for "clesphys.h", removed "include" statement.

-- In "dyn3dpar/etat0_netcdf.F". Added argument "tau_ratqs" in call to
"conf_phys" (this bug was already corrected in "dyn3d"). Moved call to
"inifilr" after call to "infotrac_init" (as in "dyn3d").

4) Other peripheral modifications:

-- Added procedures "nf95_get_att" and "nf95_def_var_scalar" in
NetCDF95 interface. Overloaded "nf95_put_var" with three more
procedures: "nf95_put_var_FourByteReal", "nf95_put_var_FourByteInt",
"nf95_put_var_1D_FourByteInt".

-- Overloaded "regr1_step_av" with one more procedure:
"regr14_step_av". Overloaded "regr3_lint" with one more procedure:
"regr34_lint".

-- Corrected call to "Init_Phys_lmdz" in "dyn3d/create_etat0_limit.F":
the last argument should be an array, not a scalar.

-- Encapsulated "conf_phys" in a module.

-- Splitted module "regr_pr" into "regr_pr_av_m" and "regr_pr_int_m".

5) Tests:

This revision was compared to revision 1259, with optimization options
"debug" and "dev", parallelization options "none", "mpi", "omp" and
"mpi_omp", 1 and 2 MPI processes, 1 and 2 OpenMP threads, with the
compiler "FORTRAN90/SX Version 2.0 for SX-8". Both programs
"create_etat0_limit" and "gcm" were tested. In all cases,
parallelization does not change the results. With "read_climoz = 0" in
the ".def" files, the results of revision 1259 and of this revision
are the same.

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F (modified) (16 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F

@@ -13 +13 @@
 
       USE ioipsl, only: histbeg, histvert, histdef, histend, histsync,
-     $     histwrite, ju2ymds, ymds2ju
-!JG     $     histwrite, ju2ymds, ymds2ju, ioget_year_len
+     $     histwrite, ju2ymds, ymds2ju, ioget_year_len
       USE comgeomphy
       USE phys_cal_mod
@@ -33 +32 @@
       USE phys_output_mod
       use open_climoz_m, only: open_climoz ! ozone climatology from a file
-      use regr_pr, only: regr_pr_av
+      use regr_pr_av_m, only: regr_pr_av
       use netcdf95, only: nf95_close
       use mod_phys_lmdz_mpi_data, only: is_mpi_root
       USE aero_mod
       use ozonecm_m, only: ozonecm ! ozone of J.-F. Royer
+      use conf_phys_m, only: conf_phys
+      use radlwsw_m, only: radlwsw
 
       IMPLICIT none
@@ -1093 +1094 @@
       integer iunit
 
-      logical, save::  read_climoz ! read ozone climatology from a file
-      integer, save:: ncid_climoz ! NetCDF file containing ozone climatology
+      integer, save::  read_climoz ! read ozone climatology
+C     Allowed values are 0, 1 and 2
+C     0: do not read an ozone climatology
+C     1: read a single ozone climatology that will be used day and night
+C     2: read two ozone climatologies, the average day and night
+C     climatology and the daylight climatology
+
+      integer, save:: ncid_climoz ! NetCDF file containing ozone climatologies
 
       real, pointer, save:: press_climoz(:)
-!     edges of pressure intervals for ozone climatology, in Pa, in strictly
+!     edges of pressure intervals for ozone climatologies, in Pa, in strictly
 !     ascending order
 
-      integer,save:: co3i = 0 ! time index in NetCDF file of current ozone field
+      integer, save:: co3i = 0
+!     time index in NetCDF file of current ozone fields
 c$OMP THREADPRIVATE(co3i) 
 
       integer ro3i
-!     required time index in NetCDF file for the ozone field, between 1 and 360
+!     required time index in NetCDF file for the ozone fields, between 1
+!     and 360
 
 #include "YOMCST.h"
@@ -1165 +1174 @@
       print*, 'Allocation des variables locales et sauvegardees'
       call phys_local_var_init
-      call phys_state_var_init
+c     appel a la lecture du run.def physique
+      call conf_phys(ok_journe, ok_mensuel,
+     .     ok_instan, ok_hf,
+     .     ok_LES,
+     .     solarlong0,seuil_inversion,
+     .     fact_cldcon, facttemps,ok_newmicro,iflag_radia,
+     .     iflag_cldcon,iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs,
+     .     ok_ade, ok_aie, aerosol_couple, 
+     .     flag_aerosol, new_aod,
+     .     bl95_b0, bl95_b1,
+     .     iflag_thermals,nsplit_thermals,tau_thermals,
+     .     iflag_thermals_ed,iflag_thermals_optflux,
+c     nv flags pour la convection et les poches froides
+     .     iflag_coupl,iflag_clos,iflag_wake, read_climoz)
+      call phys_state_var_init(read_climoz)
       print*, '================================================='
 
@@ -1181 +1204 @@
         first=.false.
 
-      endif  ! fisrt
+      endif  ! first
 
        modname = 'physiq'
@@ -1209 +1232 @@
 c
        IF (debut) THEN
-C
 !rv
 cCRinitialisation de wght_th et lalim_conv pour la definition de la couche alimentation 
@@ -1220 +1242 @@
          rain_con(:)=0.
          snow_con(:)=0.
-         bl95_b0=0.
-         bl95_b1=0.
          topswai(:)=0.
          topswad(:)=0.
@@ -1247 +1267 @@
          IF (ip_ebil_phy.ge.1) d_h_vcol_phy=0.
 c
-c appel a la lecture du run.def physique
-c
-         call conf_phys(ok_journe, ok_mensuel,
-     .                  ok_instan, ok_hf,
-     .                  ok_LES,
-     .                  solarlong0,seuil_inversion,
-     .                  fact_cldcon, facttemps,ok_newmicro,iflag_radia,
-     .         iflag_cldcon,iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs,
-     .                  ok_ade, ok_aie, aerosol_couple, 
-     .                  flag_aerosol, new_aod,
-     .                  bl95_b0, bl95_b1,
-     .                  iflag_thermals,nsplit_thermals,tau_thermals,
-     .                  iflag_thermals_ed,iflag_thermals_optflux,
-cnv flags pour la convection et les poches froides
-     .                   iflag_coupl,iflag_clos,iflag_wake, read_climoz)
-
       print*,'iflag_coupl,iflag_clos,iflag_wake',
      .   iflag_coupl,iflag_clos,iflag_wake
@@ -1453 +1457 @@
      &                        ctetaSTD,dtime,ok_veget,
      &                        type_ocean,iflag_pbl,ok_mensuel,ok_journe,
-     &                        ok_hf,ok_instan,ok_LES,ok_ade,ok_aie)
+     &              ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, read_climoz)
 c$OMP END MASTER
 c$OMP BARRIER
@@ -1544 +1548 @@
 
 C$omp single
-      if (read_climoz) then
+      if (read_climoz >= 1) then
          call open_climoz(ncid_climoz, press_climoz)
       END IF
@@ -1711 +1715 @@
 c Prescrire l'ozone et calculer l'albedo sur l'ocean.
 c
-      if (read_climoz) then
+      if (read_climoz >= 1) then
 C        Ozone from a file
 !        Update required ozone index:
-!JG : ioget_year_len n'existe pas dans les versions officiels d'ioipsl
-!JG         ro3i = int((days_elapsed + jh_cur - jh_1jan)
-!JG     $        / ioget_year_len(year_cur) * 360.) + 1
-         CALL abort_gcm(modname,
-     $     'JG : read_climoz temporary desactivated',1)
-
+         ro3i = int((days_elapsed + jh_cur - jh_1jan)
+     $        / ioget_year_len(year_cur) * 360.) + 1
          if (ro3i == 361) ro3i = 360
 C        (This should never occur, except perhaps because of roundup
@@ -1725 +1725 @@
          if (ro3i /= co3i) then
 C           Update ozone field:
-            call regr_pr_av(ncid_climoz, "tro3", julien=ro3i, 
-     &           press_in_edg=press_climoz, paprs=paprs, v3=wo)
+            if (read_climoz == 1) then
+               call regr_pr_av(ncid_climoz, (/"tro3"/), julien=ro3i,
+     $              press_in_edg=press_climoz, paprs=paprs, v3=wo)
+            else
+C              read_climoz == 2
+               call regr_pr_av(ncid_climoz,
+     $              (/"tro3         ", "tro3_daylight"/),
+     $              julien=ro3i, press_in_edg=press_climoz, paprs=paprs,
+     $              v3=wo)
+            end if
 !           Convert from mole fraction of ozone to column density of ozone in a
 !           cell, in kDU:
-            wo = wo * rmo3 / rmd * zmasse / dobson_u / 1e3
+            forall (l = 1: read_climoz) wo(:, :, l) = wo(:, :, l)
+     $           * rmo3 / rmd * zmasse / dobson_u / 1e3
 C           (By regridding ozone values for LMDZ only once every 360th of
 C           year, we have already neglected the variation of pressure in one
@@ -1738 +1747 @@
       elseif (MOD(itap-1,lmt_pas) == 0) THEN
 C        Once per day, update ozone from Royer:
-         wo = ozonecm(rlat, paprs, rjour=real(days_elapsed+1))
+         wo(:, :, 1) = ozonecm(rlat, paprs, rjour=real(days_elapsed+1))
       ENDIF
 c
@@ -2836 +2845 @@
      $                          u,
      $                          v,
-     $                          wo,
+     $                          wo(:, :, 1),
      $                          q_seri,
      $                          zxtsol,
@@ -2919 +2928 @@
      e        (kdlon,kflev,dist, rmu0, fract, solaire,
      e        paprs, pplay,zxtsol,albsol1, albsol2, t_seri,q_seri,
-     e        wo,
+     e        wo(:, :, 1),
      e        cldfra, cldemi, cldtau,
      s        heat,heat0,cool,cool0,radsol,albpla,
@@ -3529 +3538 @@
 !         close(97)
 C$OMP MASTER
-         if (read_climoz) then
+         if (read_climoz >= 1) then
             if (is_mpi_root) then
                call nf95_close(ncid_climoz)