Changeset 223

Timestamp:

Jun 6, 2001, 4:27:20 PM (23 years ago)

Author:

lmdzadmin

Message:

Impression debug dans clmain MAFO
Introduction du physiq.def LF
LF

Location:

LMDZ.3.3/branches/rel-LF/libf/phylmd

Files:

• clmain.F (modified) (3 diffs)
• conf_phys.F90 (modified) (3 diffs)
• interface_surf.F90 (modified) (8 diffs)
• physiq.F (modified) (7 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/clmain.F

@@ -28 +28 @@
 cAA il faudra sortir ces memes champs en leur ajoutant une dimension, 
 cAA c'est a dire nbsrf (nbre de subsurface).
- 
+      USE ioipsl
+      USE interface_surf
       IMPLICIT none
 c======================================================================
@@ -174 +175 @@
       REAL zx_alf1, zx_alf2 !valeur ambiante par extrapola.
 c======================================================================
-
+c
+c maf pour sorties IOISPL en cas de debugagage
+c
+      CHARACTER*80 cldebug
+      SAVE cldebug
+      CHARACTER*8 cl_surf(nbsrf)
+      SAVE cl_surf
+      INTEGER nhoridbg, nidbg
+      SAVE nhoridbg, nidbg
+      INTEGER ndexbg(iim*(jjm+1))
+      REAL zx_lon(iim,jjm+1), zx_lat(iim,jjm+1), zjulian
+      REAL tabindx(klon)
+      REAL debugtab(iim,jjm+1)
+      LOGICAL first_appel
+      SAVE first_appel
+      DATA first_appel/.true./
+      LOGICAL debugindex
+      SAVE debugindex
+      DATA debugindex/.true./
+#include "temps.h"
+      
+      IF (first_appel) THEN
+          first_appel=.false.
+!
+! initialisation sorties netcdf
+!
+          CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
+          zjulian = zjulian + day_ini
+          CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlon,zx_lon)
+          DO i = 1, iim
+            zx_lon(i,1) = rlon(i+1)
+            zx_lon(i,jjm+1) = rlon(i+1)
+          ENDDO
+          CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlat,zx_lat)
+          cldebug='sous_index'
+          CALL histbeg(cldebug, iim,zx_lon,jjm+1,zx_lat,1,iim,1,jjm
+     $        +1, 0,zjulian,dtime,nhoridbg,nidbg) 
+! no vertical axis
+          cl_surf(1)='ter'
+          cl_surf(2)='lic'
+          cl_surf(3)='oce'
+          cl_surf(4)='sic'
+          DO nsrf=1,nbsrf
+            CALL histdef(nidbg, cl_surf(nsrf),cl_surf(nsrf), "-",iim,
+     $          jjm+1,nhoridbg, 1, 1, 1, -99, 32, "inst", dtime,dtime) 
+          END DO
+          CALL histend(nidbg)
+          CALL histsync(nidbg)
+      ENDIF 
+          
       DO k = 1, klev   ! epaisseur de couche
       DO i = 1, klon
@@ -294 +344 @@
 c
       write(*,*)'CLMAIN, nsrf, knon =',nsrf, knon
+c
+c variables pour avoir une sortie IOIPSL des INDEX
+c
+      IF (debugindex) THEN 
+          tabindx(:)=0.
+c          tabindx(1:knon)=(/FLOAT(i),i=1:knon/)
+          DO i=1,knon
+            tabindx(1:knon)=FLOAT(i)
+          END DO 
+          debugtab(:,:)=0.
+          ndexbg(:)=0
+          CALL gath2cpl(tabindx,debugtab,klon,knon,iim,jjm,ni)
+          CALL histwrite(nidbg,cl_surf(nsrf),itap,debugtab,iim*(jjm+1)
+     $        ,ndexbg)
+      ENDIF 
       IF (knon.EQ.0) GOTO 99999
       DO j = 1, knon

LMDZ.3.3/branches/rel-LF/libf/phylmd/conf_phys.F90

@@ -22 +22 @@
 
   use IOIPSL     
-
+  implicit none
 !
 ! ocean:      type d'ocean (force, slab, couple)
@@ -36 +36 @@
   logical              :: ok_veget
   logical              :: ok_journe, ok_mensuel, ok_instan
+! Local
+  integer              :: numout = 6
 
 !
@@ -89 +91 @@
 !
 
-    write(numout,*)' Config ocean = ', ocean
-    write(numout,*)' Config veget = ', ok_veget
-    write(numout,*)' Sortie journaliere = ', ok_journe
-    write(numout,*)' Sortie mensuelle = ', ok_mensuel
-    write(numout,*)' Sortie instantanee = ', ok_instan
+  write(numout,*)' ##############################################'
+  write(numout,*)' Configuration des parametres de la physique: '
+  write(numout,*)' Config ocean = ', ocean
+  write(numout,*)' Config veget = ', ok_veget
+  write(numout,*)' Sortie journaliere = ', ok_journe
+  write(numout,*)' Sortie mensuelle = ', ok_mensuel
+  write(numout,*)' Sortie instantanee = ', ok_instan
 
 

LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90

@@ -1 @@
-! $Header
+! $Header$
 
   MODULE interface_surf
@@ -26 +26 @@
 
   PRIVATE
-  PUBLIC :: interfsurf,interfsurf_hq 
+  PUBLIC :: interfsurf,interfsurf_hq,gath2cpl 
 
   INTERFACE interfsurf
@@ -197 +197 @@
   real, DIMENSION(klon):: zfra
   logical              :: cumul = .false.
-  logical              :: scatter = .false.
+  logical,save         :: scatter = .false.
 
   if (check) write(*,*) 'Entree ', modname
@@ -345 +345 @@
 !  appel a sechiba
 !
-      scatter= .true.
+      scatter = .true.
+      call getin('SCATTER',scatter)
       if (.not. scatter) then
         call interfsol(itime, klon, dtime, date0, nisurf, knon, &
@@ -753 +754 @@
   real, allocatable, dimension (:,:), save :: lon_scat, lat_scat  
 
-  logical, save             :: lrestart_read = .true. , lrestart_write = .true.
+  logical, save          :: lrestart_read = .true. , lrestart_write = .false.
 
   real, dimension(klon):: qsurf
@@ -928 +929 @@
 ! Appel a la routine sols continentaux
 !
+  if (lafin) lrestart_write = .true.
+  if (check) write(*,*)'lafin ',lafin,lrestart_write
 
   petA_orc = petBcoef * dtime
@@ -960 +963 @@
 
   alb_new(:) = albedo_out(:,1)
+
+! LF essai sensible
+  fluxsens = -1. * fluxsens
+  fluxlat  = -1. * fluxlat
+
+  if (debut) lrestart_read = .false.
 
   END SUBROUTINE interfsol
@@ -2642 +2651 @@
      & * p1lay(i)/(RD*t1lay(i))
   ENDDO
+
+
 ! === Calcul de la temperature de surface ===
 ! 

LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

@@ -81 +81 @@
       EXTERNAL fromcpl, intocpl, inicma
 c      ocean = type de modele ocean a utiliser: force, slab, couple
-      character *6 ocean
-      parameter (ocean = 'force ')
+      character*6 ocean
+      SAVE ocean
+
+c      parameter (ocean = 'force ')
 c     parameter (ocean = 'couple')
       logical ok_ocean
@@ -94 +96 @@
 ccc      PARAMETER (soil_model=.FALSE.)
       logical ok_veget
-      parameter (ok_veget = .true.)
+      save ok_veget
+c     parameter (ok_veget = .true.)
 c      parameter (ok_veget = .false.)
 c======================================================================
@@ -112 +115 @@
 c======================================================================
       LOGICAL ok_journe ! sortir le fichier journalier
-      PARAMETER (ok_journe=.true.)
+      save ok_journe
+c      PARAMETER (ok_journe=.true.)
 c
       LOGICAL ok_mensuel ! sortir le fichier mensuel
-      PARAMETER (ok_mensuel=.true.)
+      save ok_mensuel
+c      PARAMETER (ok_mensuel=.true.)
 c
       LOGICAL ok_instan ! sortir le fichier instantane
-      PARAMETER (ok_instan=.false.)
+      save ok_instan
+c      PARAMETER (ok_instan=.true.)
 c
       LOGICAL ok_region ! sortir le fichier regional
@@ -527 +533 @@
        IF (debut) THEN
 
+c
+c appel a la lecture du run.def physique
+c
+         call conf_phys(ocean, ok_veget, ok_journe, ok_mensuel,
+     .                  ok_instan)
+
          DO k = 2, nvm          ! pas de vegetation
             DO i = 1, klon
@@ -786 +798 @@
 C
            call histdef(nid_day, "albe_"//clnsurf(nsrf), 
-     $         "Latent heat flux "//clnsurf(nsrf), "W/m2",  
+     $         "Albedo surf. "//clnsurf(nsrf), "W/m2",  
      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
      $         "ave(X)", zsto,zout)
@@ -1023 +1035 @@
 c
            call histdef(nid_mth, "albe_"//clnsurf(nsrf), 
-     $         "Latent heat flux "//clnsurf(nsrf), "W/m2",  
+     $         "Albedo surf. "//clnsurf(nsrf), "W/m2",  
      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
      $         "ave(X)", zsto,zout)
@@ -1371 +1383 @@
 c
            call histdef(nid_ins, "albe_"//clnsurf(nsrf), 
-     $         "Albedo "//clnsurf(nsrf), "-",  
+     $         "Albedo surf. "//clnsurf(nsrf), "-",  
      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
      $         "inst(X)", zsto,zout)