/[lmdze]/trunk/dyn3d/Guide/conf_guide.f90
ViewVC logotype

Diff of /trunk/dyn3d/Guide/conf_guide.f90

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

trunk/dyn3d/conf_guide.f90 revision 76 by guez, Fri Nov 15 18:45:49 2013 UTC trunk/Sources/dyn3d/Guide/conf_guide.f revision 221 by guez, Thu Apr 20 14:44:47 2017 UTC
# Line 1  Line 1 
1  module conf_guide_m  module conf_guide_m
2    
3      USE dimens_m, ONLY: iim, jjm
4    
5    IMPLICIT NONE    IMPLICIT NONE
6    
7    REAL tau_min_u, tau_max_u    LOGICAL:: ini_anal = .false. ! Initial = analyse
8    REAL tau_min_v, tau_max_v    LOGICAL:: guide_u = .false. ! guidage de u
9    REAL tau_min_t, tau_max_t    LOGICAL:: guide_v = .false. ! gvidage de v
10    REAL tau_min_q, tau_max_q    LOGICAL:: guide_t = .false. ! guidage de T
11    REAL tau_min_p, tau_max_p    LOGICAL:: guide_q = .false. ! guidage de q
12    LOGICAL ncep, ini_anal  
13    LOGICAL guide_u, guide_v, guide_t, guide_q, guide_p    logical, save:: ok_guide ! guidage
14    INTEGER online  
15      ! alpha détermine la part des injections de données à chaque étape
16      ! alpha=0 signifie pas d'injection
17      ! alpha=1 signifie injection totale
18      REAL, save:: alpha_q(iim + 1, jjm + 1)
19      REAL, save:: alpha_t(iim + 1, jjm + 1)
20      REAL, save:: alpha_u(iim + 1, jjm + 1), alpha_v(iim + 1, jjm)
21    
22      private iim, jjm
23    
24  contains  contains
25    
# Line 18  contains Line 28  contains
28      ! 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
29      !  Parametres de controle du run:      !  Parametres de controle du run:
30    
31      use getparam, only: ini_getparam, getpar, fin_getparam      use comconst, only: daysec, dtvr
32      use tau2alpha_m, only: lat_max_guide, lat_min_guide      use conf_gcm_m, only: day_step, iperiod
33        use dynetat0_m, only: grossismx, grossismy, rlatu, rlatv
34        use init_tau2alpha_m, only: init_tau2alpha
35        use nr_util, only: assert, pi
36        use tau2alpha_m, only: tau2alpha
37        use unit_nml_m, only: unit_nml
38        use writefield_m, only: writefield
39    
40        ! Local:
41    
42        ! Constantes de rappel, en jours :
43        REAL tau_min_u, tau_max_u, tau_min_v, tau_max_v, tau_min_t, tau_max_t
44        real tau_min_q, tau_max_q
45    
46        REAL lat_min_guide_deg, lat_max_guide_deg ! in degrees
47    
48        ! Dans le cas où on n'a les analyses que sur une bande de latitudes :
49        REAL lat_min_guide ! minimum latitude for nudging, in rad
50        real lat_max_guide ! maximum latitude for nudging, in rad
51    
52        REAL factt ! pas de temps entre deux appels au guidage, en jours
53        REAL dxdys(iim + 1, jjm + 1), dxdyu(iim + 1, jjm + 1), dxdyv(iim + 1, jjm)
54        integer i, j
55    
56        namelist /conf_guide_nml/ ini_anal, guide_u, guide_v, guide_t, guide_q, &
57             tau_min_u, tau_max_u, tau_min_v, tau_max_v, tau_min_t, tau_max_t, &
58             tau_min_q, tau_max_q, lat_min_guide_deg, lat_max_guide_deg
59    
60      !-----------------------------------------------------------------------      !-----------------------------------------------------------------------
61    
62      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 ')  
63    
64      CALL fin_getparam      ! Default values:
65        tau_min_u = 0.03
66        tau_max_u = 10.
67        tau_min_v = 0.03
68        tau_max_v = 10.
69        tau_min_t = 0.03
70        tau_max_t = 10.
71        tau_min_q = 0.03
72        tau_max_q = 10.
73        lat_min_guide_deg = -90.
74        lat_max_guide_deg = 90.
75    
76        print *, "Enter namelist 'conf_guide_nml'."
77        read(unit=*, nml=conf_guide_nml)
78        write(unit_nml, nml=conf_guide_nml)
79    
80        ok_guide = any((/guide_u, guide_v, guide_t, guide_q/))
81    
82        if (ok_guide) then
83           call assert(mod(day_step, 4 * iperiod) == 0, &
84                "conf_guide ok_guide day_step iperiod")
85           lat_min_guide = lat_min_guide_deg / 180. * pi
86           lat_max_guide = lat_max_guide_deg / 180. * pi
87           factt = dtvr * iperiod / daysec
88           print *, "factt = ", factt
89    
90           IF (abs(grossismx - 1.) < 0.1 .OR. abs(grossismy - 1.) < 0.1) THEN
91              ! grille regulière
92              if (guide_u) alpha_u = 1. - exp(- factt / tau_max_u)
93              if (guide_v) alpha_v = 1. - exp(- factt / tau_max_v)
94              if (guide_t) alpha_t = 1. - exp(- factt / tau_max_t)
95              if (guide_q) alpha_q = 1. - exp(- factt / tau_max_q)
96           else
97              call init_tau2alpha(dxdys)
98    
99              if (guide_u) then
100                 DO j = 1, jjm + 1
101                    DO i = 1, iim
102                       dxdyu(i, j) = 0.5 * (dxdys(i, j) + dxdys(i + 1, j))
103                    END DO
104                    dxdyu(iim + 1, j) = dxdyu(1, j)
105                 END DO
106    
107                 CALL tau2alpha(lat_min_guide, lat_max_guide, factt, dxdyu, rlatu, &
108                      tau_min_u, tau_max_u, alpha_u)
109                 CALL writefield("alpha_u", alpha_u)
110              end if
111    
112              if (guide_v) then
113                 DO j = 1, jjm
114                    DO i = 1, iim + 1
115                       dxdyv(i, j) = 0.5 * (dxdys(i, j) + dxdys(i, j + 1))
116                    END DO
117                 END DO
118    
119                 CALL tau2alpha(lat_min_guide, lat_max_guide, factt, dxdyv, rlatv, &
120                      tau_min_v, tau_max_v, alpha_v)
121                 CALL writefield("alpha_v", alpha_v)
122              end if
123    
124              if (guide_t) then
125                 CALL tau2alpha(lat_min_guide, lat_max_guide, factt, dxdys, rlatu, &
126                      tau_min_t, tau_max_t, alpha_t)
127                 CALL writefield("alpha_t", alpha_t)
128              end if
129    
130              if (guide_q)  then
131                 CALL tau2alpha(lat_min_guide, lat_max_guide, factt, dxdys, rlatu, &
132                      tau_min_q, tau_max_q, alpha_q)
133                 CALL writefield("alpha_q", alpha_q)
134              end if
135           end IF
136        end if
137    
138    end SUBROUTINE conf_guide    end SUBROUTINE conf_guide
139    
Legend:
Removed from v.76  
changed lines
  Added in v.221
  ViewVC Help
Powered by ViewVC 1.1.21