Changeset 1279 for LMDZ4/trunk/libf/dyn3d/etat0_netcdf.F

Timestamp:

Dec 10, 2009, 10:02:56 AM (15 years ago)

Author:

Laurent Fairhead

Message:

Merged LMDZ4-dev branch changes r1241:1278 into the trunk
Running trunk and LMDZ4-dev in LMDZOR configuration on local
machine (sequential) and SX8 (4-proc) yields identical results
(restart and restartphy are identical binarily)
Log history from r1241 to r1278 is available by switching to
source:LMDZ4/branches/LMDZ4-dev-20091210

Location:

LMDZ4/trunk

Files:

• . (modified) (1 prop)
• libf/dyn3d/etat0_netcdf.F (modified) (16 diffs)

LMDZ4/trunk/libf/dyn3d/etat0_netcdf.F

@@ -14 +14 @@
       USE phys_state_var_mod
       USE filtreg_mod
+      use regr_lat_time_climoz_m, only: regr_lat_time_climoz
+      use conf_phys_m, only: conf_phys
 #endif
 !#endif of #ifdef CPP_EARTH
+      use netcdf, only: nf90_open, NF90_NOWRITE, nf90_close
       !
       IMPLICIT NONE
       !
-#include "netcdf.inc"
 #include "dimensions.h"
 #include "paramet.h"
@@ -49 +51 @@
       REAL :: vvent(iip1, jjm, llm)
       REAL :: t3d(iip1, jjp1, llm), tpot(iip1, jjp1, llm)
-      REAL, ALLOCATABLE, DIMENSION(:,:,:,:) :: q3d
       REAL :: qsat(iip1, jjp1, llm)
+      REAL,ALLOCATABLE :: q3d(:, :, :,:)
       REAL :: tsol(klon), qsol(klon), sn(klon)
-      REAL :: tsolsrf(klon,nbsrf), qsolsrf(klon,nbsrf),snsrf(klon,nbsrf) 
+!!      REAL :: tsolsrf(klon,nbsrf)
+      real qsolsrf(klon,nbsrf),snsrf(klon,nbsrf) 
       REAL :: albe(klon,nbsrf), evap(klon,nbsrf)
       REAL :: alblw(klon,nbsrf)
@@ -72 +75 @@
       !
 
-      CHARACTER*80 :: varname
+      CHARACTER(len=80) :: varname
       !
       INTEGER :: i,j, ig, l, ji,ii1,ii2
@@ -102 +105 @@
       REAL :: w(ip1jmp1,llm)
       REAL ::phystep
-      REAL :: rugsrel(iip1*jjp1)
+CC      REAL :: rugsrel(iip1*jjp1)
       REAL :: fder(klon)
-      real zrel(iip1*jjp1),chmin,chmax
-
-      CHARACTER*80 :: visu_file
+!!      real zrel(iip1*jjp1),chmin,chmax
+
+!!      CHARACTER(len=80) :: visu_file
       INTEGER :: visuid
 
@@ -126 +129 @@
       logical              :: ok_journe, ok_mensuel, ok_instan, ok_hf
       logical              :: ok_LES
-      LOGICAL              :: ok_ade, ok_aie, aerosol_couple
+      LOGICAL              :: ok_ade, ok_aie, aerosol_couple, new_aod
+      INTEGER              :: flag_aerosol
       REAL                 :: bl95_b0, bl95_b1
       real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut
+      real                 :: tau_ratqs
       integer              :: iflag_cldcon
       integer              :: iflag_ratqs
@@ -140 +145 @@
       real :: seuil_inversion
 
+      integer  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
+
       !
       !   Constantes 
@@ -162 +174 @@
 !      CALL defrun_new(99,.TRUE.,clesphy0)
       CALL conf_gcm( 99, .TRUE. , clesphy0 )
-      call conf_phys(ok_journe, ok_mensuel, ok_instan,                  &
-     &                 ok_hf, ok_LES,                                   &
+      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,                  &
+     &                 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,        &
-     &                 iflag_coupl,iflag_clos,iflag_wake )
+     &                 iflag_coupl,iflag_clos,iflag_wake, read_climoz )
+
+! co2_ppm0 : initial value of atmospheric CO2 from .def file (co2_ppm value)
+      co2_ppm0 = co2_ppm
+
       dtvr   = daysec/FLOAT(day_step)
       print*,'dtvr',dtvr
 
-
-
       CALL iniconst()
       CALL inigeom()
 
 ! Initialisation pour traceurs
