Changeset 927 for LMDZ4

Timestamp:

Mar 14, 2008, 1:12:38 PM (17 years ago)

Author:

lmdzadmin

Message:

Ajout variables zmax0, f0 dans le startphy.nc FH
IM

Location:

LMDZ4/trunk/libf

Files:

• dyn3d/etat0_netcdf.F (modified) (3 diffs)
• phylmd/calltherm.F90 (modified) (3 diffs)
• phylmd/phyetat0.F (modified) (5 diffs)
• phylmd/phyredem.F (modified) (4 diffs)
• phylmd/phys_state_var_mod.F90 (modified) (3 diffs)
• phylmd/physiq.F (modified) (3 diffs)
• phylmd/thermcell_main.F90 (modified) (5 diffs)

LMDZ4/trunk/libf/dyn3d/etat0_netcdf.F

@@ -120 +120 @@
 
       REAL, DIMENSION(klon,klev,nbsrf) :: pbl_tke ! turb kinetic energy
+      REAL, DIMENSION(klon) :: zmax0, f0 
       REAL :: dummy(klon, nbsrf)
       !
@@ -687 +688 @@
       ocean = "slab"
       pbl_tke(:,:,:) = 1.e-8 
+      zmax0(:) = 40.
+      f0(:) = 1.e-5
 
       call fonte_neige_init(run_off_lic_0)
@@ -698 +701 @@
      $    radsol, 
      $    zmea, zstd, zsig, zgam, zthe, zpic, zval, rugsrel, 
-     $    t_ancien, q_ancien, rnebcon, ratqs, clwcon, pbl_tke )
+     $    t_ancien, q_ancien, rnebcon, ratqs, clwcon, pbl_tke,
+     $    zmax0, f0 )
 
 

LMDZ4/trunk/libf/phylmd/calltherm.F90

@@ -7 +7 @@
      &      ,d_u_ajs,d_v_ajs,d_t_ajs,d_q_ajs  &
      &      ,fm_therm,entr_therm,zqasc,clwcon0,lmax,ratqscth,  &
-     &       ratqsdiff,zqsatth,Ale_bl,Alp_bl,lalim_conv,wght_th)
+     &       ratqsdiff,zqsatth,Ale_bl,Alp_bl,lalim_conv,wght_th, &
+     &       zmax0,f0)
 
       implicit none
@@ -59 +60 @@
       real Alp(klon)
 !RC
+      !on garde le zmax du pas de temps precedent
+      real zmax0(klon), f0(klon)
 !********************************************************
 
@@ -171 +174 @@
      &      ,ratqscth,ratqsdiff,zqsatth  &
      &      ,r_aspect_thermals,l_mix_thermals,w2di_thermals  &
-     &      ,tho_thermals,Ale,Alp,lalim_conv,wght_th)
+     &      ,tho_thermals,Ale,Alp,lalim_conv,wght_th &
+     &      ,zmax0,f0)
          endif
 

LMDZ4/trunk/libf/phylmd/phyetat0.F

@@ -13 +13 @@
      .           rugsrel_p,tabcntr0,
      .           t_ancien_p,q_ancien_p,ancien_ok_p, rnebcon_p, ratqs_p,
-     .           clwcon_p,pbl_tke_p)
+     .           clwcon_p,pbl_tke_p, zmax0_p, f0_p)
 
       USE dimphy
@@ -48 +48 @@
       REAL tsol_p(klon,nbsrf)
       REAL pbl_tke_p(klon,klev,nbsrf)
+      REAL zmax0_p(klon), f0_p(klon)
       REAL tsoil_p(klon,nsoilmx,nbsrf)
       REAL tslab_p(klon), seaice_p(klon)
@@ -85 +86 @@
       REAL tsol(klon_glo,nbsrf)
       REAL pbl_tke(klon_glo,klev,nbsrf)
+      REAL zmax0(klon), f0(klon)
       REAL tsoil(klon_glo,nsoilmx,nbsrf)
 cIM "slab" ocean
@@ -1538 +1540 @@
          ENDDO
       ENDIF
-
+c
+c zmax0
+      ierr = NF_INQ_VARID (nid, "ZMAX0", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, "phyetat0: Le champ <ZMAX0> est absent"
+         PRINT*, "Depart legerement fausse. Mais je continue"
+         zmax0=40.
+      ELSE
+#ifdef NC_DOUBLE
+         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zmax0)
+#else
+         ierr = NF_GET_VAR_REAL(nid, nvarid, zmax0)
+#endif
+         IF (ierr.NE.NF_NOERR) THEN
+            PRINT*, "phyetat0: Lecture echouee pour <ZMAX0>"
+            CALL abort
+         ENDIF
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      xmin = MINval(zmax0)
+      xmax = MAXval(zmax0)
+      PRINT*,'(ecart-type) zmax0:', xmin, xmax
+c
+c           f0(ig)=1.e-5
+c f0
+      ierr = NF_INQ_VARID (nid, "f0", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, "phyetat0: Le champ <f0> est absent"
+         PRINT*, "Depart legerement fausse. Mais je continue"
+         f0=1.e-5
+      ELSE
+#ifdef NC_DOUBLE
+         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, f0)
+#else
+         ierr = NF_GET_VAR_REAL(nid, nvarid, f0)
+#endif
+         IF (ierr.NE.NF_NOERR) THEN
+            PRINT*, "phyetat0: Lecture echouee pour <f0>"
+            CALL abort
+         ENDIF
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      xmin = MINval(f0)
+      xmax = MAXval(f0)
+      PRINT*,'(ecart-type) f0:', xmin, xmax
 c
 c Fermer le fichier:
