/[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.f revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC trunk/Sources/dyn3d/Guide/conf_guide.f revision 211 by guez, Tue Dec 13 17:23:09 2016 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 abort_gcm_m, only: abort_gcm
32      use tau2alpha_m, only: lat_max_guide, lat_min_guide      use comconst, only: daysec, dtvr
33        use conf_gcm_m, only: day_step, iperiod
34        use dynetat0_m, only: grossismx, grossismy, rlatu, rlatv
35        use init_tau2alpha_m, only: init_tau2alpha
36        use nr_util, only: assert, pi
37        use tau2alpha_m, only: tau2alpha
38        use unit_nml_m, only: unit_nml
39        use writefield_m, only: writefield
40    
41        ! Local:
42    
43        ! Constantes de rappel, en jours :
44        REAL tau_min_u, tau_max_u, tau_min_v, tau_max_v, tau_min_t, tau_max_t
45        real tau_min_q, tau_max_q
46    
47        REAL lat_min_guide_deg, lat_max_guide_deg ! in degrees
48    
49        ! Dans le cas où on n'a les analyses que sur une bande de latitudes :
50        REAL lat_min_guide ! minimum latitude for nudging, in rad
51        real lat_max_guide ! maximum latitude for nudging, in rad
52    
53        REAL factt ! pas de temps entre deux appels au guidage, en jours
54        REAL dxdys(iim + 1, jjm + 1), dxdyu(iim + 1, jjm + 1), dxdyv(iim + 1, jjm)
55        integer i, j
56    
57        namelist /conf_guide_nml/ ini_anal, guide_u, guide_v, guide_t, guide_q, &
58             tau_min_u, tau_max_u, tau_min_v, tau_max_v, tau_min_t, tau_max_t, &
59             tau_min_q, tau_max_q, lat_min_guide_deg, lat_max_guide_deg
60    
61      !-----------------------------------------------------------------------      !-----------------------------------------------------------------------
62    
63      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 ')  
64    
65      CALL fin_getparam      ! Default values:
66        tau_min_u = 0.03
67        tau_max_u = 10.
68        tau_min_v = 0.03
69        tau_max_v = 10.
70        tau_min_t = 0.03
71        tau_max_t = 10.
72        tau_min_q = 0.03
73        tau_max_q = 10.
74        lat_min_guide_deg = -90.
75        lat_max_guide_deg = 90.
76    
77        print *, "Enter namelist 'conf_guide_nml'."
78        read(unit=*, nml=conf_guide_nml)
79        write(unit_nml, nml=conf_guide_nml)
80    
81        ok_guide = any((/guide_u, guide_v, guide_t, guide_q/))
82    
83        if (ok_guide) then
84           call assert(mod(day_step, 4 * iperiod) == 0, &
85                "conf_guide ok_guide day_step iperiod")
86           lat_min_guide = lat_min_guide_deg / 180. * pi
87           lat_max_guide = lat_max_guide_deg / 180. * pi
88           factt = dtvr * iperiod / daysec
89           print *, "factt = ", factt
90    
91           IF (abs(grossismx - 1.) < 0.1 .OR. abs(grossismy - 1.) < 0.1) THEN
92              ! grille regulière
93              if (guide_u) alpha_u = 1. - exp(- factt / tau_max_u)
94              if (guide_v) alpha_v = 1. - exp(- factt / tau_max_v)
95              if (guide_t) alpha_t = 1. - exp(- factt / tau_max_t)
96              if (guide_q) alpha_q = 1. - exp(- factt / tau_max_q)
97           else
98              call init_tau2alpha(dxdys)
99    
100              if (guide_u) then
101                 DO j = 1, jjm + 1
102                    DO i = 1, iim
103                       dxdyu(i, j) = 0.5 * (dxdys(i, j) + dxdys(i + 1, j))
104                    END DO
105                    dxdyu(iim + 1, j) = dxdyu(1, j)
106                 END DO
107    
108                 CALL tau2alpha(lat_min_guide, lat_max_guide, factt, dxdyu, rlatu, &
109                      tau_min_u, tau_max_u, alpha_u)
110                 CALL writefield("alpha_u", alpha_u)
111              end if
112    
113              if (guide_v) then
114                 DO j = 1, jjm
115                    DO i = 1, iim + 1
116                       dxdyv(i, j) = 0.5 * (dxdys(i, j) + dxdys(i, j + 1))
117                    END DO
118                 END DO
119    
120                 CALL tau2alpha(lat_min_guide, lat_max_guide, factt, dxdyv, rlatv, &
121                      tau_min_v, tau_max_v, alpha_v)
122                 CALL writefield("alpha_v", alpha_v)
123              end if
124    
125              if (guide_t) then
126                 CALL tau2alpha(lat_min_guide, lat_max_guide, factt, dxdys, rlatu, &
127                      tau_min_t, tau_max_t, alpha_t)
128                 CALL writefield("alpha_t", alpha_t)
129              end if
130    
131              if (guide_q)  then
132                 CALL tau2alpha(lat_min_guide, lat_max_guide, factt, dxdys, rlatu, &
133                      tau_min_q, tau_max_q, alpha_q)
134                 CALL writefield("alpha_q", alpha_q)
135              end if
136           end IF
137        end if
138    
139    end SUBROUTINE conf_guide    end SUBROUTINE conf_guide
140    
Legend:
Removed from v.82  
changed lines
  Added in v.211
  ViewVC Help
Powered by ViewVC 1.1.21