trunk/dyn3d/conf_gcm.f

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:: iperiod= 5
  ! periode pour le pas Matsuno (en pas de temps)

  integer:: iapp_tracvl= 5
  ! Should normally be equal to "iperiod"
  ! frequence du groupement des flux (en pas de temps) 

  integer:: iconser= 240
  ! periode de sortie des variables de controle
  ! (En pas de temps)

  integer:: iecri= 1 ! période d'écriture du fichier "dyn_hist.nc" (en jours)

  integer:: idissip= 10 ! periode de la dissipation (en pas)

  integer:: iphysiq= 5
  ! Help = Periode de la physique en pas de temps de la dynamique.

  integer:: dayref = 1 ! jour de l'année de l'état initial
  ! (= 350 si 20 décembre par exemple)

  integer:: anneeref = 1998 ! Annee de l'etat initial (avec 4 chiffres)

  integer:: raz_date = 0 ! Remise a zero de la date initiale
  ! 0 pas de remise a zero, on garde la date du fichier restart
  ! 1 prise en compte de la date de gcm.def avec remise a zero
  ! des compteurs de pas de temps

  REAL:: periodav= 1.
  ! periode de stockage fichier histmoy (en jour) 

  logical, save:: offline

contains

  SUBROUTINE conf_gcm(clesphy0)

    ! Auteurs : L. Fairhead, P. Le Van
    ! Version du 29/04/97

    ! Nouveaux paramètres nitergdiv, nitergrot, niterh, tetagdiv, tetagrot, 
    ! tetatemp ajoutés pour la dissipation.

    ! Autre paramètre ajouté en fin de liste de tapedef : fxyhypb
    ! Si fxyhypb = .TRUE., choix de la fonction à dérivée tangente
    ! hyperbolique
    ! Sinon, choix de fxynew, à dérivée sinusoïdale

    ! On ne compare pas les valeurs des paramètres du zoom, grossismx, 
    ! grossismy, clon, clat, fxyhypb lues sur le fichier start avec
    ! celles passées par run.def, au début du gcm. 
    ! Ces paramètres définissent entre autres la grille et doivent être
    ! cohérents, sinon il y aura divergence du gcm.

    use IOIPSL, only: getin
    use dimens_m
    use paramet_m
    use comdissnew, only: read_comdissnew
    use logic, only: iflag_phys, fxyhypb, ysinus, purmats, ok_guide
    use serre, only: clon, clat, grossismx, grossismy, alphax, alphay, &
         dzoomx, dzoomy, taux, tauy
    use clesphys, only: ok_limitvrai, ok_orolf, iflag_con, nbapp_rad, &
         ok_orodr, cycle_diurne, new_oliq, soil_model
    use iniprint, only: lunout, read_iniprint

    INTEGER, PARAMETER:: longcles = 20

    REAL, intent(out), optional:: clesphy0(longcles)

    namelist /conf_gcm_nml/dayref, anneeref, raz_date, nday, day_step, &
         iperiod, iapp_tracvl, iconser, iecri, periodav, idissip, purmats, &
         ok_guide, iflag_phys, iphysiq, cycle_diurne, soil_model, new_oliq, &
         ok_orodr, ok_orolf, ok_limitvrai, nbapp_rad, iflag_con, clon, clat, &
         grossismx, grossismy, dzoomx, dzoomy, taux, tauy

    !------------------------------------

    print *, "Call sequence information: conf_gcm"
    call read_iniprint

    print *, "Enter namelist 'conf_gcm_nml'."
    read(unit=*, nml=conf_gcm_nml)
    write(unit=*, nml=conf_gcm_nml)

    call read_comdissnew

    IF (lunout /= 5 .and. lunout /= 6) THEN
       OPEN(lunout, FILE='lmdz.out')
    ENDIF

    if (present(clesphy0)) then
       clesphy0(:) = 0.
       clesphy0(1) = REAL(iflag_con)
       clesphy0(2) = REAL(nbapp_rad)

       IF (cycle_diurne) clesphy0(3) = 1.
       IF (soil_model) clesphy0(4) = 1.
       IF (new_oliq) clesphy0(5) = 1.
       IF (ok_orodr) clesphy0(6) = 1.
       IF (ok_orolf) clesphy0(7) = 1.
       IF (ok_limitvrai) clesphy0(8) = 1.
    end if

    IF (grossismx < 1.) THEN
       PRINT *, 'Error: grossismx < 1'
       STOP 1
    ELSE
       alphax = 1. - 1. / grossismx
    ENDIF
    IF (grossismy < 1.) THEN
       PRINT *, 'Error: grossismy < 1'
       STOP 1
    ELSE
       alphay = 1. - 1. / grossismy
    ENDIF
    PRINT *, 'alphax = ', alphax
    PRINT *, 'alphay = ', alphay

    ! Key = fxyhypb
    ! Desc = Fonction hyperbolique
    ! Def = y
    ! Help = Fonction f(y) hyperbolique si = .true. 
    ! sinon sinusoidale
    fxyhypb = .TRUE.
    CALL getin('fxyhypb', fxyhypb)

    ! Key = ysinus
    ! IF = !fxyhypb
    ! Desc = Fonction en Sinus
    ! Def = y
    ! Help = Fonction f(y) avec y = Sin(latit.) si = .true. 
    ! sinon y = latit.
    ysinus = .TRUE.
    CALL getin('ysinus', ysinus)

    ! Key = offline
    ! Desc = Nouvelle eau liquide
    ! Def = n
    ! Help = Permet de mettre en route la
    ! nouvelle parametrisation de l'eau liquide
    offline = .FALSE.
    CALL getin('offline', offline)
    write(lunout, *)' offline = ', offline

  END SUBROUTINE conf_gcm

