/[lmdze]/trunk/Sources/dyn3d/Guide/conf_guide.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/Guide/conf_guide.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/libf/dyn3d/conf_guide.f90 revision 44 by guez, Wed Apr 13 12:29:18 2011 UTC trunk/Sources/dyn3d/Guide/conf_guide.f revision 210 by guez, Tue Dec 13 16:02:23 2016 UTC
# Line 2  module conf_guide_m Line 2  module conf_guide_m
2    
3    IMPLICIT NONE    IMPLICIT NONE
4    
5    REAL tau_min_u, tau_max_u    ! Constantes de rappel, en jours :
6    REAL tau_min_v, tau_max_v    REAL:: tau_min_u = 0.03
7    REAL tau_min_t, tau_max_t    REAL:: tau_max_u = 10.
8    REAL tau_min_q, tau_max_q    REAL:: tau_min_v = 0.03
9    REAL tau_min_p, tau_max_p    REAL:: tau_max_v = 10.
10    LOGICAL ncep, ini_anal    REAL:: tau_min_t = 0.03
11    LOGICAL guide_u, guide_v, guide_t, guide_q, guide_p    REAL:: tau_max_t = 10.
12    INTEGER online    REAL:: tau_min_q = 0.03
13      REAL:: tau_max_q = 10.
14    
15      LOGICAL:: ini_anal = .false. ! Initial = analyse
16      LOGICAL:: guide_u = .false. ! guidage de u
17      LOGICAL:: guide_v = .false. ! gvidage de v
18      LOGICAL:: guide_t = .false. ! guidage de T
19      LOGICAL:: guide_q = .false. ! guidage de q
20    
21      logical:: online = .true. ! contr\^ole du guidage
22      ! hors-ligne: x = x_rea
23    
24      ! Dans le cas où on n'a les analyses que sur une bande de latitudes :
25      REAL, save:: lat_min_guide ! minimum latitude for nudging, in rad
26      real, save:: lat_max_guide ! maximum latitude for nudging, in rad
27    
28      logical, save:: ok_guide ! guidage
29      REAL, save:: factt ! pas de temps entre deux appels au guidage, en jours
30    
31  contains  contains
32    
# Line 18  contains Line 35  contains
35      ! From LMDZ4/libf/dyn3d/conf_guide.F, version 1.1.1.1 2004/05/19 12:53:07      ! From LMDZ4/libf/dyn3d/conf_guide.F, version 1.1.1.1 2004/05/19 12:53:07
36      !  Parametres de controle du run:      !  Parametres de controle du run:
37    
38      use getparam, only: ini_getparam, getpar, fin_getparam      use abort_gcm_m, only: abort_gcm
39      use tau2alpha_m, only: lat_max_guide, lat_min_guide      use comconst, only: daysec, dtvr
40        use conf_gcm_m, only: day_step, iperiod
41        use nr_util, only: assert, pi
42        use unit_nml_m, only: unit_nml
43    
44        ! Local:
45    
46        REAL:: lat_min_guide_deg = -90. ! in degrees
47        real:: lat_max_guide_deg = 90. ! in degrees
48    
49        namelist /conf_guide_nml/ ini_anal, guide_u, guide_v, guide_t, guide_q, &
50             online, tau_min_u, tau_max_u, tau_min_v, tau_max_v, tau_min_t, &
51             tau_max_t, tau_min_q, tau_max_q, lat_min_guide_deg, lat_max_guide_deg
52    
53      !-----------------------------------------------------------------------      !-----------------------------------------------------------------------
54    
55      print *, "Call sequence information: conf_guide"      print *, "Call sequence information: conf_guide"
     call ini_getparam('guide.eff')  
   
     call getpar('online',1,online,'Index de controle du guide')  
     CALL getpar('ncep',.false.,ncep,'Coordonnee vert NCEP ou ECMWF')  
     CALL getpar('ini_anal',.false.,ini_anal,'Initial = analyse')  
   
     CALL getpar('guide_u',.true.,guide_u,'guidage de u')  
     CALL getpar('guide_v',.true.,guide_v,'guidage de v')  
     CALL getpar('guide_T',.true.,guide_T,'guidage de T')  
     CALL getpar('guide_P',.true.,guide_P,'guidage de P')  
     CALL getpar('guide_Q',.true.,guide_Q,'guidage de Q')  
   
     !   Constantes de rappel. Unite : fraction de jour  
     CALL getpar('tau_min_u',0.02,tau_min_u,'Cste de rappel min, u')  
     CALL getpar('tau_max_u', 10.,tau_max_u,'Cste de rappel max, u')  
     CALL getpar('tau_min_v',0.02,tau_min_v,'Cste de rappel min, v')  
     CALL getpar('tau_max_v', 10.,tau_max_v,'Cste de rappel max, v')  
     CALL getpar('tau_min_T',0.02,tau_min_T,'Cste de rappel min, T')  
     CALL getpar('tau_max_T', 10.,tau_max_T,'Cste de rappel max, T')  
     CALL getpar('tau_min_Q',0.02,tau_min_Q,'Cste de rappel min, Q')  
     CALL getpar('tau_max_Q', 10.,tau_max_Q,'Cste de rappel max, Q')  
     CALL getpar('tau_min_P',0.02,tau_min_P,'Cste de rappel min, P')  
     CALL getpar('tau_max_P', 10.,tau_max_P,'Cste de rappel max, P')  
   
     !   Latitude min et max pour le rappel.  
     !   dans le cas ou on 'a les analyses que sur une bande de latitudes.  
     CALL getpar('lat_min_guide',-90.,lat_min_guide &  
          ,'Latitude minimum pour le guidage ')  
     CALL getpar('lat_max_guide', 90.,lat_max_guide &  
          ,'Latitude maximum pour le guidage ')  
56    
57      CALL fin_getparam      print *, "Enter namelist 'conf_guide_nml'."
58        read(unit=*, nml=conf_guide_nml)
59        write(unit_nml, nml=conf_guide_nml)
60    
61        lat_min_guide = lat_min_guide_deg / 180. * pi
62        lat_max_guide = lat_max_guide_deg / 180. * pi
63    
64        ok_guide = any((/guide_u, guide_v, guide_t, guide_q/))
65        if (ok_guide .and. mod(day_step, 4 * iperiod) /= 0) call &
66             abort_gcm("conf_guide", 'ok_guide day_step iperiod')
67    
68        if (ok_guide .and. online) then
69           factt = dtvr * iperiod / daysec
70           print *, "factt = ", factt
71        end if
72    
73    end SUBROUTINE conf_guide    end SUBROUTINE conf_guide
74    
Legend:
Removed from v.44  
changed lines
  Added in v.210
  ViewVC Help
Powered by ViewVC 1.1.21