Index: LMDZ5/branches/testing/libf/dyn3d/calfis.F =================================================================== --- LMDZ5/branches/testing/libf/dyn3d/calfis.F (revision 2024) +++ LMDZ5/branches/testing/libf/dyn3d/calfis.F (revision 2056) @@ -163,9 +163,4 @@ REAL unskap, pksurcp c -cIM diagnostique PVteta, Amip2 - INTEGER,PARAMETER :: ntetaSTD=3 - REAL,SAVE :: rtetaSTD(ntetaSTD)=(/350.,380.,405./) ! Earth-specific, beware !! - REAL PVteta(ngridmx,ntetaSTD) -c REAL flxwfi(ngridmx,llm) ! Flux de masse verticale sur la grille physiq c @@ -431,14 +426,4 @@ ENDDO -c - if (planet_type=="earth") then -#ifdef CPP_PHYS -! PVtheta calls tetalevel, which is in the physics -cIM calcul PV a teta=350, 380, 405K - CALL PVtheta(ngridmx,llm,pucov,pvcov,pteta, - $ ztfi,zplay,zplev, - $ ntetaSTD,rtetaSTD,PVteta) -#endif - endif c c On change de grille, dynamique vers physiq, pour le flux de masse verticale @@ -491,7 +476,5 @@ . zdqfi, . zdpsrf, -cIM diagnostique PVteta, Amip2 - . pducov, - . PVteta) + . pducov) else if ( planet_type=="generic" ) then Index: LMDZ5/branches/testing/libf/dyn3d/conf_gcm.F =================================================================== --- LMDZ5/branches/testing/libf/dyn3d/conf_gcm.F (revision 2024) +++ LMDZ5/branches/testing/libf/dyn3d/conf_gcm.F (revision 2056) @@ -2,8 +2,8 @@ ! $Id$ ! -c -c +! +! SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 ) -c +! USE control_mod #ifdef CPP_IOIPSL @@ -17,14 +17,14 @@ IMPLICIT NONE -c----------------------------------------------------------------------- -c Auteurs : L. Fairhead , P. Le Van . -c -c Arguments : -c -c tapedef : -c etatinit : = TRUE , on ne compare pas les valeurs des para- -c -metres du zoom avec celles lues sur le fichier start . -c clesphy0 : sortie . -c +!----------------------------------------------------------------------- +! Auteurs : L. Fairhead , P. Le Van . +! +! Arguments : +! +! tapedef : +! etatinit : = TRUE , on ne compare pas les valeurs des para- +! -metres du zoom avec celles lues sur le fichier start . +! clesphy0 : sortie . +! LOGICAL etatinit INTEGER tapedef @@ -33,7 +33,7 @@ PARAMETER( longcles = 20 ) REAL clesphy0( longcles ) -c -c Declarations : -c -------------- +! +! Declarations : +! -------------- #include "dimensions.h" #include "paramet.h" @@ -47,8 +47,8 @@ ! #include "clesphys.h" #include "iniprint.h" -c -c -c local: -c ------ +! +! +! local: +! ------ CHARACTER ch1*72,ch2*72,ch3*72,ch4*12 @@ -58,32 +58,32 @@ INTEGER i LOGICAL use_filtre_fft -c -c ------------------------------------------------------------------- -c -c ......... Version du 29/04/97 .......... -c -c Nouveaux parametres nitergdiv,nitergrot,niterh,tetagdiv,tetagrot, -c tetatemp ajoutes pour la dissipation . -c -c Autre parametre ajoute en fin de liste de tapedef : ** fxyhypb ** -c -c Si fxyhypb = .TRUE. , choix de la fonction a derivee tangente hyperb. -c Sinon , choix de fxynew , a derivee sinusoidale .. -c -c ...... etatinit = . TRUE. si defrun est appele dans ETAT0_LMD ou -c LIMIT_LMD pour l'initialisation de start.dat (dic) et -c de limit.dat ( dic) ........... -c Sinon etatinit = . FALSE . -c -c Donc etatinit = .F. si on veut comparer les valeurs de grossismx , -c grossismy,clon,clat, fxyhypb lues sur le fichier start avec -c celles passees par run.def , au debut du gcm, apres l'appel a -c lectba . -c Ces parmetres definissant entre autres la grille et doivent etre -c pareils et coherents , sinon il y aura divergence du gcm . -c -c----------------------------------------------------------------------- -c initialisations: -c ---------------- +! +! ------------------------------------------------------------------- +! +! ......... Version du 29/04/97 .......... +! +! Nouveaux parametres nitergdiv,nitergrot,niterh,tetagdiv,tetagrot, +! tetatemp ajoutes pour la dissipation . +! +! Autre parametre ajoute en fin de liste de tapedef : ** fxyhypb ** +! +! Si fxyhypb = .TRUE. , choix de la fonction a derivee tangente hyperb. +! Sinon , choix de fxynew , a derivee sinusoidale .. +! +! ...... etatinit = . TRUE. si defrun est appele dans ETAT0_LMD ou +! LIMIT_LMD pour l'initialisation de start.dat (dic) et +! de limit.dat ( dic) ........... +! Sinon etatinit = . FALSE . +! +! Donc etatinit = .F. si on veut comparer les valeurs de grossismx , +! grossismy,clon,clat, fxyhypb lues sur le fichier start avec +! celles passees par run.def , au debut du gcm, apres l'appel a +! lectba . +! Ces parmetres definissant entre autres la grille et doivent etre +! pareils et coherents , sinon il y aura divergence du gcm . +! +!----------------------------------------------------------------------- +! initialisations: +! ---------------- !Config Key = lunout @@ -95,5 +95,5 @@ CALL getin('lunout', lunout) IF (lunout /= 5 .and. lunout /= 6) THEN - OPEN(UNIT=lunout,FILE='lmdz.out',ACTION='write', + OPEN(UNIT=lunout,FILE='lmdz.out',ACTION='write', & & STATUS='unknown',FORM='formatted') ENDIF @@ -107,7 +107,7 @@ CALL getin('prt_level',prt_level) -c----------------------------------------------------------------------- -c Parametres de controle du run: -c----------------------------------------------------------------------- +!----------------------------------------------------------------------- +! Parametres de controle du run: +!----------------------------------------------------------------------- !Config Key = planet_type !Config Desc = planet type ("earth", "mars", "venus", ...) @@ -232,6 +232,6 @@ CALL getin('dissip_period',dissip_period) -ccc .... P. Le Van , modif le 29/04/97 .pour la dissipation ... -ccc +!cc .... P. Le Van , modif le 29/04/97 .pour la dissipation ... +!cc !Config Key = lstardis @@ -348,5 +348,5 @@ CALL getin('ok_guide',ok_guide) -c ............................................................... +! ............................................................... !Config Key = read_start @@ -390,7 +390,7 @@ ENDDO -ccc .... P. Le Van , ajout le 7/03/95 .pour le zoom ... -c ......... ( modif le 17/04/96 ) ......... -c +!cc .... P. Le Van , ajout le 7/03/95 .pour le zoom ... +! ......... ( modif le 17/04/96 ) ......... +! IF( etatinit ) GO TO 100 @@ -411,8 +411,8 @@ CALL getin('clat',clatt) -c -c +! +! IF( ABS(clat - clatt).GE. 0.001 ) THEN - write(lunout,*)'conf_gcm: La valeur de clat passee par run.def', + write(lunout,*)'conf_gcm: La valeur de clat passee par run.def', & & ' est differente de celle lue sur le fichier start ' STOP @@ -429,5 +429,5 @@ IF( ABS(grossismx - grossismxx).GE. 0.001 ) THEN - write(lunout,*)'conf_gcm: La valeur de grossismx passee par ', + write(lunout,*)'conf_gcm: La valeur de grossismx passee par ', & & 'run.def est differente de celle lue sur le fichier start ' STOP @@ -443,5 +443,5 @@ IF( ABS(grossismy - grossismyy).GE. 0.001 ) THEN - write(lunout,*)'conf_gcm: La valeur de grossismy passee par ', + write(lunout,*)'conf_gcm: La valeur de grossismy passee par ', & & 'run.def est differente de celle lue sur le fichier start ' STOP @@ -449,5 +449,5 @@ IF( grossismx.LT.1. ) THEN - write(lunout,*) + write(lunout,*) & & 'conf_gcm: *** ATTENTION !! grossismx < 1 . *** ' STOP @@ -458,5 +458,5 @@ IF( grossismy.LT.1. ) THEN - write(lunout,*) + write(lunout,*) & & 'conf_gcm: *** ATTENTION !! grossismy < 1 . *** ' STOP @@ -466,8 +466,8 @@ write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay -c -c alphax et alphay sont les anciennes formulat. des grossissements -c -c +! +! alphax et alphay sont les anciennes formulat. des grossissements +! +! !Config Key = fxyhypb @@ -482,6 +482,6 @@ IF( fxyhypbb ) THEN write(lunout,*)' ******** PBS DANS CONF_GCM ******** ' - write(lunout,*)' *** fxyhypb lu sur le fichier start est ', - * 'F alors qu il est T sur run.def ***' + write(lunout,*)' *** fxyhypb lu sur le fichier start est ', & + & 'F alors qu il est T sur run.def ***' STOP ENDIF @@ -489,10 +489,10 @@ IF( .NOT.fxyhypbb ) THEN write(lunout,*)' ******** PBS DANS CONF_GCM ******** ' - write(lunout,*)' *** fxyhypb lu sur le fichier start est ', - * 'T alors qu il est F sur run.def **** ' + write(lunout,*)' *** fxyhypb lu sur le fichier start est ', & + & 'T alors qu il est F sur run.def **** ' STOP ENDIF ENDIF -c +! !Config Key = dzoomx !Config Desc = extension en longitude @@ -505,6 +505,6 @@ IF( fxyhypb ) THEN IF( ABS(dzoomx - dzoomxx).GE. 0.001 ) THEN - write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ', - * 'run.def est differente de celle lue sur le fichier start ' + write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ', & + & 'run.def est differente de celle lue sur le fichier start ' STOP ENDIF @@ -521,6 +521,6 @@ IF( fxyhypb ) THEN IF( ABS(dzoomy - dzoomyy).GE. 0.001 ) THEN - write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ', - * 'run.def est differente de celle lue sur le fichier start ' + write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ', & + & 'run.def est differente de celle lue sur le fichier start ' STOP ENDIF @@ -536,6 +536,6 @@ IF( fxyhypb ) THEN IF( ABS(taux - tauxx).GE. 0.001 ) THEN - write(lunout,*)'conf_gcm: La valeur de taux passee par ', - * 'run.def est differente de celle lue sur le fichier start ' + write(lunout,*)'conf_gcm: La valeur de taux passee par ', & + & 'run.def est differente de celle lue sur le fichier start ' STOP ENDIF @@ -551,11 +551,11 @@ IF( fxyhypb ) THEN IF( ABS(tauy - tauyy).GE. 0.001 ) THEN - write(lunout,*)'conf_gcm: La valeur de tauy passee par ', - * 'run.def est differente de celle lue sur le fichier start ' + write(lunout,*)'conf_gcm: La valeur de tauy passee par ', & + & 'run.def est differente de celle lue sur le fichier start ' STOP ENDIF ENDIF -cc +!c IF( .NOT.fxyhypb ) THEN @@ -572,6 +572,6 @@ IF( ysinuss ) THEN write(lunout,*)' ******** PBS DANS CONF_GCM ******** ' - write(lunout,*)' *** ysinus lu sur le fichier start est F', - * ' alors qu il est T sur run.def ***' + write(lunout,*)' *** ysinus lu sur le fichier start est F', & + & ' alors qu il est T sur run.def ***' STOP ENDIF @@ -579,11 +579,11 @@ IF( .NOT.ysinuss ) THEN write(lunout,*)' ******** PBS DANS CONF_GCM ******** ' - write(lunout,*)' *** ysinus lu sur le fichier start est T', - * ' alors qu il est F sur run.def **** ' + write(lunout,*)' *** ysinus lu sur le fichier start est T', & + & ' alors qu il est F sur run.def **** ' STOP ENDIF ENDIF ENDIF ! of IF( .NOT.fxyhypb ) -c +! !Config Key = offline !Config Desc = Nouvelle eau liquide @@ -682,6 +682,6 @@ RETURN -c ............................................... -c +! ............................................... +! 100 CONTINUE !Config Key = clon @@ -718,6 +718,5 @@ IF( grossismx.LT.1. ) THEN - write(lunout,*) - & 'conf_gcm: *** ATTENTION !! grossismx < 1 . *** ' + write(lunout,*)'conf_gcm: ***ATTENTION !! grossismx < 1 . *** ' STOP ELSE @@ -727,6 +726,5 @@ IF( grossismy.LT.1. ) THEN - write(lunout,*) - & 'conf_gcm: *** ATTENTION !! grossismy < 1 . *** ' + write(lunout,*) 'conf_gcm: ***ATTENTION !! grossismy < 1 . *** ' STOP ELSE @@ -735,8 +733,8 @@ write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay -c -c alphax et alphay sont les anciennes formulat. des grossissements -c -c +! +! alphax et alphay sont les anciennes formulat. des grossissements +! +! !Config Key = fxyhypb @@ -786,5 +784,5 @@ ysinus = .TRUE. CALL getin('ysinus',ysinus) -c +! !Config Key = offline !Config Desc = Nouvelle eau liquide @@ -864,6 +862,6 @@ vert_prof_dissip = merge(1, 0, ok_strato .and. llm==39) CALL getin('vert_prof_dissip', vert_prof_dissip) - call assert(vert_prof_dissip == 0 .or. vert_prof_dissip == 1, - $ "bad value for vert_prof_dissip") + call assert(vert_prof_dissip == 0 .or. vert_prof_dissip == 1, & + & "bad value for vert_prof_dissip") !Config Key = ok_gradsfile @@ -892,5 +890,5 @@ write(lunout,*)' #########################################' - write(lunout,*)' Configuration des parametres de cel0' + write(lunout,*)' Configuration des parametres de cel0' & & //'_limit: ' write(lunout,*)' planet_type = ', planet_type @@ -937,5 +935,5 @@ write(lunout,*)' ok_limit = ', ok_limit write(lunout,*)' ok_etat0 = ', ok_etat0 -c +! RETURN END Index: LMDZ5/branches/testing/libf/dyn3d/gcm.F =================================================================== --- LMDZ5/branches/testing/libf/dyn3d/gcm.F (revision 2024) +++ LMDZ5/branches/testing/libf/dyn3d/gcm.F (revision 2056) @@ -105,7 +105,4 @@ REAL ps(ip1jmp1) ! pression au sol REAL p (ip1jmp1,llmp1 ) ! pression aux interfac.des couches - REAL pks(ip1jmp1) ! exner au sol - REAL pk(ip1jmp1,llm) ! exner au milieu des couches - REAL pkf(ip1jmp1,llm) ! exner filt.au milieu des couches REAL masse(ip1jmp1,llm) ! masse d'air REAL phis(ip1jmp1) ! geopotentiel au sol @@ -131,5 +128,4 @@ data call_iniphys/.true./ - REAL alpha(ip1jmp1,llm),beta(ip1jmp1,llm) c+jld variables test conservation energie c REAL ecin(ip1jmp1,llm),ecin0(ip1jmp1,llm) @@ -466,5 +462,9 @@ - day_end = day_ini + nday + if (nday>=0) then + day_end = day_ini + nday + else + day_end = day_ini - nday/day_step + endif WRITE(lunout,300)day_ini,day_end 300 FORMAT('1'/,15x,'run du jour',i7,2x,'au jour',i7//) Index: LMDZ5/branches/testing/libf/dyn3d/guide_mod.F90 =================================================================== --- LMDZ5/branches/testing/libf/dyn3d/guide_mod.F90 (revision 2024) +++ LMDZ5/branches/testing/libf/dyn3d/guide_mod.F90 (revision 2056) @@ -593,4 +593,6 @@ SUBROUTINE guide_interp(psi,teta) + use exner_hyb_m, only: exner_hyb + use exner_milieu_m, only: exner_milieu IMPLICIT NONE @@ -614,6 +616,5 @@ REAL, DIMENSION (iip1,jjm,llm) :: pbary ! Variables pour fonction Exner (P milieu couche) - REAL, DIMENSION (iip1,jjp1,llm) :: pk, pkf - REAL, DIMENSION (iip1,jjp1,llm) :: alpha, beta + REAL, DIMENSION (iip1,jjp1,llm) :: pk REAL, DIMENSION (iip1,jjp1) :: pks REAL :: prefkap,unskap @@ -680,7 +681,7 @@ CALL pression( ip1jmp1, ap, bp, psi, p ) if (pressure_exner) then - CALL exner_hyb(ip1jmp1,psi,p,alpha,beta,pks,pk,pkf) + CALL exner_hyb(ip1jmp1,psi,p,pks,pk) else - CALL exner_milieu(ip1jmp1,psi,p,beta,pks,pk,pkf) + CALL exner_milieu(ip1jmp1,psi,p,pks,pk) endif ! .... Calcul de pls , pression au milieu des couches ,en Pascals Index: LMDZ5/branches/testing/libf/dyn3d/iniacademic.F90 =================================================================== --- LMDZ5/branches/testing/libf/dyn3d/iniacademic.F90 (revision 2024) +++ LMDZ5/branches/testing/libf/dyn3d/iniacademic.F90 (revision 2056) @@ -14,4 +14,6 @@ #endif USE Write_Field + use exner_hyb_m, only: exner_hyb + use exner_milieu_m, only: exner_milieu ! Author: Frederic Hourdin original: 15/01/93 @@ -54,5 +56,4 @@ REAL pks(ip1jmp1) ! exner au sol REAL pk(ip1jmp1,llm) ! exner au milieu des couches - REAL pkf(ip1jmp1,llm) ! exner filt.au milieu des couches REAL phi(ip1jmp1,llm) ! geopotentiel REAL ddsin,zsig,tetapv,w_pv ! variables auxiliaires @@ -70,5 +71,5 @@ integer idum - REAL alpha(ip1jmp1,llm),beta(ip1jmp1,llm),zdtvr + REAL zdtvr character(len=*),parameter :: modname="iniacademic" @@ -223,7 +224,7 @@ CALL pression ( ip1jmp1, ap, bp, ps, p ) if (pressure_exner) then - CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf ) - else - call exner_milieu(ip1jmp1,ps,p,beta,pks,pk,pkf) + CALL exner_hyb( ip1jmp1, ps, p, pks, pk) + else + call exner_milieu(ip1jmp1,ps,p,pks,pk) endif CALL massdair(p,masse) Index: LMDZ5/branches/testing/libf/dyn3d/leapfrog.F =================================================================== --- LMDZ5/branches/testing/libf/dyn3d/leapfrog.F (revision 2024) +++ LMDZ5/branches/testing/libf/dyn3d/leapfrog.F (revision 2056) @@ -19,4 +19,7 @@ & iecri, ip_ebil_dyn, ok_dynzon, ok_dyn_ins, & periodav, ok_dyn_ave, output_grads_dyn + use exner_hyb_m, only: exner_hyb + use exner_milieu_m, only: exner_milieu + IMPLICIT NONE @@ -158,5 +161,4 @@ character*10 string10 - REAL alpha(ip1jmp1,llm),beta(ip1jmp1,llm) REAL :: flxw(ip1jmp1,llm) ! flux de masse verticale @@ -196,5 +198,9 @@ - itaufin = nday*day_step + if (nday>=0) then + itaufin = nday*day_step + else + itaufin = -nday + endif itaufinp1 = itaufin +1 itau = 0 @@ -217,7 +223,7 @@ CALL pression ( ip1jmp1, ap, bp, ps, p ) if (pressure_exner) then - CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf ) + CALL exner_hyb( ip1jmp1, ps, p, pks, pk, pkf ) else - CALL exner_milieu( ip1jmp1, ps, p, beta, pks, pk, pkf ) + CALL exner_milieu( ip1jmp1, ps, p, pks, pk, pkf ) endif @@ -373,8 +379,12 @@ CALL pression ( ip1jmp1, ap, bp, ps, p ) if (pressure_exner) then - CALL exner_hyb( ip1jmp1, ps, p,alpha,beta,pks, pk, pkf ) + CALL exner_hyb( ip1jmp1, ps, p,pks, pk, pkf ) else - CALL exner_milieu( ip1jmp1, ps, p, beta, pks, pk, pkf ) + CALL exner_milieu( ip1jmp1, ps, p, pks, pk, pkf ) endif + +! Appel a geopot ajoute le 2014/05/08 pour garantir la convergence numerique +! avec dyn3dmem + CALL geopot ( ip1jmp1, teta , pk , pks, phis , phi ) ! rdaym_ini = itau * dtvr / daysec @@ -448,7 +458,7 @@ CALL massdair(p,masse) if (pressure_exner) then - CALL exner_hyb(ip1jmp1,ps,p,alpha,beta,pks,pk,pkf) + CALL exner_hyb(ip1jmp1,ps,p,pks,pk,pkf) else - CALL exner_milieu(ip1jmp1,ps,p,beta,pks,pk,pkf) + CALL exner_milieu(ip1jmp1,ps,p,pks,pk,pkf) endif @@ -506,7 +516,7 @@ CALL pression ( ip1jmp1, ap, bp, ps, p ) if (pressure_exner) then - CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf ) + CALL exner_hyb( ip1jmp1, ps, p, pks, pk, pkf ) else - CALL exner_milieu( ip1jmp1, ps, p, beta, pks, pk, pkf ) + CALL exner_milieu( ip1jmp1, ps, p, pks, pk, pkf ) endif CALL massdair(p,masse)