end module conf_gcm_m
1	guez	3	module conf_gcm_m
2
3			! This module is clean: no C preprocessor directive, no include line
4
5			IMPLICIT NONE
6
7			INTEGER:: nday= 10
8			! Nombre de jours d'integration
9			! On pourait aussi permettre des mois ou des annees !
10
11			integer:: day_step= 240 ! nombre de pas par jour, multiple de iperiod
12
13			integer:: iperiod= 5
14			! periode pour le pas Matsuno (en pas de temps)
15
16			integer:: iapp_tracvl= 5
17			! Should normally be equal to "iperiod"
18			! frequence du groupement des flux (en pas de temps)
19
20			integer:: iconser= 240
21			! periode de sortie des variables de controle
22			! (En pas de temps)
23
24			integer:: iecri= 1 ! période d'écriture du fichier "dyn_hist.nc" (en jours)
25
26			integer:: idissip= 10 ! periode de la dissipation (en pas)
27
28			integer:: iphysiq= 5
29			! Help = Periode de la physique en pas de temps de la dynamique.
30
31			integer:: dayref = 1 ! jour de l'année de l'état initial
32			! (= 350 si 20 décembre par exemple)
33
34			integer:: anneeref = 1998 ! Annee de l'etat initial (avec 4 chiffres)
35
36			integer:: raz_date = 0 ! Remise a zero de la date initiale
37			! 0 pas de remise a zero, on garde la date du fichier restart
38			! 1 prise en compte de la date de gcm.def avec remise a zero
39			! des compteurs de pas de temps
40
41			REAL:: periodav= 1.
42			! periode de stockage fichier histmoy (en jour)
43
44			logical, save:: offline
45
46			contains
47
48			SUBROUTINE conf_gcm(clesphy0)
49
50			! Auteurs : L. Fairhead, P. Le Van
51			! Version du 29/04/97
52
53			! Nouveaux paramètres nitergdiv, nitergrot, niterh, tetagdiv, tetagrot,
54			! tetatemp ajoutés pour la dissipation.
55
56			! Autre paramètre ajouté en fin de liste de tapedef : fxyhypb
57			! Si fxyhypb = .TRUE., choix de la fonction à dérivée tangente
58			! hyperbolique
59			! Sinon, choix de fxynew, à dérivée sinusoïdale
60
61			! On ne compare pas les valeurs des paramètres du zoom, grossismx,
62			! grossismy, clon, clat, fxyhypb lues sur le fichier start avec
63			! celles passées par run.def, au début du gcm.
64			! Ces paramètres définissent entre autres la grille et doivent être
65			! cohérents, sinon il y aura divergence du gcm.
66
67			use IOIPSL, only: getin
68			use dimens_m
69			use paramet_m
70			use comdissnew, only: read_comdissnew
71			use logic, only: iflag_phys, fxyhypb, ysinus, purmats, ok_guide
72			use serre, only: clon, clat, grossismx, grossismy, alphax, alphay, &
73			dzoomx, dzoomy, taux, tauy
74			use clesphys, only: ok_limitvrai, ok_orolf, iflag_con, nbapp_rad, &
75			ok_orodr, cycle_diurne, new_oliq, soil_model
76			use iniprint, only: lunout, read_iniprint
77
78			INTEGER, PARAMETER:: longcles = 20
79
80			REAL, intent(out), optional:: clesphy0(longcles)
81
82			namelist /conf_gcm_nml/dayref, anneeref, raz_date, nday, day_step, &
83			iperiod, iapp_tracvl, iconser, iecri, periodav, idissip, purmats, &
84			ok_guide, iflag_phys, iphysiq, cycle_diurne, soil_model, new_oliq, &
85			ok_orodr, ok_orolf, ok_limitvrai, nbapp_rad, iflag_con, clon, clat, &
86			grossismx, grossismy, dzoomx, dzoomy, taux, tauy
87
88			!------------------------------------
89
90			print *, "Call sequence information: conf_gcm"
91			call read_iniprint
92
93			print *, "Enter namelist 'conf_gcm_nml'."
94			read(unit=*, nml=conf_gcm_nml)
95			write(unit=*, nml=conf_gcm_nml)
96
97			call read_comdissnew
98
99			IF (lunout /= 5 .and. lunout /= 6) THEN
100			OPEN(lunout, FILE='lmdz.out')
101			ENDIF
102
103			if (present(clesphy0)) then
104			clesphy0(:) = 0.
105			clesphy0(1) = REAL(iflag_con)
106			clesphy0(2) = REAL(nbapp_rad)
107
108			IF (cycle_diurne) clesphy0(3) = 1.
109			IF (soil_model) clesphy0(4) = 1.
110			IF (new_oliq) clesphy0(5) = 1.
111			IF (ok_orodr) clesphy0(6) = 1.
112			IF (ok_orolf) clesphy0(7) = 1.
113			IF (ok_limitvrai) clesphy0(8) = 1.
114			end if
115
116			IF (grossismx < 1.) THEN
117			PRINT *, 'Error: grossismx < 1'
118			STOP 1
119			ELSE
120			alphax = 1. - 1. / grossismx
121			ENDIF
122			IF (grossismy < 1.) THEN
123			PRINT *, 'Error: grossismy < 1'
124			STOP 1
125			ELSE
126			alphay = 1. - 1. / grossismy
127			ENDIF
128			PRINT *, 'alphax = ', alphax
129			PRINT *, 'alphay = ', alphay
130
131			! Key = fxyhypb
132			! Desc = Fonction hyperbolique
133			! Def = y
134			! Help = Fonction f(y) hyperbolique si = .true.
135			! sinon sinusoidale
136			fxyhypb = .TRUE.
137			CALL getin('fxyhypb', fxyhypb)
138
139			! Key = ysinus
140			! IF = !fxyhypb
141			! Desc = Fonction en Sinus
142			! Def = y
143			! Help = Fonction f(y) avec y = Sin(latit.) si = .true.
144			! sinon y = latit.
145			ysinus = .TRUE.
146			CALL getin('ysinus', ysinus)
147
148			! Key = offline
149			! Desc = Nouvelle eau liquide
150			! Def = n
151			! Help = Permet de mettre en route la
152			! nouvelle parametrisation de l'eau liquide
153			offline = .FALSE.
154			CALL getin('offline', offline)
155			write(lunout, *)' offline = ', offline
156
157			END SUBROUTINE conf_gcm
158
159			end module conf_gcm_m