Changeset 1154 for LMDZ4/branches/LMDZ4-dev/libf/phylmd

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".

Location:

LMDZ4/branches/LMDZ4-dev/libf/phylmd

Files:

• clesphys.h (modified) (2 diffs)
• conf_phys.F90 (modified) (4 diffs)
• indicesol.h (modified) (1 diff)
• open_climoz_m.f90 (added)
• physiq.F (modified) (8 diffs)
• regr_lat_time_climoz_m.f90 (added)
• regr_pr.f90 (added)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/clesphys.h

@@ -13 +13 @@
        INTEGER nbapp_rad, iflag_con
        REAL co2_ppm, solaire
-       REAL*8 RCO2, RCH4, RN2O, RCFC11, RCFC12  
-       REAL*8 CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt
+       REAL(kind=8) RCO2, RCH4, RN2O, RCFC11, RCFC12  
+       REAL(kind=8) CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt
 
 !OM ---> correction du bilan d'eau global
@@ -39 +39 @@
 !IM lev_histmth : niveau sorties mensuelles
        INTEGER lev_histhf, lev_histday, lev_histmth
-       CHARACTER*4 type_run
+       CHARACTER(len=4) type_run
 ! aer_type: pour utiliser un fichier constant dans readaerosol 
        CHARACTER*8 :: aer_type 

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

@@ -7 +7 @@
 
   subroutine 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, &
- &                     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, &
- &                     iflag_coupl,iflag_clos,iflag_wake )
+                       ok_LES,&
+                       solarlong0,seuil_inversion, &
+                       fact_cldcon, facttemps,ok_newmicro,iflag_radia,&
+                       iflag_cldcon, &
+                       iflag_ratqs,ratqsbas,ratqshaut, &
+                       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, &
+                       iflag_coupl,iflag_clos,iflag_wake, read_climoz)
 
    use IOIPSL
@@ -140 +140 @@
   LOGICAL,SAVE :: ok_strato_omp
   LOGICAL,SAVE :: ok_hines_omp
+
+  logical, intent(out):: read_climoz ! read ozone climatology
 !
 
@@ -1221 +1223 @@
   call getin('ecrit_LES', ecrit_LES_omp)
 !
+  read_climoz = .false. ! default value
+  call getin('read_climoz', read_climoz)
 
 !$OMP END MASTER
@@ -1467 +1471 @@
 
   write(numout,*)' lonmin lonmax latmin latmax bilKP_ins =',&
- & lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
+   lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
   write(numout,*)' ecrit_ hf, day, mth, reg, tra, ISCCP, LES',&
- & ecrit_hf, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra, ecrit_ISCCP, ecrit_LES
+   ecrit_hf, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra, ecrit_ISCCP, ecrit_LES
 
   write(numout,*) 'ok_strato = ', ok_strato
   write(numout,*) 'ok_hines = ',  ok_hines
+  write(numout,*) 'read_climoz = ', read_climoz
   
 !$OMP END MASTER

LMDZ4/branches/LMDZ4-dev/libf/phylmd/indicesol.h

@@ -24 +24 @@
       PARAMETER (epsfra=1.0E-05)
 !
-      CHARACTER *3 clnsurf(nbsrf)
+      CHARACTER(len=3) clnsurf(nbsrf)
       DATA clnsurf/'ter', 'lic', 'oce', 'sic'/
       SAVE clnsurf

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