-      CALL infotrac_init
-      ALLOCATE(q3d(iip1,jjp1,llm,nqtot))
-
+      call infotrac_init
+      ALLOCATE(q3d(iip1, jjp1, llm, nqtot))
 
       CALL inifilr()
-      CALL phys_state_var_init()
+      CALL phys_state_var_init(read_climoz)
       !
       latfi(1) = ASIN(1.0)
@@ -244 +257 @@
 
       write(*,*)'Essai de lecture masque ocean'
-      iret = nf_open("o2a.nc", NF_NOWRITE, nid_o2a)
+      iret = nf90_open("o2a.nc", NF90_NOWRITE, nid_o2a)
       if (iret .ne. 0) then
         write(*,*)'ATTENTION!! pas de fichier o2a.nc trouve'
@@ -263 +276 @@
       else
         couple = .true.
-        iret = nf_close(nid_o2a)
+        iret = nf90_close(nid_o2a)
         call flininfo("o2a.nc", iml_omask, jml_omask, llm_tmp, ttm_tmp
      $    , nid_o2a)
@@ -400 +413 @@
      .                           maxval(qsat(:,:,:))
       !
-      WRITE(*,*) 'QSAT :', qsat(10,20,:)
+CC      WRITE(*,*) 'QSAT :', qsat(10,20,:)
       !
       varname = 'q'
@@ -411 +424 @@
       q3d(:,:,:,1) = qd(:,:,:)
       !
+
+!     Ozone climatology:
+      if (read_climoz >= 1) call regr_lat_time_climoz(read_climoz)
+
       varname = 'tsol'
       ! This line needs to be replaced by a call to restget to get the values in the restart file
@@ -475 +492 @@
      .     jjm, rlonu, rlatv , interbar )
 c
-      rugsrel(:) = 0.0
-      IF(ok_orodr)  THEN
-        DO i = 1, iip1* jjp1
-         rugsrel(i) = MAX( 1.e-05, zstd(i)* zsig(i) /2. )
-        ENDDO
-      ENDIF
+cc      rugsrel(:) = 0.0
+cc      IF(ok_orodr)  THEN
+cc        DO i = 1, iip1* jjp1
+cc         rugsrel(i) = MAX( 1.e-05, zstd(i)* zsig(i) /2. )
+cc        ENDDO
+cc      ENDIF
 
 
@@ -650 +667 @@
       itau_phy = 0
       iday = dayref +itau/day_step
-      time = FLOAT(itau-(iday-dayref)*day_step)/day_step
+      time = real(itau-(iday-dayref)*day_step)/day_step
 c     
       IF(time.GT.1)  THEN
@@ -714 +731 @@
       q_ancien = 0.
       agesno = 0.
+c
       frugs(1:klon,is_oce) = rugmer(1:klon)
       frugs(1:klon,is_ter) = MAX(1.0e-05, zstd(1:klon)*zsig(1:klon)/2.0)
@@ -750 +768 @@
 
 C     Sortie Visu pour les champs dynamiques
-      if (1.eq.0 ) then
-      print*,'sortie visu'
-      time_step = 1.
-      t_ops = 2.
-      t_wrt = 2.
-      itau = 2.
-      visu_file='Etat0_visu.nc'
-      CALL initdynav(visu_file,dayref,anneeref,time_step,
-     .              t_ops, t_wrt, visuid)
-      CALL writedynav(visuid, itau,vvent ,
-     .                uvent,tpot,pk,phi,q3d,masse,psol,phis)
-      else
+cc      if (1.eq.0 ) then
+cc      print*,'sortie visu'
+cc      time_step = 1.
+cc      t_ops = 2.
+cc      t_wrt = 2.
+cc      itau = 2.
+cc      visu_file='Etat0_visu.nc'
+cc      CALL initdynav(visu_file,dayref,anneeref,time_step,
+cc     .              t_ops, t_wrt, visuid)
+cc      CALL writedynav(visuid, itau,vvent ,
+cc     .                uvent,tpot,pk,phi,q3d,masse,psol,phis)
+cc      else
          print*,'CCCCCCCCCCCCCCCCCC REACTIVER SORTIE VISU DANS ETAT0'
-      endif
+cc      endif
       print*,'entree histclo'
       CALL histclo
-
-      DEALLOCATE(q3d)
 
 #endif 
@@ -774 +790 @@
       !
       END SUBROUTINE etat0_netcdf
-