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

Annotation of /trunk/Sources/dyn3d/inithist.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 144 - (hide annotations)
Wed Jun 10 16:46:46 2015 UTC (8 years, 11 months ago) by guez
File size: 4109 byte(s) 
In procedure fxhyp, the convoluted computation of tanh(fa/fb) occurred
three times. Extracted it into a function. Also, the computation of
xmoy and fxm was repeated. So stored the values in arrays instead.

In procedure fxhyp, in the computation of fhyp, there were tests
xtild(i) == 0. and xtild(i) == pi_d. No use to do these tests at each
iteration. We now they are true for i == nmax and i == 2 * nmax,
respectively, and we know they are false for other values of
"i". Similarly, in the computations of ffdx and xxpr, there were the
tests xmoy == 0. and xmoy == pi_d, these could not be true.

Moved files from bibio to dyn3d, following LMDZ.
1 guez 3 module inithist_m
2    
3     implicit none
4    
5     contains
6    
7 guez 129 subroutine inithist(tstep, nq, t_ops, t_wrt)
8 guez 3
9 guez 39 ! From inithist.F, version 1.1.1.1 2004/05/19 12:53:05
10 guez 61 ! L. Fairhead, LMD, 03/99
11 guez 3
12 guez 61 ! Routine d'initialisation des écritures des fichiers histoires
13     ! LMDZ au format IOIPSL.
14 guez 3
15 guez 61 USE com_io_dyn, ONLY: histid, histuid, histvid
16 guez 92 USE dimens_m, ONLY: jjm, llm
17     USE disvert_m, ONLY: presnivs
18 guez 139 use dynetat0_m, only: day_ref, annee_ref, rlatu, rlatv, rlonu, rlonv
19 guez 61 USE histbeg_totreg_m, ONLY : histbeg_totreg
20     USE histdef_m, ONLY : histdef
21     USE histend_m, ONLY : histend
22     USE histvert_m, ONLY : histvert
23 guez 92 USE iniadvtrac_m, ONLY: ttext
24     USE nr_util, ONLY: pi
25 guez 61 USE paramet_m, ONLY: iip1, jjp1
26     USE temps, ONLY: itau_dyn
27 guez 92 USE ymds2ju_m, ONLY: ymds2ju
28 guez 3
29 guez 61 real, intent(in):: tstep ! durée du pas de temps en secondes
30     integer, intent(in):: nq ! nombre de traceurs
31     real, intent(in):: t_ops ! fréquence de l'opération pour IOIPSL
32     real, intent(in):: t_wrt ! fréquence d'écriture sur le fichier
33 guez 3
34 guez 61 ! Variables locales:
35 guez 3 real zjulian
36     integer iq
37 guez 39 real rlong(iip1, jjp1), rlat(iip1, jjp1)
38 guez 3 integer uhoriid, vhoriid, thoriid, zvertiid
39 guez 39 integer ii, jj
40 guez 3
41     !-----------------------------------------------------------------------
42    
43 guez 129 CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
44 guez 3
45     do jj = 1, jjp1
46     do ii = 1, iip1
47 guez 39 rlong(ii, jj) = rlonu(ii) * 180. / pi
48     rlat(ii, jj) = rlatu(jj) * 180. / pi
49 guez 3 enddo
50     enddo
51    
52 guez 56 call histbeg_totreg("dyn_histu.nc", rlong(:,1), rlat(1,:), 1, iip1, 1, &
53     jjp1, itau_dyn, zjulian, tstep, uhoriid, histuid)
54 guez 3
55     do jj = 1, jjm
56     do ii = 1, iip1
57 guez 39 rlong(ii, jj) = rlonv(ii) * 180. / pi
58     rlat(ii, jj) = rlatv(jj) * 180. / pi
59 guez 3 enddo
60     enddo
61    
62 guez 56 ! Creation du fichier histoire pour la grille en V (oblige pour l'instant,
63     ! IOIPSL ne permet pas de grilles avec des nombres de point differents dans
64     ! un meme fichier)
65    
66 guez 61 call histbeg_totreg('dyn_histv.nc', rlong(:, 1), rlat(1, :jjm), 1, iip1, &
67     1, jjm, itau_dyn, zjulian, tstep, vhoriid, histvid)
68    
69 guez 3 do jj = 1, jjp1
70     do ii = 1, iip1
71 guez 39 rlong(ii, jj) = rlonv(ii) * 180. / pi
72     rlat(ii, jj) = rlatu(jj) * 180. / pi
73 guez 3 enddo
74     enddo
75    
76 guez 56 call histbeg_totreg("dyn_hist.nc", rlong(:, 1), rlat(1, :), 1, iip1, 1, &
77     jjp1, itau_dyn, zjulian, tstep, thoriid, histid)
78 guez 61
79 guez 39 ! Appel a histvert pour la grille verticale
80 guez 61
81 guez 67 call histvert(histid, 'presnivs', 'Niveaux pression','mb', presnivs/100., &
82     zvertiid,'down')
83     call histvert(histvid, 'presnivs', 'Niveaux pression','mb', &
84 guez 61 presnivs/100., zvertiid,'down')
85 guez 67 call histvert(histuid, 'presnivs', 'Niveaux pression','mb', &
86 guez 61 presnivs/100., zvertiid,'down')
87    
88 guez 39 ! Appels a histdef pour la definition des variables a sauvegarder
89 guez 3
90 guez 61 call histdef(histuid, 'u', 'vent u', 'm/s', iip1, jjp1, uhoriid, llm, 1, &
91     llm, zvertiid, 'inst(X)', t_ops, t_wrt)
92     call histdef(histvid, 'v', 'vent v', 'm/s', iip1, jjm, vhoriid, llm, 1, &
93     llm, zvertiid, 'inst(X)', t_ops, t_wrt)
94     call histdef(histid, 'teta', 'temperature potentielle', '-', iip1, jjp1, &
95     thoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
96     call histdef(histid, 'phi', 'geopotentiel', '-', iip1, jjp1, thoriid, &
97     llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
98    
99 guez 39 ! Traceurs
100     DO iq=1, nq
101 guez 61 call histdef(histid, ttext(iq), ttext(iq), '-', iip1, jjp1, thoriid, &
102     llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
103 guez 3 enddo
104 guez 61
105     call histdef(histid, 'masse', 'masse', 'kg', iip1, jjp1, thoriid, llm, 1, &
106     llm, zvertiid, 'inst(X)', t_ops, t_wrt)
107     call histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', iip1, jjp1, &
108     thoriid, 1, 1, 1, -99, 'inst(X)', t_ops, t_wrt)
109     call histdef(histid, 'phis', 'geopotentiel au sol', '-', iip1, jjp1, &
110     thoriid, 1, 1, 1, -99, 'inst(X)', t_ops, t_wrt)
111    
112 guez 56 call histend(histid)
113     call histend(histuid)
114     call histend(histvid)
115 guez 3
116     end subroutine inithist
117    
118     end module inithist_m
  ViewVC Help
Powered by ViewVC 1.1.21