Changeset 996 for LMDZ4/trunk/libf/phylmd/physiq.F

Timestamp:

Sep 9, 2008, 3:22:23 PM (16 years ago)

Author:

lsce

Message:

• Modifications liées au calcul des nouveau sous-fractions
• Nettoyage de ocean slab : il reste uniquement la version avec glace de mer forcé
• Nouveaux variables pour distiguer la version et type d'ocean : type_ocean=force/slab/couple, version_ocean=opa8/nemo pour couplé ou version_ocean=sicOBS pour slab

JG

File:

• LMDZ4/trunk/libf/phylmd/physiq.F (modified) (13 diffs)

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -24 +24 @@
       USE vampir
       USE pbl_surface_mod, ONLY : pbl_surface
-      USE surface_data,     ONLY : ocean, ok_veget
+      USE change_srf_frac_mod
+      USE surface_data,     ONLY : type_ocean, ok_veget
       USE phys_local_var_mod ! Variables internes non sauvegardees de la physique
       USE phys_state_var_mod ! Variables sauvegardees de la physique
-
-
-      USE ocean_slab_mod, ONLY   : ocean_slab_get_vars
-      USE ocean_cpl_mod, ONLY    : ocean_cpl_get_vars
-      USE ocean_forced_mod, ONLY : ocean_forced_get_vars
       USE fonte_neige_mod, ONLY  : fonte_neige_get_vars
       USE phys_output_mod
@@ -120 +116 @@
       LOGICAL, SAVE :: rnpb=.TRUE.
 c$OMP THREADPRIVATE(rnpb)
-cIM "slab" ocean
-      REAL tslab(klon)    !Temperature du slab-ocean
-      REAL seaice(klon)   !glace de mer (kg/m2) 
-      REAL fluxo(klon)    !flux turbulents ocean-glace de mer 
-      REAL fluxg(klon)    !flux turbulents ocean-atmosphere
       REAL amn, amx
       INTEGER igout
@@ -699 +690 @@
       REAL wfbils(klon,nbsrf) ! bilan de chaleur au sol, pour chaque
 C                             ! type de sous-surface et pondere par la fraction
+      REAL slab_wfbils(klon)  ! bilan de chaleur au sol pour le cas de slab, sur les points d'ocean
+
       REAL fder(klon)         
       REAL ve(klon) ! integr. verticale du transport meri. de l'energie
@@ -715 +708 @@
       SAVE lmt_pas                ! frequence de mise a jour
 c$OMP THREADPRIVATE(lmt_pas) 
-cIM
-      REAL pctsrf_new(klon,nbsrf) !pourcentage surfaces issus d'ORCHIDEE
-
-cym      SAVE pctsrf                 ! sous-fraction du sol
 
 cIM sorties
@@ -1239 +1228 @@
 c appel a la lecture du run.def physique
 c
-         call conf_phys(ocean, ok_veget, ok_journe, ok_mensuel,
+         call conf_phys(ok_journe, ok_mensuel,
      .                  ok_instan, ok_hf,
      .                  solarlong0,seuil_inversion,
@@ -1279 +1268 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
-         CALL phyetat0 ("startphy.nc",ocean, ok_veget,clesphy0,tabcntr0)
+         CALL phyetat0 ("startphy.nc",clesphy0,tabcntr0)
 cIM begin
           print*,'physiq: clwcon rnebcon ratqs',clwcon(1,1),rnebcon(1,1)
@@ -1288 +1277 @@
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-
-
-       DO i=1,klon
-         IF ( abs( pctsrf(i, is_ter) + pctsrf(i, is_lic) + 
-     $       pctsrf(i, is_oce) + pctsrf(i, is_sic)  - 1.) .GT. EPSFRA) 
-     $       THEN 
-            WRITE(*,*)
-     $   'physiq apres lecture de restart: pb sous surface au point ',
-     $   i, pctsrf(i, 1 : nbsrf)
-         ENDIF 
-      ENDDO
-
 c
 C on remet le calendrier a zero
@@ -1518 +1495 @@
        call phys_output_open(jjmp1,nqmax,nlevSTD,clevSTD,nbteta, 
      &                        ctetaSTD,dtime,presnivs,ok_veget,
-     &                        ocean,iflag_pbl,ok_mensuel,ok_journe,
+     &                        type_ocean,iflag_pbl,ok_mensuel,ok_journe,
      &                        ok_hf,ok_instan,nid_files)
 c$OMP END MASTER
@@ -1581 +1558 @@
 c
       ENDIF
-c
-c   ****************     Fin  de   IF ( debut  )   ***************
-c
+!
+!   ****************     Fin  de   IF ( debut  )   ***************
+!
+!
+! Incrementer le compteur de la physique
+!
+      itap   = itap + 1
+      julien = MOD(NINT(xjour),360)
+      if (julien .eq. 0) julien = 360
+
+!
+! Update fraction of the sub-surfaces (pctsrf) and 
+! initialize, where a new fraction has appeared, all variables depending 
+! on the surface fraction.
+!
+      CALL change_srf_frac(itap, dtime, julien, 
+     *     pctsrf, falb1, falb2, ftsol, u10m, v10m, pbl_tke)
+
 
 ! Tendances bidons pour les processus qui n'affectent pas certaines
@@ -1725 +1717 @@
 cIM END
 c
-c Incrementer le compteur de la physique
-c
-      itap   = itap + 1
-      julien = MOD(NINT(xjour),360)
-      if (julien .eq. 0) julien = 360
-c
 c Mettre en action les conditions aux limites (albedo, sst, etc.).
 c Prescrire l'ozone et calculer l'albedo sur l'ocean.
@@ -1836 +1822 @@
      s     zxtsol,    zxfluxlat, zt2m,    qsat2m, 
      s     d_t_vdf,   d_q_vdf,   d_u_vdf, d_v_vdf,
-     s     ycoefh,    pctsrf_new,                
+     s     ycoefh,    slab_wfbils,                
      d     qsol,      zq2m,      s_pblh,  s_lcl,
      d     s_capCL,   s_oliqCL,  s_cteiCL,s_pblT,
@@ -1847 +1833 @@
      -     dsens,     devap,     zxsnow,
      -     zxfluxt,   zxfluxq,   q2m,     fluxq, pbl_tke )
-c
-c
-      pctsrf(:,:) = pctsrf_new(:,:)
+
       
 !-----------------------------------------------------------------------------------------
@@ -3332 +3316 @@
      .     zxfqcalving, zxfqfonte, zxffonte)
 
-      IF (ocean == 'slab') THEN
-         ! Get some variables from module ocean_slab_mod
-         CALL ocean_slab_get_vars(tslab, seaice, fluxo, fluxg)
-      ELSEIF (ocean == 'couple') THEN
-         ! Get some variables from module ocean_cpl_mod
-         CALL ocean_cpl_get_vars(fluxo, fluxg)
-      ELSE 
-         ! Get some variables from module ocean_forced_mod
-         CALL ocean_forced_get_vars(fluxo, fluxg)
-      ENDIF
-
 
 c Commente par abderrahmane le 11 2 08