Index: LMDZ5/trunk/libf/dyn3d/conf_gcm.F =================================================================== --- LMDZ5/trunk/libf/dyn3d/conf_gcm.F (revision 2140) +++ LMDZ5/trunk/libf/dyn3d/conf_gcm.F (revision 2141) @@ -355,6 +355,4 @@ CALL getin('ok_guide',ok_guide) -! ............................................................... - !Config Key = read_start !Config Desc = Initialize model using a 'start.nc' file @@ -381,4 +379,7 @@ CALL getin('iphysiq', iphysiq) + if (mod(iphysiq, iperiod) /= 0) call abort_gcm("conf_gcm", + $ "iphysiq must be a multiple of iperiod", 1) + !Config Key = ip_ebil_dyn !Config Desc = PRINT level for energy conserv. diag. @@ -421,5 +422,5 @@ ! 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 @@ -436,5 +437,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 @@ -450,5 +451,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 @@ -456,5 +457,5 @@ IF( grossismx.LT.1. ) THEN - write(lunout,*) & + write(lunout,*) & 'conf_gcm: *** ATTENTION !! grossismx < 1 . *** ' STOP @@ -465,5 +466,5 @@ IF( grossismy.LT.1. ) THEN - write(lunout,*) & + write(lunout,*) & 'conf_gcm: *** ATTENTION !! grossismy < 1 . *** ' STOP @@ -489,6 +490,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 @@ -496,6 +497,6 @@ 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 @@ -512,6 +513,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 @@ -528,6 +529,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 @@ -543,6 +544,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 @@ -558,6 +559,6 @@ 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 @@ -579,6 +580,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 @@ -586,6 +587,6 @@ 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 @@ -725,5 +726,6 @@ IF( grossismx.LT.1. ) THEN - write(lunout,*)'conf_gcm: ***ATTENTION !! grossismx < 1 . *** ' + write(lunout,*) + & 'conf_gcm: *** ATTENTION !! grossismx < 1 . *** ' STOP ELSE @@ -869,6 +871,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 @@ -897,5 +899,5 @@ write(lunout,*)' #########################################' - write(lunout,*)' Configuration des parametres de cel0' & + write(lunout,*)' Configuration des parametres de cel0' & //'_limit: ' write(lunout,*)' planet_type = ', planet_type Index: LMDZ5/trunk/libf/dyn3dmem/conf_gcm.F =================================================================== --- LMDZ5/trunk/libf/dyn3dmem/conf_gcm.F (revision 2140) +++ LMDZ5/trunk/libf/dyn3dmem/conf_gcm.F (revision 2141) @@ -2,8 +2,8 @@ ! $Id$ ! -c -c +! +! SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 ) -c +! USE control_mod #ifdef CPP_IOIPSL @@ -21,14 +21,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 @@ -37,7 +37,7 @@ PARAMETER( longcles = 20 ) REAL clesphy0( longcles ) -c -c Declarations : -c -------------- +! +! Declarations : +! -------------- #include "dimensions.h" #include "paramet.h" @@ -51,8 +51,8 @@ ! #include "clesphys.h" #include "iniprint.h" -c -c -c local: -c ------ +! +! +! local: +! ------ CHARACTER ch1*72,ch2*72,ch3*72,ch4*12 @@ -66,33 +66,33 @@ integer,external :: OMP_GET_NUM_THREADS #endif -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 !Config Desc = unite de fichier pour les impressions @@ -134,7 +134,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", ...) @@ -266,6 +266,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 @@ -409,4 +409,7 @@ CALL getin('iphysiq', iphysiq) + if (mod(iphysiq, iperiod) /= 0) call abort_gcm("conf_gcm", + $ "iphysiq must be a multiple of iperiod", 1) + !Config Key = ip_ebil_dyn !Config Desc = PRINT level for energy conserv. diag. @@ -422,7 +425,7 @@ -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 @@ -443,6 +446,6 @@ 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', @@ -498,8 +501,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 @@ -526,5 +529,5 @@ ENDIF ENDIF -c +! !Config Key = dzoomx !Config Desc = extension en longitude @@ -589,5 +592,5 @@ ENDIF -cc +!c IF( .NOT.fxyhypb ) THEN @@ -617,5 +620,5 @@ ENDIF ENDIF ! of IF( .NOT.fxyhypb ) -c +! !Config Key = offline !Config Desc = Nouvelle eau liquide @@ -722,6 +725,6 @@ RETURN -c ............................................... -c +! ............................................... +! 100 CONTINUE !Config Key = clon @@ -767,6 +770,5 @@ IF( grossismy.LT.1. ) THEN - write(lunout,*) - & 'conf_gcm: *** ATTENTION !! grossismy < 1 . *** ' + write(lunout,*) 'conf_gcm: ***ATTENTION !! grossismy < 1 . *** ' STOP ELSE @@ -775,8 +777,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 @@ -826,5 +828,5 @@ ysinus = .TRUE. CALL getin('ysinus',ysinus) -c +! !Config Key = offline !Config Desc = Nouvelle eau liquide @@ -891,5 +893,5 @@ !Config Def = false !Config Help = permet d'activer l'utilisation des FFT pour effectuer -!Config le filtrage aux poles. +!Config le filtrage aux poles. use_filtre_fft=.FALSE. CALL getin('use_filtre_fft',use_filtre_fft) @@ -1001,5 +1003,5 @@ write(lunout,*)' ok_limit = ', ok_limit write(lunout,*)' ok_etat0 = ', ok_etat0 -c +! RETURN END Index: LMDZ5/trunk/libf/dyn3dpar/conf_gcm.F =================================================================== --- LMDZ5/trunk/libf/dyn3dpar/conf_gcm.F (revision 2140) +++ LMDZ5/trunk/libf/dyn3dpar/conf_gcm.F (revision 2141) @@ -2,8 +2,8 @@ ! $Id$ ! -c -c +! +! SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 ) -c +! USE control_mod #ifdef CPP_IOIPSL @@ -20,14 +20,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 @@ -36,7 +36,7 @@ PARAMETER( longcles = 20 ) REAL clesphy0( longcles ) -c -c Declarations : -c -------------- +! +! Declarations : +! -------------- #include "dimensions.h" #include "paramet.h" @@ -50,8 +50,8 @@ ! #include "clesphys.h" #include "iniprint.h" -c -c -c local: -c ------ +! +! +! local: +! ------ CHARACTER ch1*72,ch2*72,ch3*72,ch4*12 @@ -65,33 +65,33 @@ integer,external :: OMP_GET_NUM_THREADS #endif -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 !Config Desc = unite de fichier pour les impressions @@ -133,7 +133,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", ...) @@ -265,6 +265,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 @@ -408,4 +408,7 @@ CALL getin('iphysiq', iphysiq) + if (mod(iphysiq, iperiod) /= 0) call abort_gcm("conf_gcm", + $ "iphysiq must be a multiple of iperiod", 1) + !Config Key = ip_ebil_dyn !Config Desc = PRINT level for energy conserv. diag. @@ -421,7 +424,7 @@ -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 @@ -442,6 +445,6 @@ 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', @@ -497,8 +500,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 @@ -525,5 +528,5 @@ ENDIF ENDIF -c +! !Config Key = dzoomx !Config Desc = extension en longitude @@ -588,5 +591,5 @@ ENDIF -cc +!c IF( .NOT.fxyhypb ) THEN @@ -616,5 +619,5 @@ ENDIF ENDIF ! of IF( .NOT.fxyhypb ) -c +! !Config Key = offline !Config Desc = Nouvelle eau liquide @@ -721,6 +724,6 @@ RETURN -c ............................................... -c +! ............................................... +! 100 CONTINUE !Config Key = clon @@ -766,6 +769,5 @@ IF( grossismy.LT.1. ) THEN - write(lunout,*) - & 'conf_gcm: *** ATTENTION !! grossismy < 1 . *** ' + write(lunout,*) 'conf_gcm: ***ATTENTION !! grossismy < 1 . *** ' STOP ELSE @@ -774,8 +776,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 @@ -825,5 +827,5 @@ ysinus = .TRUE. CALL getin('ysinus',ysinus) -c +! !Config Key = offline !Config Desc = Nouvelle eau liquide @@ -890,5 +892,5 @@ !Config Def = false !Config Help = permet d'activer l'utilisation des FFT pour effectuer -!Config le filtrage aux poles. +!Config le filtrage aux poles. use_filtre_fft=.FALSE. CALL getin('use_filtre_fft',use_filtre_fft) @@ -999,5 +1001,5 @@ write(lunout,*)' ok_limit = ', ok_limit write(lunout,*)' ok_etat0 = ', ok_etat0 -c +! RETURN END