/[lmdze]/trunk/Sources/dyn3d/initfluxsto.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/initfluxsto.f

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

trunk/libf/bibio/initfluxsto.f90 revision 31 by guez, Thu Apr 1 14:59:19 2010 UTC trunk/Sources/bibio/initfluxsto.f revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC
# Line 23  SUBROUTINE initfluxsto(tstep, t_ops, t_w Line 23  SUBROUTINE initfluxsto(tstep, t_ops, t_w
23    
24    !   L. Fairhead, 03/99                                                !   L. Fairhead, 03/99                                            
25    
   USE calendar  
   USE histcom  
   USE dimens_m  
   USE paramet_m  
26    USE comconst    USE comconst
   USE comvert  
   USE logic  
27    USE comgeom    USE comgeom
28      use conf_gcm_m
29      USE dimens_m
30      USE disvert_m
31      use dynetat0_m, only: day_ref, annee_ref
32      USE histbeg_totreg_m, ONLY : histbeg_totreg
33      USE histdef_m, ONLY : histdef
34      USE histend_m, ONLY : histend
35      use histhori_regular_m, only: histhori_regular
36      use histsync_m, only: histsync
37      USE histvert_m, ONLY : histvert
38      USE nr_util, ONLY : pi
39      USE paramet_m
40    USE serre    USE serre
41    USE temps, ONLY : annee_ref, day_ref, itau_dyn    USE temps, ONLY : itau_dyn
42    USE ener    use ymds2ju_m, only: ymds2ju
43    
44    IMPLICIT NONE    IMPLICIT NONE
45    
# Line 47  SUBROUTINE initfluxsto(tstep, t_ops, t_w Line 53  SUBROUTINE initfluxsto(tstep, t_ops, t_w
53    
54    !   Variables locales                                                      !   Variables locales                                                  
55    REAL zjulian    REAL zjulian
56    CHARACTER*3 str    CHARACTER(len=3) str
57    CHARACTER*10 ctrac    CHARACTER(len=10) ctrac
58    INTEGER iq    INTEGER iq
59    REAL rlong(iip1, jjp1), rlat(iip1, jjp1)    REAL rlong(iip1, jjp1), rlat(iip1, jjp1)
60    INTEGER uhoriid, vhoriid, thoriid, zvertiid, dhoriid, dvertiid    INTEGER uhoriid, vhoriid, thoriid, zvertiid, dhoriid, dvertiid
61    INTEGER ii, jj    INTEGER ii, jj, l
   INTEGER zan, idayref  
62    LOGICAL ok_sync    LOGICAL ok_sync
63    
64    !---------------------------------------------------------    !---------------------------------------------------------
65    
66    !  Initialisations                                                          !  Initialisations                                                      
   
   pi = 4.*atan(1.)  
67    str = 'q  '    str = 'q  '
68    ctrac = 'traceur   '    ctrac = 'traceur   '
69    ok_sync = .TRUE.    ok_sync = .TRUE.
70    
71    !  Appel a histbeg: creation du fichier netcdf et initialisations diverses    !  Appel a histbeg: creation du fichier netcdf et initialisations diverses
72    
73    zan = annee_ref    CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
   idayref = day_ref  
   CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)  
74    
75    DO jj = 1, jjp1    DO jj = 1, jjp1
76       DO ii = 1, iip1       DO ii = 1, iip1
# Line 112  SUBROUTINE initfluxsto(tstep, t_ops, t_w Line 113  SUBROUTINE initfluxsto(tstep, t_ops, t_w
113    
114    !  Appel a histvert pour la grille verticale                                !  Appel a histvert pour la grille verticale                            
115    
116    CALL histvert(fileid, 'sig_s', 'Niveaux sigma', 'sigma_level', llm, &    CALL histvert(fileid, 'sig_s', 'Niveaux sigma', 'sigma_level', &
117         nivsigs, zvertiid)         (/(real(l), l = 1, llm)/), zvertiid)
118    ! Pour le fichier V                                                        ! Pour le fichier V                                                    
119    CALL histvert(filevid, 'sig_s', 'Niveaux sigma', 'sigma_level', llm, &    CALL histvert(filevid, 'sig_s', 'Niveaux sigma', 'sigma_level', &
120         nivsigs, zvertiid)         (/(real(l), l = 1, llm)/), zvertiid)
121    ! pour le fichier def                                                      ! pour le fichier def                                                  
122    nivd(1) = 1    nivd(1) = 1
123    CALL histvert(filedid, 'sig_s', 'Niveaux sigma', 'sigma_level', 1, nivd, &    CALL histvert(filedid, 'sig_s', 'Niveaux sigma', 'sigma_level', nivd, &
124         dvertiid)         dvertiid)
125    
126    !  Appels a histdef pour la definition des variables a sauvegarder          !  Appels a histdef pour la definition des variables a sauvegarder      
Legend:
Removed from v.31  
changed lines
  Added in v.134
  ViewVC Help
Powered by ViewVC 1.1.21