Index: trunk/LMDZ.COMMON/libf/dyn3d/conf_gcm.F =================================================================== --- trunk/LMDZ.COMMON/libf/dyn3d/conf_gcm.F (revision 500) +++ trunk/LMDZ.COMMON/libf/dyn3d/conf_gcm.F (revision 617) @@ -34,4 +34,5 @@ #include "serre.h" #include "comdissnew.h" +#include "iniprint.h" #include "temps.h" #include "comconst.h" @@ -39,5 +40,4 @@ ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique ! #include "clesphys.h" -#include "iniprint.h" c c @@ -410,8 +410,223 @@ CALL getin('ip_ebil_dyn',ip_ebil_dyn) +!Config Key = offline +!Config Desc = Nouvelle eau liquide +!Config Def = n +!Config Help = Permet de mettre en route la +!Config nouvelle parametrisation de l'eau liquide ! + offline = .FALSE. + CALL getin('offline',offline) + +!Config Key = type_trac +!Config Desc = Choix de couplage avec model de chimie INCA ou REPROBUS +!Config Def = lmdz +!Config Help = +!Config 'lmdz' = pas de couplage, pur LMDZ +!Config 'inca' = model de chime INCA +!Config 'repr' = model de chime REPROBUS + type_trac = 'lmdz' + CALL getin('type_trac',type_trac) + +!Config Key = config_inca +!Config Desc = Choix de configuration de INCA +!Config Def = none +!Config Help = Choix de configuration de INCA : +!Config 'none' = sans INCA +!Config 'chem' = INCA avec calcul de chemie +!Config 'aero' = INCA avec calcul des aerosols + config_inca = 'none' + CALL getin('config_inca',config_inca) + +!Config Key = ok_dynzon +!Config Desc = calcul et sortie des transports +!Config Def = n +!Config Help = Permet de mettre en route le calcul des transports +!Config + ok_dynzon = .FALSE. + CALL getin('ok_dynzon',ok_dynzon) + +!Config Key = ok_dyn_ins +!Config Desc = sorties instantanees dans la dynamique +!Config Def = n +!Config Help = +!Config + ok_dyn_ins = .FALSE. + CALL getin('ok_dyn_ins',ok_dyn_ins) + +!Config Key = ok_dyn_ave +!Config Desc = sorties moyennes dans la dynamique +!Config Def = n +!Config Help = +!Config + ok_dyn_ave = .FALSE. + CALL getin('ok_dyn_ave',ok_dyn_ave) + +!Config Key = use_filtre_fft +!Config Desc = flag d'activation des FFT pour le filtre +!Config Def = false +!Config Help = permet d'activer l'utilisation des FFT pour effectuer +!Config le filtrage aux poles. +! Le filtre fft n'est pas implemente dans dyn3d + use_filtre_fft=.FALSE. + CALL getin('use_filtre_fft',use_filtre_fft) + + IF (use_filtre_fft) THEN + write(lunout,*)'STOP !!!' + write(lunout,*)'use_filtre_fft n est pas implemente dans dyn3d' + STOP + ENDIF + +!Config key = ok_strato +!Config Desc = activation de la version strato +!Config Def = .FALSE. +!Config Help = active la version stratosphérique de LMDZ de F. Lott + + ok_strato=.TRUE. + CALL getin('ok_strato',ok_strato) + +!Config Key = ok_gradsfile +!Config Desc = activation des sorties grads du guidage +!Config Def = n +!Config Help = active les sorties grads du guidage + + ok_gradsfile = .FALSE. + CALL getin('ok_gradsfile',ok_gradsfile) + +!Config Key = ok_limit +!Config Desc = creation des fichiers limit dans create_etat0_limit +!Config Def = y +!Config Help = production du fichier limit.nc requise + + ok_limit = .TRUE. + CALL getin('ok_limit',ok_limit) + +!Config Key = ok_etat0 +!Config Desc = creation des fichiers etat0 dans create_etat0_limit +!Config Def = y +!Config Help = production des fichiers start.nc, startphy.nc requise + + ok_etat0 = .TRUE. + CALL getin('ok_etat0',ok_etat0) + +!Config Key = grilles_gcm_netcdf +!Config Desc = creation de fichier grilles_gcm.nc dans create_etat0_limit +!Config Def = n + grilles_gcm_netcdf = .FALSE. + CALL getin('grilles_gcm_netcdf',grilles_gcm_netcdf) + +c---------------------------------------- ccc .... P. Le Van , ajout le 7/03/95 .pour le zoom ... c ......... ( modif le 17/04/96 ) ......... c - IF( etatinit ) GO TO 100 +C ZOOM PARAMETERS ... the ones read in start.nc prevail anyway ! (SL, 2012) +c +c---------------------------------------- + IF( etatinit ) then + +!Config Key = clon +!Config Desc = centre du zoom, longitude +!Config Def = 0 +!Config Help = longitude en degres du centre +!Config du zoom + clon = 0. + CALL getin('clon',clon) + +!Config Key = clat +!Config Desc = centre du zoom, latitude +!Config Def = 0 +!Config Help = latitude en degres du centre du zoom +!Config + clat = 0. + CALL getin('clat',clat) + +!Config Key = grossismx +!Config Desc = zoom en longitude +!Config Def = 1.0 +!Config Help = facteur de grossissement du zoom, +!Config selon la longitude + grossismx = 1.0 + CALL getin('grossismx',grossismx) + +!Config Key = grossismy +!Config Desc = zoom en latitude +!Config Def = 1.0 +!Config Help = facteur de grossissement du zoom, +!Config selon la latitude + grossismy = 1.0 + CALL getin('grossismy',grossismy) + + IF( grossismx.LT.1. ) THEN + write(lunout,*) + & 'conf_gcm: *** ATTENTION !! grossismx < 1 . *** ' + STOP + ELSE + alphax = 1. - 1./ grossismx + ENDIF + + + IF( grossismy.LT.1. ) THEN + write(lunout,*) + & 'conf_gcm: *** ATTENTION !! grossismy < 1 . *** ' + STOP + ELSE + alphay = 1. - 1./ grossismy + ENDIF + + write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay +c +c alphax et alphay sont les anciennes formulat. des grossissements +c +c + +!Config Key = fxyhypb +!Config Desc = Fonction hyperbolique +!Config Def = y +!Config Help = Fonction f(y) hyperbolique si = .true. +!Config sinon sinusoidale + fxyhypb = .TRUE. + CALL getin('fxyhypb',fxyhypb) + +!Config Key = dzoomx +!Config Desc = extension en longitude +!Config Def = 0 +!Config Help = extension en longitude de la zone du zoom +!Config ( fraction de la zone totale) + dzoomx = 0.0 + CALL getin('dzoomx',dzoomx) + +!Config Key = dzoomy +!Config Desc = extension en latitude +!Config Def = 0 +!Config Help = extension en latitude de la zone du zoom +!Config ( fraction de la zone totale) + dzoomy = 0.0 + CALL getin('dzoomy',dzoomy) + +!Config Key = taux +!Config Desc = raideur du zoom en X +!Config Def = 3 +!Config Help = raideur du zoom en X + taux = 3.0 + CALL getin('taux',taux) + +!Config Key = tauy +!Config Desc = raideur du zoom en Y +!Config Def = 3 +!Config Help = raideur du zoom en Y + tauy = 3.0 + CALL getin('tauy',tauy) + +!Config Key = ysinus +!Config IF = !fxyhypb +!Config Desc = Fonction en Sinus +!Config Def = y +!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true. +!Config sinon y = latit. + ysinus = .TRUE. + CALL getin('ysinus',ysinus) +c +c---------------------------------------- + else ! etatinit=false +c---------------------------------------- !Config Key = clon @@ -605,319 +820,11 @@ ENDIF ENDIF ! of IF( .NOT.fxyhypb ) -c -!Config Key = offline -!Config Desc = Nouvelle eau liquide -!Config Def = n -!Config Help = Permet de mettre en route la -!Config nouvelle parametrisation de l'eau liquide ! - offline = .FALSE. - CALL getin('offline',offline) - -!Config Key = type_trac -!Config Desc = Choix de couplage avec model de chimie INCA ou REPROBUS -!Config Def = lmdz -!Config Help = -!Config 'lmdz' = pas de couplage, pur LMDZ -!Config 'inca' = model de chime INCA -!Config 'repr' = model de chime REPROBUS - type_trac = 'lmdz' - CALL getin('type_trac',type_trac) - -!Config Key = config_inca -!Config Desc = Choix de configuration de INCA -!Config Def = none -!Config Help = Choix de configuration de INCA : -!Config 'none' = sans INCA -!Config 'chem' = INCA avec calcul de chemie -!Config 'aero' = INCA avec calcul des aerosols - config_inca = 'none' - CALL getin('config_inca',config_inca) - -!Config Key = ok_dynzon -!Config Desc = calcul et sortie des transports -!Config Def = n -!Config Help = Permet de mettre en route le calcul des transports -!Config - ok_dynzon = .FALSE. - CALL getin('ok_dynzon',ok_dynzon) - -!Config Key = ok_dyn_ins -!Config Desc = sorties instantanees dans la dynamique -!Config Def = n -!Config Help = -!Config - ok_dyn_ins = .FALSE. - CALL getin('ok_dyn_ins',ok_dyn_ins) - -!Config Key = ok_dyn_ave -!Config Desc = sorties moyennes dans la dynamique -!Config Def = n -!Config Help = -!Config - ok_dyn_ave = .FALSE. - CALL getin('ok_dyn_ave',ok_dyn_ave) + + endif ! etatinit +c---------------------------------------- write(lunout,*)' #########################################' - write(lunout,*)' Configuration des parametres du gcm: ' - write(lunout,*)' planet_type = ', planet_type - write(lunout,*)' calend = ', calend - write(lunout,*)' dayref = ', dayref - write(lunout,*)' anneeref = ', anneeref - write(lunout,*)' nday = ', nday - if (less1day) then - write(lunout,*)' Run only for a fraction of day ! ' - write(lunout,*)' fractday = ', fractday - endif - write(lunout,*)' day_step = ', day_step - write(lunout,*)' iperiod = ', iperiod - write(lunout,*)' iconser = ', iconser - write(lunout,*)' iecri = ', iecri - write(lunout,*)' periodav = ', periodav - write(lunout,*)' output_grads_dyn = ', output_grads_dyn - write(lunout,*)' dissip_period = ', dissip_period - write(lunout,*)' lstardis = ', lstardis - write(lunout,*)' nitergdiv = ', nitergdiv - write(lunout,*)' nitergrot = ', nitergrot - write(lunout,*)' niterh = ', niterh - write(lunout,*)' tetagdiv = ', tetagdiv - write(lunout,*)' tetagrot = ', tetagrot - write(lunout,*)' tetatemp = ', tetatemp - write(lunout,*)' coefdis = ', coefdis - write(lunout,*)' purmats = ', purmats - write(lunout,*)' read_start = ', read_start - write(lunout,*)' iflag_phys = ', iflag_phys - write(lunout,*)' iphysiq = ', iphysiq - write(lunout,*)' iflag_trac = ', iflag_trac - write(lunout,*)' clonn = ', clonn - write(lunout,*)' clatt = ', clatt - write(lunout,*)' grossismx = ', grossismx - write(lunout,*)' grossismy = ', grossismy - write(lunout,*)' fxyhypbb = ', fxyhypbb - write(lunout,*)' dzoomxx = ', dzoomxx - write(lunout,*)' dzoomy = ', dzoomyy - write(lunout,*)' tauxx = ', tauxx - write(lunout,*)' tauyy = ', tauyy - write(lunout,*)' offline = ', offline - write(lunout,*)' type_trac = ', type_trac - write(lunout,*)' config_inca = ', config_inca - write(lunout,*)' ok_dynzon = ', ok_dynzon - write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins - write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave - - RETURN -c ............................................... -c -100 CONTINUE -!Config Key = clon -!Config Desc = centre du zoom, longitude -!Config Def = 0 -!Config Help = longitude en degres du centre -!Config du zoom - clon = 0. - CALL getin('clon',clon) - -!Config Key = clat -!Config Desc = centre du zoom, latitude -!Config Def = 0 -!Config Help = latitude en degres du centre du zoom -!Config - clat = 0. - CALL getin('clat',clat) - -!Config Key = grossismx -!Config Desc = zoom en longitude -!Config Def = 1.0 -!Config Help = facteur de grossissement du zoom, -!Config selon la longitude - grossismx = 1.0 - CALL getin('grossismx',grossismx) - -!Config Key = grossismy -!Config Desc = zoom en latitude -!Config Def = 1.0 -!Config Help = facteur de grossissement du zoom, -!Config selon la latitude - grossismy = 1.0 - CALL getin('grossismy',grossismy) - - IF( grossismx.LT.1. ) THEN - write(lunout,*) - & 'conf_gcm: *** ATTENTION !! grossismx < 1 . *** ' - STOP - ELSE - alphax = 1. - 1./ grossismx - ENDIF - - - IF( grossismy.LT.1. ) THEN - write(lunout,*) - & 'conf_gcm: *** ATTENTION !! grossismy < 1 . *** ' - STOP - ELSE - alphay = 1. - 1./ grossismy - ENDIF - - write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay -c -c alphax et alphay sont les anciennes formulat. des grossissements -c -c - -!Config Key = fxyhypb -!Config Desc = Fonction hyperbolique -!Config Def = y -!Config Help = Fonction f(y) hyperbolique si = .true. -!Config sinon sinusoidale - fxyhypb = .TRUE. - CALL getin('fxyhypb',fxyhypb) - -!Config Key = dzoomx -!Config Desc = extension en longitude -!Config Def = 0 -!Config Help = extension en longitude de la zone du zoom -!Config ( fraction de la zone totale) - dzoomx = 0.0 - CALL getin('dzoomx',dzoomx) - -!Config Key = dzoomy -!Config Desc = extension en latitude -!Config Def = 0 -!Config Help = extension en latitude de la zone du zoom -!Config ( fraction de la zone totale) - dzoomy = 0.0 - CALL getin('dzoomy',dzoomy) - -!Config Key = taux -!Config Desc = raideur du zoom en X -!Config Def = 3 -!Config Help = raideur du zoom en X - taux = 3.0 - CALL getin('taux',taux) - -!Config Key = tauy -!Config Desc = raideur du zoom en Y -!Config Def = 3 -!Config Help = raideur du zoom en Y - tauy = 3.0 - CALL getin('tauy',tauy) - -!Config Key = ysinus -!Config IF = !fxyhypb -!Config Desc = Fonction en Sinus -!Config Def = y -!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true. -!Config sinon y = latit. - ysinus = .TRUE. - CALL getin('ysinus',ysinus) -c -!Config Key = offline -!Config Desc = Nouvelle eau liquide -!Config Def = n -!Config Help = Permet de mettre en route la -!Config nouvelle parametrisation de l'eau liquide ! - offline = .FALSE. - CALL getin('offline',offline) - -!Config Key = type_trac -!Config Desc = Choix de couplage avec model de chimie INCA ou REPROBUS -!Config Def = lmdz -!Config Help = -!Config 'lmdz' = pas de couplage, pur LMDZ -!Config 'inca' = model de chime INCA -!Config 'repr' = model de chime REPROBUS - type_trac = 'lmdz' - CALL getin('type_trac',type_trac) - -!Config Key = config_inca -!Config Desc = Choix de configuration de INCA -!Config Def = none -!Config Help = Choix de configuration de INCA : -!Config 'none' = sans INCA -!Config 'chem' = INCA avec calcul de chemie -!Config 'aero' = INCA avec calcul des aerosols - config_inca = 'none' - CALL getin('config_inca',config_inca) - -!Config Key = ok_dynzon -!Config Desc = sortie des transports zonaux dans la dynamique -!Config Def = n -!Config Help = -!Config - ok_dynzon = .FALSE. - CALL getin('ok_dynzon',ok_dynzon) - -!Config Key = ok_dyn_ins -!Config Desc = sorties instantanees dans la dynamique -!Config Def = n -!Config Help = -!Config - ok_dyn_ins = .FALSE. - CALL getin('ok_dyn_ins',ok_dyn_ins) - -!Config Key = ok_dyn_ave -!Config Desc = sorties moyennes dans la dynamique -!Config Def = n -!Config Help = -!Config - ok_dyn_ave = .FALSE. - CALL getin('ok_dyn_ave',ok_dyn_ave) - -!Config Key = use_filtre_fft -!Config Desc = flag d'activation des FFT pour le filtre -!Config Def = false -!Config Help = permet d'activer l'utilisation des FFT pour effectuer -!Config le filtrage aux poles. -! Le filtre fft n'est pas implemente dans dyn3d - use_filtre_fft=.FALSE. - CALL getin('use_filtre_fft',use_filtre_fft) - - IF (use_filtre_fft) THEN - write(lunout,*)'STOP !!!' - write(lunout,*)'use_filtre_fft n est pas implemente dans dyn3d' - STOP - ENDIF - -!Config key = ok_strato -!Config Desc = activation de la version strato -!Config Def = .FALSE. -!Config Help = active la version stratosphérique de LMDZ de F. Lott - - ok_strato=.TRUE. - CALL getin('ok_strato',ok_strato) - -!Config Key = ok_gradsfile -!Config Desc = activation des sorties grads du guidage -!Config Def = n -!Config Help = active les sorties grads du guidage - - ok_gradsfile = .FALSE. - CALL getin('ok_gradsfile',ok_gradsfile) - -!Config Key = ok_limit -!Config Desc = creation des fichiers limit dans create_etat0_limit -!Config Def = y -!Config Help = production du fichier limit.nc requise - - ok_limit = .TRUE. - CALL getin('ok_limit',ok_limit) - -!Config Key = ok_etat0 -!Config Desc = creation des fichiers etat0 dans create_etat0_limit -!Config Def = y -!Config Help = production des fichiers start.nc, startphy.nc requise - - ok_etat0 = .TRUE. - CALL getin('ok_etat0',ok_etat0) - -!Config Key = grilles_gcm_netcdf -!Config Desc = creation de fichier grilles_gcm.nc dans create_etat0_limit -!Config Def = n - grilles_gcm_netcdf = .FALSE. - CALL getin('grilles_gcm_netcdf',grilles_gcm_netcdf) - - write(lunout,*)' #########################################' - write(lunout,*)' Configuration des parametres de cel0' - & //'_limit: ' + write(lunout,*)' Configuration des parametres lus via run.def ' write(lunout,*)' planet_type = ', planet_type write(lunout,*)' calend = ', calend Index: trunk/LMDZ.COMMON/libf/dyn3dpar/conf_gcm.F =================================================================== --- trunk/LMDZ.COMMON/libf/dyn3dpar/conf_gcm.F (revision 500) +++ trunk/LMDZ.COMMON/libf/dyn3dpar/conf_gcm.F (revision 617) @@ -38,5 +38,4 @@ #include "serre.h" #include "comdissnew.h" -!#include "clesphys.h" #include "iniprint.h" #include "temps.h" @@ -44,4 +43,5 @@ ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique +! #include "clesphys.h" c c @@ -201,8 +201,7 @@ !Config Key = nsplit_phys -!Config Desc = nombre d'iteration de la physique -!Config Def = 240 -!Config Help = nombre d'itration de la physique -! +!Config Desc = nombre de subdivisions par pas physique +!Config Def = 1 +!Config Help = nombre de subdivisions par pas physique nsplit_phys = 1 CALL getin('nsplit_phys',nsplit_phys) @@ -422,203 +421,5 @@ ip_ebil_dyn = 0 CALL getin('ip_ebil_dyn',ip_ebil_dyn) -! - - -ccc .... P. Le Van , ajout le 7/03/95 .pour le zoom ... -c ......... ( modif le 17/04/96 ) ......... -c - IF( etatinit ) GO TO 100 - -!Config Key = clon -!Config Desc = centre du zoom, longitude -!Config Def = 0 -!Config Help = longitude en degres du centre -!Config du zoom - clonn = 0. - CALL getin('clon',clonn) - -!Config Key = clat -!Config Desc = centre du zoom, latitude -!Config Def = 0 -!Config Help = latitude en degres du centre du zoom -!Config - clatt = 0. - 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', - & ' est differente de celle lue sur le fichier start ' - STOP - ENDIF - -!Config Key = grossismx -!Config Desc = zoom en longitude -!Config Def = 1.0 -!Config Help = facteur de grossissement du zoom, -!Config selon la longitude - grossismxx = 1.0 - CALL getin('grossismx',grossismxx) - - - IF( ABS(grossismx - grossismxx).GE. 0.001 ) THEN - write(lunout,*)'conf_gcm: La valeur de grossismx passee par ', - & 'run.def est differente de celle lue sur le fichier start ' - STOP - ENDIF - -!Config Key = grossismy -!Config Desc = zoom en latitude -!Config Def = 1.0 -!Config Help = facteur de grossissement du zoom, -!Config selon la latitude - grossismyy = 1.0 - CALL getin('grossismy',grossismyy) - - IF( ABS(grossismy - grossismyy).GE. 0.001 ) THEN - write(lunout,*)'conf_gcm: La valeur de grossismy passee par ', - & 'run.def est differente de celle lue sur le fichier start ' - STOP - ENDIF - - IF( grossismx.LT.1. ) THEN - write(lunout,*) - & 'conf_gcm: *** ATTENTION !! grossismx < 1 . *** ' - STOP - ELSE - alphax = 1. - 1./ grossismx - ENDIF - - - IF( grossismy.LT.1. ) THEN - write(lunout,*) - & 'conf_gcm: *** ATTENTION !! grossismy < 1 . *** ' - STOP - ELSE - alphay = 1. - 1./ grossismy - ENDIF - - write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay -c -c alphax et alphay sont les anciennes formulat. des grossissements -c -c - -!Config Key = fxyhypb -!Config Desc = Fonction hyperbolique -!Config Def = y -!Config Help = Fonction f(y) hyperbolique si = .true. -!Config sinon sinusoidale - fxyhypbb = .TRUE. - CALL getin('fxyhypb',fxyhypbb) - - IF( .NOT.fxyhypb ) THEN - 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 ***' - STOP - ENDIF - ELSE - 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 **** ' - STOP - ENDIF - ENDIF -c -!Config Key = dzoomx -!Config Desc = extension en longitude -!Config Def = 0 -!Config Help = extension en longitude de la zone du zoom -!Config ( fraction de la zone totale) - dzoomxx = 0.0 - CALL getin('dzoomx',dzoomxx) - - 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 ' - STOP - ENDIF - ENDIF - -!Config Key = dzoomy -!Config Desc = extension en latitude -!Config Def = 0 -!Config Help = extension en latitude de la zone du zoom -!Config ( fraction de la zone totale) - dzoomyy = 0.0 - CALL getin('dzoomy',dzoomyy) - - 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 ' - STOP - ENDIF - ENDIF - -!Config Key = taux -!Config Desc = raideur du zoom en X -!Config Def = 3 -!Config Help = raideur du zoom en X - tauxx = 3.0 - CALL getin('taux',tauxx) - - 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 ' - STOP - ENDIF - ENDIF - -!Config Key = tauyy -!Config Desc = raideur du zoom en Y -!Config Def = 3 -!Config Help = raideur du zoom en Y - tauyy = 3.0 - CALL getin('tauy',tauyy) - - 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 ' - STOP - ENDIF - ENDIF - -cc - IF( .NOT.fxyhypb ) THEN - -!Config Key = ysinus -!Config IF = !fxyhypb -!Config Desc = Fonction en Sinus -!Config Def = y -!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true. -!Config sinon y = latit. - ysinuss = .TRUE. - CALL getin('ysinus',ysinuss) - - IF( .NOT.ysinus ) THEN - 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 ***' - STOP - ENDIF - ELSE - 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 **** ' - STOP - ENDIF - ENDIF - ENDIF ! of IF( .NOT.fxyhypb ) -c + !Config Key = offline !Config Desc = Nouvelle eau liquide @@ -681,56 +482,87 @@ CALL getin('ok_dyn_ave',ok_dyn_ave) - write(lunout,*)' #########################################' - write(lunout,*)' Configuration des parametres du gcm: ' - write(lunout,*)' planet_type = ', planet_type - write(lunout,*)' calend = ', calend - write(lunout,*)' dayref = ', dayref - write(lunout,*)' anneeref = ', anneeref - write(lunout,*)' nday = ', nday - if (less1day) then - write(lunout,*)' Run only for a fraction of day ! ' - write(lunout,*)' fractday = ', fractday - endif - write(lunout,*)' day_step = ', day_step - write(lunout,*)' iperiod = ', iperiod - write(lunout,*)' nsplit_phys = ', nsplit_phys - write(lunout,*)' iconser = ', iconser - write(lunout,*)' iecri = ', iecri - write(lunout,*)' periodav = ', periodav - write(lunout,*)' output_grads_dyn = ', output_grads_dyn - write(lunout,*)' dissip_period = ', dissip_period - write(lunout,*)' lstardis = ', lstardis - write(lunout,*)' nitergdiv = ', nitergdiv - write(lunout,*)' nitergrot = ', nitergrot - write(lunout,*)' niterh = ', niterh - write(lunout,*)' tetagdiv = ', tetagdiv - write(lunout,*)' tetagrot = ', tetagrot - write(lunout,*)' tetatemp = ', tetatemp - write(lunout,*)' coefdis = ', coefdis - write(lunout,*)' purmats = ', purmats - write(lunout,*)' read_start = ', read_start - write(lunout,*)' iflag_phys = ', iflag_phys - write(lunout,*)' iphysiq = ', iphysiq - write(lunout,*)' iflag_trac = ', iflag_trac - write(lunout,*)' clonn = ', clonn - write(lunout,*)' clatt = ', clatt - write(lunout,*)' grossismx = ', grossismx - write(lunout,*)' grossismy = ', grossismy - write(lunout,*)' fxyhypbb = ', fxyhypbb - write(lunout,*)' dzoomxx = ', dzoomxx - write(lunout,*)' dzoomy = ', dzoomyy - write(lunout,*)' tauxx = ', tauxx - write(lunout,*)' tauyy = ', tauyy - write(lunout,*)' offline = ', offline - write(lunout,*)' type_trac = ', type_trac - write(lunout,*)' config_inca = ', config_inca - write(lunout,*)' ok_dynzon = ', ok_dynzon - write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins - write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave - - RETURN -c ............................................... -c -100 CONTINUE +!Config Key = use_filtre_fft +!Config Desc = flag d'activation des FFT pour le filtre +!Config Def = false +!Config Help = permet d'activer l'utilisation des FFT pour effectuer +!Config le filtrage aux poles. + use_filtre_fft=.FALSE. + CALL getin('use_filtre_fft',use_filtre_fft) + + IF (use_filtre_fft .AND. grossismx /= 1.0) THEN + write(lunout,*)'WARNING !!! ' + write(lunout,*)"Le zoom en longitude est incompatible", + & " avec l'utilisation du filtre FFT ", + & "---> filtre FFT désactivé " + use_filtre_fft=.FALSE. + ENDIF + +!Config Key = use_mpi_alloc +!Config Desc = Utilise un buffer MPI en memoire globale +!Config Def = false +!Config Help = permet d'activer l'utilisation d'un buffer MPI +!Config en memoire globale a l'aide de la fonction MPI_ALLOC. +!Config Cela peut ameliorer la bande passante des transferts MPI +!Config d'un facteur 2 + use_mpi_alloc=.FALSE. + CALL getin('use_mpi_alloc',use_mpi_alloc) + +!Config Key = omp_chunk +!Config Desc = taille des blocs openmp +!Config Def = 1 +!Config Help = defini la taille des packets d'iteration openmp +!Config distribuee a chaque tache lors de l'entree dans une +!Config boucle parallelisee + + omp_chunk=1 + CALL getin('omp_chunk',omp_chunk) + +!Config key = ok_strato +!Config Desc = activation de la version strato +!Config Def = .FALSE. +!Config Help = active la version stratosphérique de LMDZ de F. Lott + + ok_strato=.TRUE. + CALL getin('ok_strato',ok_strato) + +!Config Key = ok_gradsfile +!Config Desc = activation des sorties grads du guidage +!Config Def = n +!Config Help = active les sorties grads du guidage + + ok_gradsfile = .FALSE. + CALL getin('ok_gradsfile',ok_gradsfile) + +!Config Key = ok_limit +!Config Desc = creation des fichiers limit dans create_etat0_limit +!Config Def = y +!Config Help = production du fichier limit.nc requise + + ok_limit = .TRUE. + CALL getin('ok_limit',ok_limit) + +!Config Key = ok_etat0 +!Config Desc = creation des fichiers etat0 dans create_etat0_limit +!Config Def = y +!Config Help = production des fichiers start.nc, startphy.nc requise + + ok_etat0 = .TRUE. + CALL getin('ok_etat0',ok_etat0) + +!Config Key = grilles_gcm_netcdf +!Config Desc = creation de fichier grilles_gcm.nc dans create_etat0_limit +!Config Def = n + grilles_gcm_netcdf = .FALSE. + CALL getin('grilles_gcm_netcdf',grilles_gcm_netcdf) + +c---------------------------------------- +ccc .... P. Le Van , ajout le 7/03/95 .pour le zoom ... +c ......... ( modif le 17/04/96 ) ......... +c +C ZOOM PARAMETERS ... the ones read in start.nc prevail anyway ! (SL, 2012) +c +c---------------------------------------- + IF( etatinit ) then + !Config Key = clon !Config Desc = centre du zoom, longitude @@ -835,143 +667,205 @@ CALL getin('ysinus',ysinus) c -!Config Key = offline -!Config Desc = Nouvelle eau liquide -!Config Def = n -!Config Help = Permet de mettre en route la -!Config nouvelle parametrisation de l'eau liquide ! - offline = .FALSE. - CALL getin('offline',offline) - IF (offline .AND. adjust) THEN - WRITE(lunout,*) - & 'WARNING : option offline does not work with adjust=y :' - WRITE(lunout,*) 'the files defstoke.nc, fluxstoke.nc ', - & 'and fluxstokev.nc will not be created' - WRITE(lunout,*) - & 'only the file phystoke.nc will still be created ' - END IF - -!Config Key = type_trac -!Config Desc = Choix de couplage avec model de chimie INCA ou REPROBUS -!Config Def = lmdz -!Config Help = -!Config 'lmdz' = pas de couplage, pur LMDZ -!Config 'inca' = model de chime INCA -!Config 'repr' = model de chime REPROBUS - type_trac = 'lmdz' - CALL getin('type_trac',type_trac) - -!Config Key = config_inca -!Config Desc = Choix de configuration de INCA -!Config Def = none -!Config Help = Choix de configuration de INCA : -!Config 'none' = sans INCA -!Config 'chem' = INCA avec calcul de chemie -!Config 'aero' = INCA avec calcul des aerosols - config_inca = 'none' - CALL getin('config_inca',config_inca) - -!Config Key = ok_dynzon -!Config Desc = calcul et sortie des transports -!Config Def = n -!Config Help = Permet de mettre en route le calcul des transports -!Config - ok_dynzon = .FALSE. - CALL getin('ok_dynzon',ok_dynzon) - -!Config Key = ok_dyn_ins -!Config Desc = sorties instantanees dans la dynamique -!Config Def = n -!Config Help = -!Config - ok_dyn_ins = .FALSE. - CALL getin('ok_dyn_ins',ok_dyn_ins) - -!Config Key = ok_dyn_ave -!Config Desc = sorties moyennes dans la dynamique -!Config Def = n -!Config Help = -!Config - ok_dyn_ave = .FALSE. - CALL getin('ok_dyn_ave',ok_dyn_ave) - -!Config Key = use_filtre_fft -!Config Desc = flag d'activation des FFT pour le filtre -!Config Def = false -!Config Help = permet d'activer l'utilisation des FFT pour effectuer -!Config le filtrage aux poles. - use_filtre_fft=.FALSE. - CALL getin('use_filtre_fft',use_filtre_fft) - - IF (use_filtre_fft .AND. grossismx /= 1.0) THEN - write(lunout,*)'WARNING !!! ' - write(lunout,*)"Le zoom en longitude est incompatible", - & " avec l'utilisation du filtre FFT ", - & "---> filtre FFT désactivé " - use_filtre_fft=.FALSE. +c---------------------------------------- + else ! etatinit=false +c---------------------------------------- + +!Config Key = clon +!Config Desc = centre du zoom, longitude +!Config Def = 0 +!Config Help = longitude en degres du centre +!Config du zoom + clonn = 0. + CALL getin('clon',clonn) + +!Config Key = clat +!Config Desc = centre du zoom, latitude +!Config Def = 0 +!Config Help = latitude en degres du centre du zoom +!Config + clatt = 0. + 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', + & ' est differente de celle lue sur le fichier start ' + STOP + ENDIF + +!Config Key = grossismx +!Config Desc = zoom en longitude +!Config Def = 1.0 +!Config Help = facteur de grossissement du zoom, +!Config selon la longitude + grossismxx = 1.0 + CALL getin('grossismx',grossismxx) + + + IF( ABS(grossismx - grossismxx).GE. 0.001 ) THEN + write(lunout,*)'conf_gcm: La valeur de grossismx passee par ', + & 'run.def est differente de celle lue sur le fichier start ' + STOP + ENDIF + +!Config Key = grossismy +!Config Desc = zoom en latitude +!Config Def = 1.0 +!Config Help = facteur de grossissement du zoom, +!Config selon la latitude + grossismyy = 1.0 + CALL getin('grossismy',grossismyy) + + IF( ABS(grossismy - grossismyy).GE. 0.001 ) THEN + write(lunout,*)'conf_gcm: La valeur de grossismy passee par ', + & 'run.def est differente de celle lue sur le fichier start ' + STOP ENDIF - + IF( grossismx.LT.1. ) THEN + write(lunout,*) + & 'conf_gcm: *** ATTENTION !! grossismx < 1 . *** ' + STOP + ELSE + alphax = 1. - 1./ grossismx + ENDIF + + + IF( grossismy.LT.1. ) THEN + write(lunout,*) + & 'conf_gcm: *** ATTENTION !! grossismy < 1 . *** ' + STOP + ELSE + alphay = 1. - 1./ grossismy + ENDIF + + write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay +c +c alphax et alphay sont les anciennes formulat. des grossissements +c +c + +!Config Key = fxyhypb +!Config Desc = Fonction hyperbolique +!Config Def = y +!Config Help = Fonction f(y) hyperbolique si = .true. +!Config sinon sinusoidale + fxyhypbb = .TRUE. + CALL getin('fxyhypb',fxyhypbb) + + IF( .NOT.fxyhypb ) THEN + 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 ***' + STOP + ENDIF + ELSE + 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 **** ' + STOP + ENDIF + ENDIF +c +!Config Key = dzoomx +!Config Desc = extension en longitude +!Config Def = 0 +!Config Help = extension en longitude de la zone du zoom +!Config ( fraction de la zone totale) + dzoomxx = 0.0 + CALL getin('dzoomx',dzoomxx) + + 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 ' + STOP + ENDIF + ENDIF + +!Config Key = dzoomy +!Config Desc = extension en latitude +!Config Def = 0 +!Config Help = extension en latitude de la zone du zoom +!Config ( fraction de la zone totale) + dzoomyy = 0.0 + CALL getin('dzoomy',dzoomyy) + + 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 ' + STOP + ENDIF + ENDIF -!Config Key = use_mpi_alloc -!Config Desc = Utilise un buffer MPI en m�moire globale -!Config Def = false -!Config Help = permet d'activer l'utilisation d'un buffer MPI -!Config en m�moire globale a l'aide de la fonction MPI_ALLOC. -!Config Cela peut am�liorer la bande passante des transferts MPI -!Config d'un facteur 2 - use_mpi_alloc=.FALSE. - CALL getin('use_mpi_alloc',use_mpi_alloc) - -!Config Key = omp_chunk -!Config Desc = taille des blocs openmp -!Config Def = 1 -!Config Help = defini la taille des packets d'it�ration openmp -!Config distribu�e � chaque t�che lors de l'entr�e dans une -!Config boucle parall�lis�e - - omp_chunk=1 - CALL getin('omp_chunk',omp_chunk) - -!Config key = ok_strato -!Config Desc = activation de la version strato -!Config Def = .FALSE. -!Config Help = active la version stratosphérique de LMDZ de F. Lott - - ok_strato=.TRUE. - CALL getin('ok_strato',ok_strato) - -!Config Key = ok_gradsfile -!Config Desc = activation des sorties grads du guidage -!Config Def = n -!Config Help = active les sorties grads du guidage - - ok_gradsfile = .FALSE. - CALL getin('ok_gradsfile',ok_gradsfile) - -!Config Key = ok_limit -!Config Desc = creation des fichiers limit dans create_etat0_limit +!Config Key = taux +!Config Desc = raideur du zoom en X +!Config Def = 3 +!Config Help = raideur du zoom en X + tauxx = 3.0 + CALL getin('taux',tauxx) + + 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 ' + STOP + ENDIF + ENDIF + +!Config Key = tauyy +!Config Desc = raideur du zoom en Y +!Config Def = 3 +!Config Help = raideur du zoom en Y + tauyy = 3.0 + CALL getin('tauy',tauyy) + + 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 ' + STOP + ENDIF + ENDIF + +cc + IF( .NOT.fxyhypb ) THEN + +!Config Key = ysinus +!Config IF = !fxyhypb +!Config Desc = Fonction en Sinus !Config Def = y -!Config Help = production du fichier limit.nc requise - - ok_limit = .TRUE. - CALL getin('ok_limit',ok_limit) - -!Config Key = ok_etat0 -!Config Desc = creation des fichiers etat0 dans create_etat0_limit -!Config Def = y -!Config Help = production des fichiers start.nc, startphy.nc requise - - ok_etat0 = .TRUE. - CALL getin('ok_etat0',ok_etat0) - -!Config Key = grilles_gcm_netcdf -!Config Desc = creation de fichier grilles_gcm.nc dans create_etat0_limit -!Config Def = n - grilles_gcm_netcdf = .FALSE. - CALL getin('grilles_gcm_netcdf',grilles_gcm_netcdf) +!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true. +!Config sinon y = latit. + ysinuss = .TRUE. + CALL getin('ysinus',ysinuss) + + IF( .NOT.ysinus ) THEN + 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 ***' + STOP + ENDIF + ELSE + 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 **** ' + STOP + ENDIF + ENDIF + ENDIF ! of IF( .NOT.fxyhypb ) + + endif ! etatinit +c---------------------------------------- + write(lunout,*)' #########################################' - write(lunout,*)' Configuration des parametres de cel0' - & //'_limit: ' + write(lunout,*)' Configuration des parametres lus via run.def ' write(lunout,*)' planet_type = ', planet_type write(lunout,*)' calend = ', calend