@@ -1575 +1623 @@
          call Scatter( pbl_tke,pbl_tke_p)
       endif
+      call Scatter( zmax0,zmax0_p)
+      call Scatter( f0,f0_p)
       call Scatter( tsoil,tsoil_p)
       call Scatter( tslab,tslab_p)

LMDZ4/trunk/libf/phylmd/phyredem.F

@@ -10 +10 @@
      .           zgam_p,zthe_p,zpic_p,zval_p,rugsrel_p,
      .           t_ancien_p, q_ancien_p, rnebcon_p, ratqs_p, clwcon_p,
-     .           pbl_tke_p)
+     .           pbl_tke_p, zmax0_p, f0_p)
 
       USE dimphy
@@ -39 +39 @@
       REAL tsol_p(klon,nbsrf)
       REAL pbl_tke_p(klon,klev,nbsrf)
+      REAL zmax0_p(klon), f0_p(klon)
       REAL tsoil_p(klon,nsoilmx,nbsrf)
       CHARACTER*6 ocean
@@ -73 +74 @@
       REAL tsol(klon_glo,nbsrf)
       REAL pbl_tke(klon_glo,klev,nbsrf)
+      REAL zmax0(klon_glo), f0(klon)
       REAL tsoil(klon_glo,nsoilmx,nbsrf)
       REAL tslab(klon_glo), seaice(klon_glo)
@@ -139 +141 @@
       call Gather( tsol_p,tsol)
       call Gather( pbl_tke_p,pbl_tke)
+      call Gather( zmax0_p,zmax0)
+      call Gather( f0_p,f0)
       call Gather( tsoil_p,tsoil)
       call Gather( tslab_p,tslab)

LMDZ4/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -29 +29 @@
       REAL, ALLOCATABLE, SAVE :: ratqs(:,:)
       REAL, ALLOCATABLE, SAVE :: pbl_tke(:,:,:) ! turb kinetic energy
+      REAL, ALLOCATABLE, SAVE :: zmax0(:), f0(:) ! 
 
 CONTAINS
@@ -56 +57 @@
       ALLOCATE( ratqs(klon,klev))
       ALLOCATE(pbl_tke(klon,klev+1,nbsrf))
+      ALLOCATE(zmax0(klon), f0(klon))
 
 END SUBROUTINE phys_state_var_init
@@ -72 +74 @@
       deallocate( rugoro, t_ancien, q_ancien, clwcon, rnebcon )
       deallocate( ratqs, pbl_tke )
+      deallocate( zmax0, f0 )
 
 END SUBROUTINE phys_state_var_end

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -1640 +1640 @@
      .       zmea,zstd,zsig,zgam,zthe,zpic,zval,rugoro,tabcntr0,
      .       t_ancien, q_ancien, ancien_ok, rnebcon, ratqs,clwcon,
-     .       pbl_tke)
+     .       pbl_tke, zmax0, f0)
 
 
@@ -2645 +2645 @@
      s      ,ratqsdiff,zqsatth
 con rajoute ale et alp, et les caracteristiques de la couche alim
-     s      ,Ale_bl,Alp_bl,lalim_conv,wght_th)
+     s      ,Ale_bl,Alp_bl,lalim_conv,wght_th, zmax0, f0)
          endif
 
@@ -3757 +3757 @@
      .      zmea,zstd,zsig,zgam,zthe,zpic,zval,rugoro,
      .      t_ancien, q_ancien, rnebcon, ratqs, clwcon,
-     .      pbl_tke)
+     .      pbl_tke, zmax0, f0)
          open(97,form="unformatted",file="finbin")
          write(97) u_seri,v_seri,t_seri,q_seri

LMDZ4/trunk/libf/phylmd/thermcell_main.F90

@@ -9 +9 @@
      &                  ,ratqscth,ratqsdiff,zqsatth  &
      &                  ,r_aspect,l_mix,w2di,tho &
-     &                  ,Ale_bl,Alp_bl,lalim_conv,wght_th)
+     &                  ,Ale_bl,Alp_bl,lalim_conv,wght_th &
+     &                  ,zmax0, f0)
 
       IMPLICIT NONE
@@ -74 +75 @@
 !on garde le zmax du pas de temps precedent
       real zmax0(klon)
-      save zmax0
+!FH/IM     save zmax0
 
       real zlev(klon,klev+1),zlay(klon,klev)
@@ -130 +131 @@
       real alim_star(klon,klev)
       real f(klon), f0(klon)
-      save f0
+!FH/IM     save f0
       real zlevinter(klon)
       logical debut
@@ -174 +175 @@
 !
       do ig=1,klon      
-      if ((debut).or.((.not.debut).and.(f0(ig).lt.1.e-10))) then
+!FH/IM 130308     if ((debut).or.((.not.debut).and.(f0(ig).lt.1.e-10))) then
+      if ((.not.debut).and.(f0(ig).lt.1.e-10)) then
             f0(ig)=1.e-5
             zmax0(ig)=40.
@@ -180 +182 @@
       endif
       enddo 
-