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

Timestamp:

May 4, 2009, 5:24:19 PM (16 years ago)

Author:

lguez

Message:

-- Added "NetCDF95" interface in "bibio".
-- NetCDF95 uses the module "typesizes", which is part of NetCDF, so we
exclude dependency on "typesizes" in "bld.cfg".
-- Added "assert_eq" and "assert" procedures, which are in the public
part of Numerical Recipes.
-- Added some interpolation and regridding utilities in "bibio".
-- Added the ability to read an ozone climatology from a NetCDF file.
-- Commented out unused variables and code in "etat0_netcdf".
-- Updated calls to NetCDF in "etat0_netcdf": from Fortran 77
interface to Fortran 90 interface.
-- Removed useless "deallocate" at the end of "etat0_netcdf".
-- Corrected some declarations not conforming to Fortran standard, such
as "integer*4", or obsolescent such as "character*4".
-- Replaced some calls to not-standard function "float" by calls to
"real".
-- On Brodie at IDRIS, the NetCDF library compiled with OpenMP should
be used. Changed path in "arch-SX8_BRODIE.path".
-- Added warning for incompatibility of debugging options and OpenMP
parallelization in "makelmdz_fcm".

File:

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

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

@@ -28 +28 @@
       USE fonte_neige_mod, ONLY  : fonte_neige_get_vars
       USE phys_output_mod
+      use open_climoz_m, only: open_climoz ! ozone climatology
+      use regr_pr, only: regr_pr_av
+      use netcdf95, only: nf95_close
+      use mod_phys_lmdz_mpi_data, only: is_mpi_root
 
       IMPLICIT none
@@ -710 +714 @@
       SAVE lmt_pas                ! frequence de mise a jour
 c$OMP THREADPRIVATE(lmt_pas) 
+      real zmasse(klon, llm) 
+C     (column-density of mass of air in a cell, in kg m-2)
+      real, parameter:: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
 
 cIM sorties
@@ -1093 +1100 @@
       LOGICAL,SAVE :: first=.true.
 c$OMP THREADPRIVATE(first)
+
+      logical, save::  read_climoz ! read ozone climatology
+      integer, save:: ncid_climoz ! NetCDF file containing ozone climatology
+
+      real, pointer, save:: press_climoz(:)
+!     edges of pressure intervals for ozone climatology, in Pa, in strictly
+!     ascending order
+
 #include "YOMCST.h"
 #include "YOETHF.h"
@@ -1139 +1154 @@
 
       torsfc=0.
+      forall (k=1: llm) zmasse(:, k) = (paprs(:, k)-paprs(:, k+1)) / rg
 
       if (first) then 
@@ -1234 +1250 @@
      .                  iflag_thermals_ed,iflag_thermals_optflux,
 cnv flags pour la convection et les poches froides
-     .                   iflag_coupl,iflag_clos,iflag_wake)
+     .                   iflag_coupl,iflag_clos,iflag_wake, read_climoz)
 
       print*,'iflag_coupl,iflag_clos,iflag_wake',
@@ -1555 +1571 @@
       call iniradia(klon,klev,paprs(1,1:klev+1))
 
+C$omp single
+      if (read_climoz) then
+         call open_climoz(ncid_climoz, press_climoz)
+      END IF
+C$omp end single
       ENDIF
 !
@@ -1723 +1744 @@
 c
       IF (MOD(itap-1,lmt_pas) .EQ. 0) THEN
-         if(prt_level.ge.1) WRITE(lunout,*)' PHYS cond  julien ',julien
-         CALL ozonecm( FLOAT(julien), rlat, paprs, wo)
+C        Once per day, update ozone:
+         if (read_climoz) then
+C           Ozone climatology from a NetCDF file
+            call regr_pr_av(ncid_climoz, "tro3", julien, press_climoz,
+     $           paprs, wo)
+!           Convert from mole fraction of ozone to column density of ozone in a
+!           cell, in kDU:
+            wo = wo * 48. / 29. * zmasse / dobson_u / 1e3
+C           (By regridding ozone values for LMDZ only once per day, we
+C           have already neglected the variation of pressure in one
+C           day. So do not recompute "wo" at each time step even if
+C           "zmasse" changes a little.)
+         else
+            CALL ozonecm(real(julien), rlat, paprs, wo)
+         end if
       ENDIF
 c
@@ -3496 +3530 @@
 !         write(97) u_seri,v_seri,t_seri,q_seri
 !         close(97)
+C$OMP single
+         if (read_climoz) then
+            if (is_mpi_root) then
+               call nf95_close(ncid_climoz)
+            end if
+            deallocate(press_climoz) ! pointer
+         end if
+C$OMP end single nowait
       ENDIF