- Timestamp:
- Jul 24, 2024, 1:27:51 PM (8 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/initfluxsto.f90
r5113 r5114 1 2 1 ! $Id$ 3 2 4 3 SUBROUTINE initfluxsto & 5 (infile, tstep,t_ops,t_wrt, &6 fileid, filevid,filedid)7 8 4 (infile, tstep, t_ops, t_wrt, & 5 fileid, filevid, filedid) 6 7 USE IOIPSL 9 8 USE comconst_mod, ONLY: pi 10 9 USE comvert_mod, ONLY: nivsigs 11 10 USE temps_mod, ONLY: annee_ref, day_ref, itau_dyn 11 USE lmdz_description, ONLY: descript 12 12 13 13 IMPLICIT NONE … … 43 43 include "paramet.h" 44 44 include "comgeom.h" 45 include "description.h"46 45 include "iniprint.h" 47 46 48 47 ! Arguments 49 48 ! 50 character(len =*) :: infile49 character(len = *) :: infile 51 50 real :: tstep, t_ops, t_wrt 52 integer :: fileid, filevid, filedid51 integer :: fileid, filevid, filedid 53 52 54 53 ! This routine needs IOIPSL to work … … 58 57 integer :: tau0 59 58 real :: zjulian 60 character(len =3) :: str61 character(len =10) :: ctrac59 character(len = 3) :: str 60 character(len = 10) :: ctrac 62 61 integer :: iq 63 real :: rlong(iip1, jjp1), rlat(iip1,jjp1),rl(1,1)64 integer :: uhoriid, vhoriid, thoriid, zvertiid, dhoriid,dvertiid65 integer :: ii, jj62 real :: rlong(iip1, jjp1), rlat(iip1, jjp1), rl(1, 1) 63 integer :: uhoriid, vhoriid, thoriid, zvertiid, dhoriid, dvertiid 64 integer :: ii, jj 66 65 integer :: zan, idayref 67 66 logical :: ok_sync … … 70 69 ! 71 70 pi = 4. * atan (1.) 72 str ='q '71 str = 'q ' 73 72 ctrac = 'traceur ' 74 73 ok_sync = .TRUE. … … 82 81 tau0 = itau_dyn 83 82 84 83 do jj = 1, jjp1 85 84 do ii = 1, iip1 86 rlong(ii, jj) = rlonu(ii) * 180. / pi87 rlat(ii, jj) = rlatu(jj) * 180. / pi85 rlong(ii, jj) = rlonu(ii) * 180. / pi 86 rlat(ii, jj) = rlatu(jj) * 180. / pi 88 87 enddo 89 88 enddo 90 89 91 CALL histbeg(infile, iip1, rlong(:, 1), jjp1, rlat(1,:), &92 1, iip1, 1, jjp1, &93 tau0, zjulian, tstep, uhoriid, fileid)90 CALL histbeg(infile, iip1, rlong(:, 1), jjp1, rlat(1, :), & 91 1, iip1, 1, jjp1, & 92 tau0, zjulian, tstep, uhoriid, fileid) 94 93 ! 95 94 ! Creation du fichier histoire pour la grille en V (oblige pour l'instant, … … 97 96 ! un meme fichier) 98 97 99 100 98 do jj = 1, jjm 101 99 do ii = 1, iip1 102 rlong(ii, jj) = rlonv(ii) * 180. / pi103 rlat(ii, jj) = rlatv(jj) * 180. / pi100 rlong(ii, jj) = rlonv(ii) * 180. / pi 101 rlat(ii, jj) = rlatv(jj) * 180. / pi 104 102 enddo 105 103 enddo 106 104 107 CALL histbeg('fluxstokev.nc', iip1, rlong(:, 1), jjm, rlat(1,:), &108 1, iip1, 1, jjm, &109 tau0, zjulian, tstep, vhoriid, filevid)110 111 rl(1,1) = 1.105 CALL histbeg('fluxstokev.nc', iip1, rlong(:, 1), jjm, rlat(1, :), & 106 1, iip1, 1, jjm, & 107 tau0, zjulian, tstep, vhoriid, filevid) 108 109 rl(1, 1) = 1. 112 110 CALL histbeg('defstoke.nc', 1, rl, 1, rl, & 113 1, 1, 1, 1, &114 tau0, zjulian, tstep, dhoriid, filedid)111 1, 1, 1, 1, & 112 tau0, zjulian, tstep, dhoriid, filedid) 115 113 116 114 ! … … 119 117 do jj = 1, jjp1 120 118 do ii = 1, iip1 121 rlong(ii, jj) = rlonv(ii) * 180. / pi122 rlat(ii, jj) = rlatu(jj) * 180. / pi119 rlong(ii, jj) = rlonv(ii) * 180. / pi 120 rlat(ii, jj) = rlatu(jj) * 180. / pi 123 121 enddo 124 122 enddo 125 123 126 124 CALL histhori(fileid, iip1, rlong, jjp1, rlat, 'scalar', & 127 'Grille points scalaires', thoriid)125 'Grille points scalaires', thoriid) 128 126 129 127 ! … … 131 129 ! 132 130 CALL histvert(fileid, 'sig_s', 'Niveaux sigma', & 133 'sigma_level', &134 llm, nivsigs, zvertiid)131 'sigma_level', & 132 llm, nivsigs, zvertiid) 135 133 ! Pour le fichier V 136 134 CALL histvert(filevid, 'sig_s', 'Niveaux sigma', & 137 'sigma_level', &138 llm, nivsigs, zvertiid)135 'sigma_level', & 136 llm, nivsigs, zvertiid) 139 137 ! pour le fichier def 140 138 nivd(1) = 1 141 139 CALL histvert(filedid, 'sig_s', 'Niveaux sigma', & 142 'sigma_level', &143 1, nivd, dvertiid)140 'sigma_level', & 141 1, nivd, dvertiid) 144 142 145 143 ! 146 144 ! Appels a histdef pour la definition des variables a sauvegarder 147 145 148 149 iip1, jjp1,thoriid, 1,1,1, -99, 32, &150 "once", t_ops, t_wrt) 151 152 153 iip1,jjp1,thoriid, 1,1,1, -99, 32, &154 155 156 157 1, 1,dhoriid, 1,1,1, -99, 32, &158 "once", t_ops, t_wrt) 159 160 161 1,1,dhoriid, 1,1,1, -99, 32, &162 163 164 165 1,1,dhoriid, 1,1,1, -99, 32, &166 146 CALL histdef(fileid, "phis", "Surface geop. height", "-", & 147 iip1, jjp1, thoriid, 1, 1, 1, -99, 32, & 148 "once", t_ops, t_wrt) 149 150 CALL histdef(fileid, "aire", "Grid area", "-", & 151 iip1, jjp1, thoriid, 1, 1, 1, -99, 32, & 152 "once", t_ops, t_wrt) 153 154 CALL histdef(filedid, "dtvr", "tps dyn", "s", & 155 1, 1, dhoriid, 1, 1, 1, -99, 32, & 156 "once", t_ops, t_wrt) 157 158 CALL histdef(filedid, "istdyn", "tps stock", "s", & 159 1, 1, dhoriid, 1, 1, 1, -99, 32, & 160 "once", t_ops, t_wrt) 161 162 CALL histdef(filedid, "istphy", "tps stock phy", "s", & 163 1, 1, dhoriid, 1, 1, 1, -99, 32, & 164 "once", t_ops, t_wrt) 167 165 168 166 … … 171 169 ! 172 170 CALL histdef(fileid, 'masse', 'Masse', 'kg', & 173 iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &174 32, 'inst(X)', t_ops, t_wrt)171 iip1, jjp1, thoriid, llm, 1, llm, zvertiid, & 172 32, 'inst(X)', t_ops, t_wrt) 175 173 ! 176 174 ! Pbaru 177 175 ! 178 176 CALL histdef(fileid, 'pbaru', 'flx de masse zonal', 'kg m/s', & 179 iip1, jjp1, uhoriid, llm, 1, llm, zvertiid, &180 32, 'inst(X)', t_ops, t_wrt)177 iip1, jjp1, uhoriid, llm, 1, llm, zvertiid, & 178 32, 'inst(X)', t_ops, t_wrt) 181 179 182 180 ! … … 184 182 ! 185 183 CALL histdef(filevid, 'pbarv', 'flx de masse mer', 'kg m/s', & 186 iip1, jjm, vhoriid, llm, 1, llm, zvertiid, &187 32, 'inst(X)', t_ops, t_wrt)184 iip1, jjm, vhoriid, llm, 1, llm, zvertiid, & 185 32, 'inst(X)', t_ops, t_wrt) 188 186 ! 189 187 ! w 190 188 ! 191 189 CALL histdef(fileid, 'w', 'flx de masse vert', 'kg m/s', & 192 iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &193 32, 'inst(X)', t_ops, t_wrt)190 iip1, jjp1, thoriid, llm, 1, llm, zvertiid, & 191 32, 'inst(X)', t_ops, t_wrt) 194 192 195 193 ! … … 197 195 ! 198 196 CALL histdef(fileid, 'teta', 'temperature potentielle', '-', & 199 iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &200 32, 'inst(X)', t_ops, t_wrt)197 iip1, jjp1, thoriid, llm, 1, llm, zvertiid, & 198 32, 'inst(X)', t_ops, t_wrt) 201 199 ! 202 200 … … 205 203 ! 206 204 CALL histdef(fileid, 'phi', 'geopotentiel instantane', '-', & 207 iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &208 32, 'inst(X)', t_ops, t_wrt)205 iip1, jjp1, thoriid, llm, 1, llm, zvertiid, & 206 32, 'inst(X)', t_ops, t_wrt) 209 207 ! 210 208 ! Fin
Note: See TracChangeset
for help on using the changeset viewer.