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

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

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

revision 178 by guez, Fri Mar 11 18:47:26 2016 UTC revision 210 by guez, Tue Dec 13 16:02:23 2016 UTC
# Line 1  Line 1 
1  MODULE guide_m  MODULE guide_m
2    
3    ! From dyn3d/guide.F, version 1.3 2005/05/25 13:10:09    ! From dyn3d/guide.F, version 1.3, 2005/05/25 13:10:09
4    ! and dyn3d/guide.h, version 1.1.1.1 2004/05/19 12:53:06    ! and dyn3d/guide.h, version 1.1.1.1, 2004/05/19 12:53:06
5    
6    IMPLICIT NONE    IMPLICIT NONE
7    
# Line 9  CONTAINS Line 9  CONTAINS
9    
10    SUBROUTINE guide(itau, ucov, vcov, teta, q, ps)    SUBROUTINE guide(itau, ucov, vcov, teta, q, ps)
11    
12      ! Author: F.Hourdin      ! Author: F. Hourdin
13    
14      USE comconst, ONLY: cpp, kappa      USE comconst, ONLY: cpp, kappa
15      USE conf_gcm_m, ONLY: day_step      USE conf_gcm_m, ONLY: day_step
# Line 21  CONTAINS Line 21  CONTAINS
21      use dynetat0_m, only: grossismx, grossismy, rlatu, rlatv      use dynetat0_m, only: grossismx, grossismy, rlatu, rlatv
22      USE exner_hyb_m, ONLY: exner_hyb      USE exner_hyb_m, ONLY: exner_hyb
23      use init_tau2alpha_m, only: init_tau2alpha      use init_tau2alpha_m, only: init_tau2alpha
24      USE paramet_m, ONLY: iip1, jjp1, llmp1      USE paramet_m, ONLY: iip1, jjp1
25      USE q_sat_m, ONLY: q_sat      USE q_sat_m, ONLY: q_sat
26      use read_reanalyse_m, only: read_reanalyse      use read_reanalyse_m, only: read_reanalyse
27      use tau2alpha_m, only: tau2alpha      use tau2alpha_m, only: tau2alpha
# Line 39  CONTAINS Line 39  CONTAINS
39    
40      ! Local:      ! Local:
41    
42      ! variables dynamiques pour les réanalyses      ! Variables dynamiques pour les réanalyses
43    
44      REAL, save:: ucovrea1(iim + 1, jjm + 1, llm), vcovrea1(iim + 1, jjm, llm)      REAL, save:: ucovrea1(iim + 1, jjm + 1, llm), vcovrea1(iim + 1, jjm, llm)
45      ! vents covariants reanalyses      ! vents covariants r\'eanalyses
46    
47      REAL, save:: tetarea1(iim + 1, jjm + 1, llm) ! temp pot reales      REAL, save:: tetarea1(iim + 1, jjm + 1, llm)
48      REAL, save:: qrea1(iim + 1, jjm + 1, llm) ! temp pot reales      ! potential temperture from reanalysis
49        
50        REAL, save:: qrea1(iim + 1, jjm + 1, llm)
51    
52      REAL, save:: ucovrea2(iim + 1, jjm + 1, llm), vcovrea2(iim + 1, jjm, llm)      REAL, save:: ucovrea2(iim + 1, jjm + 1, llm), vcovrea2(iim + 1, jjm, llm)
53      ! vents covariants reanalyses      ! vents covariants reanalyses
54    
55      REAL, save:: tetarea2(iim + 1, jjm + 1, llm) ! temp pot reales      REAL, save:: tetarea2(iim + 1, jjm + 1, llm)
56      REAL, save:: qrea2(iim + 1, jjm + 1, llm) ! temp pot reales      ! potential temperture from reanalysis
57        
58        REAL, save:: qrea2(iim + 1, jjm + 1, llm)
59    
60      ! alpha détermine la part des injections de données à chaque étape      ! alpha détermine la part des injections de données à chaque étape
61      ! alpha=0 signifie pas d'injection      ! alpha=0 signifie pas d'injection
# Line 64  CONTAINS Line 68  CONTAINS
68      REAL tau      REAL tau
69    
70      ! TEST SUR QSAT      ! TEST SUR QSAT
71      REAL p(iim + 1, jjm + 1, llmp1)      REAL p(iim + 1, jjm + 1, llm + 1)
72      real pk(iim + 1, jjm + 1, llm), pks(iim + 1, jjm + 1)      real pk(iim + 1, jjm + 1, llm), pks(iim + 1, jjm + 1)
73      REAL qsat(iim + 1, jjm + 1, llm)      REAL qsat(iim + 1, jjm + 1, llm)
74    
# Line 76  CONTAINS Line 80  CONTAINS
80         IF (online) THEN         IF (online) THEN
81            IF (abs(grossismx - 1.) < 0.1 .OR. abs(grossismy - 1.) < 0.1) THEN            IF (abs(grossismx - 1.) < 0.1 .OR. abs(grossismy - 1.) < 0.1) THEN
82               ! grille regulière               ! grille regulière
83               if (guide_u) alpha_u = factt / tau_max_u               if (guide_u) alpha_u = 1. - exp(- factt / tau_max_u)
84               if (guide_v) alpha_v = factt / tau_max_v               if (guide_v) alpha_v = 1. - exp(- factt / tau_max_v)
85               if (guide_t) alpha_t = factt / tau_max_t               if (guide_t) alpha_t = 1. - exp(- factt / tau_max_t)
86               if (guide_q) alpha_q = factt / tau_max_q               if (guide_q) alpha_q = 1. - exp(- factt / tau_max_q)
87            else            else
88               call init_tau2alpha(dxdys, dxdyu, dxdyv)               call init_tau2alpha(dxdys, dxdyu, dxdyv)
89    
Legend:
Removed from v.178  
changed lines
  Added in v.210
  ViewVC Help
Powered by ViewVC 1.1.21