--- trunk/libf/dyn3d/conf_gcm.f90 2010/03/25 14:29:07 27 +++ trunk/libf/dyn3d/conf_gcm.f90 2011/02/22 13:49:36 40 @@ -1,29 +1,24 @@ module conf_gcm_m - ! This module is clean: no C preprocessor directive, no include line - IMPLICIT NONE - INTEGER:: nday= 10 - ! Nombre de jours d'integration - ! On pourait aussi permettre des mois ou des annees ! - - integer:: day_step= 240 ! nombre de pas par jour, multiple de iperiod + INTEGER:: nday = 10 ! nombre de jours d'intégration + integer:: day_step = 240 ! nombre de pas par jour, multiple de iperiod - integer:: iperiod= 5 + integer:: iperiod = 5 ! periode pour le pas Matsuno (en pas de temps) - integer:: iapp_tracvl= 5 + integer:: iapp_tracvl = 5 ! Should normally be equal to "iperiod" ! frequence du groupement des flux (en pas de temps) - integer:: iconser= 240 - ! période de sortie des variables de contrôle (en pas de temps) + integer:: iconser = 240 + ! number of time steps between output of control variables - integer:: iecri= 1 ! période d'écriture du fichier "dyn_hist.nc" (en jours) + integer:: iecri = 1 ! time interval between outputs to "dyn_hist.nc" (in days) - integer:: iphysiq= 5 - ! Période de la physique en pas de temps de la dynamique. + integer:: iphysiq = 5 + ! number of time steps of dynamics between calls to physics integer:: dayref = 1 ! jour de l'année de l'état initial ! (= 350 si 20 décembre par exemple) @@ -35,12 +30,10 @@ ! gcm.def avec remise a zero des compteurs de pas de temps) ! (pas de remise a zero: on garde la date du fichier restart) - REAL:: periodav= 1. - ! periode de stockage fichier histmoy (en jour) + REAL:: periodav = 1. ! time interval between outputs to "histmoy" (in days) logical:: offline = .FALSE. - ! Nouvelle eau liquide - ! Permet de mettre en route la nouvelle parametrisation de l'eau liquide + ! permet de mettre en route la nouvelle parametrisation de l'eau liquide contains @@ -55,6 +48,7 @@ ! Ces paramètres définissent entre autres la grille et doivent être ! cohérents, sinon il y aura divergence du gcm. + use abort_gcm_m, only: abort_gcm use comdissnew, only: read_comdissnew use logic, only: read_logic use serre, only: clon, clat, grossismx, grossismy, alphax, alphay, & @@ -78,6 +72,16 @@ read(unit=*, nml=conf_gcm_nml) write(unit=*, nml=conf_gcm_nml) + IF (MOD(day_step, iperiod) /= 0) call abort_gcm(modname = "conf_gcm", & + message = & + 'Il faut choisir un nombre de pas par jour multiple de "iperiod".', & + ierr = 1) + + IF (MOD(day_step,iphysiq)/= 0) call abort_gcm(modname = "conf_gcm", & + message = & + 'Il faut choisir un nombre de pas par jour multiple de "iphysiq".', & + ierr = 1) + IF (grossismx < 1.) THEN PRINT *, 'Error: grossismx < 1' STOP 1