Changeset 2408 for LMDZ5/branches/testing/libf/phylmd/initphysto.F90
- Timestamp:
- Dec 14, 2015, 11:43:09 AM (9 years ago)
- Location:
- LMDZ5/branches/testing
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/testing
- Property svn:mergeinfo changed
/LMDZ5/trunk merged: 2293-2295,2297,2299-2302,2305-2313,2315,2317-2380,2382-2396
- Property svn:mergeinfo changed
-
LMDZ5/branches/testing/libf/phylmd/initphysto.F90
r1910 r2408 8 8 USE IOIPSL 9 9 USE iophy 10 USE control_mod11 10 USE indice_sol_mod 11 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, nbp_lev 12 USE time_phylmdz_mod, ONLY: day_ref, annee_ref 12 13 13 14 IMPLICIT NONE … … 37 38 ! 38 39 ! ===================================================================== 39 !40 ! Declarations41 INCLUDE "dimensions.h"42 INCLUDE "paramet.h"43 INCLUDE "comconst.h"44 INCLUDE "comgeom.h"45 INCLUDE "temps.h"46 INCLUDE "logic.h"47 INCLUDE "description.h"48 INCLUDE "serre.h"49 40 50 41 ! Arguments … … 58 49 INTEGER nhoriid, i 59 50 INTEGER l,k 60 REAL nivsigs( llm)51 REAL nivsigs(nbp_lev) 61 52 INTEGER tau0 62 53 REAL zjulian … … 66 57 INTEGER zan, idayref 67 58 LOGICAL ok_sync 68 REAL zx_lon( iim,jjm+1), zx_lat(iim,jjm+1)59 REAL zx_lon(nbp_lon,nbp_lat), zx_lat(nbp_lon,nbp_lat) 69 60 CHARACTER(len=12) :: nvar 70 61 71 62 ! Initialisations 72 63 ! 73 pi = 4. * ATAN (1.)74 64 ok_sync= .TRUE. 75 65 ! … … 88 78 ! Appel a histvert pour la grille verticale 89 79 ! 90 DO l=1, llm80 DO l=1,nbp_lev 91 81 nivsigs(l)=REAL(l) 92 82 ENDDO … … 94 84 CALL histvert(fileid, 'sig_s', 'Niveaux sigma', & 95 85 'sigma_level', & 96 llm, nivsigs, zvertiid)86 nbp_lev, nivsigs, zvertiid) 97 87 ! 98 88 ! Appels a histdef pour la definition des variables a sauvegarder 99 89 ! 100 90 CALL histdef(fileid, "phis", "Surface geop. height", "-", & 101 iim,jj_nb,nhoriid, 1,1,1, -99, 32, &91 nbp_lon,jj_nb,nhoriid, 1,1,1, -99, 32, & 102 92 "once", t_ops, t_wrt) 103 93 104 94 CALL histdef(fileid, "aire", "Grid area", "-", & 105 iim,jj_nb,nhoriid, 1,1,1, -99, 32, &95 nbp_lon,jj_nb,nhoriid, 1,1,1, -99, 32, & 106 96 "once", t_ops, t_wrt) 107 97 108 98 CALL histdef(fileid, "longitudes", "longitudes", "-", & 109 iim,jj_nb,nhoriid, 1,1,1, -99, 32, &99 nbp_lon,jj_nb,nhoriid, 1,1,1, -99, 32, & 110 100 "once", t_ops, t_wrt) 111 101 112 102 CALL histdef(fileid, "latitudes", "latitudes", "-", & 113 iim,jj_nb,nhoriid, 1,1,1, -99, 32, &103 nbp_lon,jj_nb,nhoriid, 1,1,1, -99, 32, & 114 104 "once", t_ops, t_wrt) 115 105 ! T 116 CALL histdef(fileid, 't', 'Temperature', 'K', iim, jj_nb, nhoriid, &117 llm, 1, llm, zvertiid, 32, 'inst(X)', t_ops, t_wrt)106 CALL histdef(fileid, 't', 'Temperature', 'K', nbp_lon, jj_nb, nhoriid, & 107 nbp_lev, 1, nbp_lev, zvertiid, 32, 'inst(X)', t_ops, t_wrt) 118 108 ! mfu 119 CALL histdef(fileid, 'mfu', 'flx m. pan. mt', 'kg m/s', iim, jj_nb, nhoriid, &120 llm, 1, llm, zvertiid,32, 'inst(X)', t_ops, t_wrt)109 CALL histdef(fileid, 'mfu', 'flx m. pan. mt', 'kg m/s',nbp_lon, jj_nb, nhoriid, & 110 nbp_lev, 1, nbp_lev, zvertiid,32, 'inst(X)', t_ops, t_wrt) 121 111 ! mfd 122 CALL histdef(fileid, 'mfd', 'flx m. pan. des', 'kg m/s', iim, jj_nb, nhoriid, &123 llm, 1, llm, zvertiid, 32, 'inst(X)', t_ops, t_wrt)112 CALL histdef(fileid, 'mfd', 'flx m. pan. des', 'kg m/s',nbp_lon, jj_nb, nhoriid, & 113 nbp_lev, 1, nbp_lev, zvertiid, 32, 'inst(X)', t_ops, t_wrt) 124 114 ! en_u 125 CALL histdef(fileid, 'en_u', 'flx ent pan mt', 'kg m/s', iim, jj_nb, nhoriid, &126 llm, 1, llm, zvertiid,32, 'inst(X)', t_ops, t_wrt)115 CALL histdef(fileid, 'en_u', 'flx ent pan mt', 'kg m/s', nbp_lon, jj_nb, nhoriid, & 116 nbp_lev, 1, nbp_lev, zvertiid,32, 'inst(X)', t_ops, t_wrt) 127 117 ! de_u 128 CALL histdef(fileid, 'de_u', 'flx det pan mt', 'kg m/s', iim, jj_nb, nhoriid, &129 llm, 1, llm, zvertiid,32, 'inst(X)', t_ops, t_wrt)118 CALL histdef(fileid, 'de_u', 'flx det pan mt', 'kg m/s',nbp_lon, jj_nb, nhoriid, & 119 nbp_lev, 1, nbp_lev, zvertiid,32, 'inst(X)', t_ops, t_wrt) 130 120 ! en_d 131 CALL histdef(fileid, 'en_d', 'flx ent pan dt', 'kg m/s', iim, jj_nb, nhoriid, &132 llm, 1, llm, zvertiid,32, 'inst(X)', t_ops, t_wrt)121 CALL histdef(fileid, 'en_d', 'flx ent pan dt', 'kg m/s', nbp_lon, jj_nb, nhoriid, & 122 nbp_lev, 1, nbp_lev, zvertiid,32, 'inst(X)', t_ops, t_wrt) 133 123 ! de_d 134 CALL histdef(fileid, 'de_d', 'flx det pan dt', 'kg m/s', iim, jj_nb, nhoriid, &135 llm, 1, llm, zvertiid, 32, 'inst(X)', t_ops, t_wrt)124 CALL histdef(fileid, 'de_d', 'flx det pan dt', 'kg m/s', nbp_lon, jj_nb, nhoriid, & 125 nbp_lev, 1, nbp_lev, zvertiid, 32, 'inst(X)', t_ops, t_wrt) 136 126 ! coefh 137 CALL histdef(fileid, "coefh", " ", " ", iim, jj_nb, nhoriid, &138 llm, 1, llm, zvertiid,32, "inst(X)", t_ops, t_wrt)127 CALL histdef(fileid, "coefh", " ", " ", nbp_lon, jj_nb, nhoriid, & 128 nbp_lev, 1, nbp_lev, zvertiid,32, "inst(X)", t_ops, t_wrt) 139 129 ! fm_th 140 CALL histdef(fileid, "fm_th", " ", " ", iim, jj_nb, nhoriid, &141 llm, 1, llm, zvertiid,32, "inst(X)", t_ops, t_wrt)130 CALL histdef(fileid, "fm_th", " ", " ",nbp_lon, jj_nb, nhoriid, & 131 nbp_lev, 1, nbp_lev, zvertiid,32, "inst(X)", t_ops, t_wrt) 142 132 ! en_th 143 CALL histdef(fileid, "en_th", " ", " ", iim, jj_nb, nhoriid, &144 llm, 1, llm, zvertiid,32, "inst(X)", t_ops, t_wrt)133 CALL histdef(fileid, "en_th", " ", " ",nbp_lon, jj_nb, nhoriid, & 134 nbp_lev, 1, nbp_lev, zvertiid,32, "inst(X)", t_ops, t_wrt) 145 135 ! frac_impa 146 CALL histdef(fileid, 'frac_impa', ' ', ' ', iim, jj_nb, nhoriid, &147 llm, 1, llm, zvertiid,32, 'inst(X)', t_ops, t_wrt)136 CALL histdef(fileid, 'frac_impa', ' ', ' ',nbp_lon, jj_nb, nhoriid, & 137 nbp_lev, 1, nbp_lev, zvertiid,32, 'inst(X)', t_ops, t_wrt) 148 138 ! frac_nucl 149 CALL histdef(fileid, 'frac_nucl', ' ', ' ', iim, jj_nb, nhoriid, &150 llm, 1, llm, zvertiid,32, 'inst(X)', t_ops, t_wrt)139 CALL histdef(fileid, 'frac_nucl', ' ', ' ',nbp_lon, jj_nb, nhoriid, & 140 nbp_lev, 1, nbp_lev, zvertiid,32, 'inst(X)', t_ops, t_wrt) 151 141 ! pyu1 152 CALL histdef(fileid, "pyu1", " ", " ", iim,jj_nb,nhoriid, &142 CALL histdef(fileid, "pyu1", " ", " ", nbp_lon,jj_nb,nhoriid, & 153 143 1,1,1, -99, 32, "inst(X)", t_ops, t_wrt) 154 144 ! pyv1 155 CALL histdef(fileid, "pyv1", " ", " ", iim,jj_nb,nhoriid, &145 CALL histdef(fileid, "pyv1", " ", " ", nbp_lon,jj_nb,nhoriid, & 156 146 1,1,1, -99, 32,"inst(X)", t_ops, t_wrt) 157 147 ! ftsol1 158 CALL histdef(fileid, "ftsol1", " ", " ", iim, jj_nb, nhoriid, &148 CALL histdef(fileid, "ftsol1", " ", " ",nbp_lon, jj_nb, nhoriid, & 159 149 1, 1,1, -99,32, "inst(X)", t_ops, t_wrt) 160 150 ! ftsol2 161 CALL histdef(fileid, "ftsol2", " ", " ", iim, jj_nb, nhoriid, &151 CALL histdef(fileid, "ftsol2", " ", " ",nbp_lon, jj_nb, nhoriid, & 162 152 1, 1,1, -99,32, "inst(X)", t_ops, t_wrt) 163 153 ! ftsol3 164 CALL histdef(fileid, "ftsol3", " ", " ", iim, jj_nb, nhoriid, &154 CALL histdef(fileid, "ftsol3", " ", " ", nbp_lon, jj_nb, nhoriid, & 165 155 1, 1,1, -99,32, "inst(X)", t_ops, t_wrt) 166 156 ! ftsol4 167 CALL histdef(fileid, "ftsol4", " ", " ", iim, jj_nb, nhoriid, &157 CALL histdef(fileid, "ftsol4", " ", " ",nbp_lon, jj_nb, nhoriid, & 168 158 1, 1,1, -99, 32, "inst(X)", t_ops, t_wrt) 169 159 ! psrf1 170 CALL histdef(fileid, "psrf1", " ", " ", iim, jj_nb, nhoriid, &160 CALL histdef(fileid, "psrf1", " ", " ",nbp_lon, jj_nb, nhoriid, & 171 161 1, 1, 1, -99,32, "inst(X)", t_ops, t_wrt) 172 162 ! psrf2 173 CALL histdef(fileid, "psrf2", " ", " ", iim, jj_nb, nhoriid, &163 CALL histdef(fileid, "psrf2", " ", " ",nbp_lon, jj_nb, nhoriid, & 174 164 1, 1, 1, -99, 32, "inst(X)", t_ops, t_wrt) 175 165 ! psrf3 176 CALL histdef(fileid, "psrf3", " ", " ", iim, jj_nb, nhoriid, &166 CALL histdef(fileid, "psrf3", " ", " ",nbp_lon, jj_nb, nhoriid, & 177 167 1, 1, 1, -99, 32, "inst(X)", t_ops, t_wrt) 178 168 ! psrf4 179 CALL histdef(fileid, "psrf4", " ", " ", iim, jj_nb, nhoriid, &169 CALL histdef(fileid, "psrf4", " ", " ", nbp_lon, jj_nb, nhoriid, & 180 170 1, 1, 1, -99,32, "inst(X)", t_ops, t_wrt) 181 171 ! sh 182 CALL histdef(fileid, 'sh', '', '', iim, jj_nb, nhoriid, &183 llm, 1, llm, zvertiid, 32, 'inst(X)', t_ops, t_wrt)172 CALL histdef(fileid, 'sh', '', '', nbp_lon, jj_nb, nhoriid, & 173 nbp_lev, 1, nbp_lev, zvertiid, 32, 'inst(X)', t_ops, t_wrt) 184 174 ! da 185 CALL histdef(fileid, 'da', '', '', iim, jj_nb, nhoriid, &186 llm, 1, llm, zvertiid, 32, 'inst(X)', t_ops, t_wrt)175 CALL histdef(fileid, 'da', '', '', nbp_lon, jj_nb, nhoriid, & 176 nbp_lev, 1, nbp_lev, zvertiid, 32, 'inst(X)', t_ops, t_wrt) 187 177 ! mp 188 CALL histdef(fileid, 'mp', '', '', iim, jj_nb, nhoriid, &189 llm, 1, llm, zvertiid, 32, 'inst(X)', t_ops, t_wrt)178 CALL histdef(fileid, 'mp', '', '', nbp_lon, jj_nb, nhoriid, & 179 nbp_lev, 1, nbp_lev, zvertiid, 32, 'inst(X)', t_ops, t_wrt) 190 180 ! upwd 191 CALL histdef(fileid, 'upwd', '', '', iim, jj_nb, nhoriid, &192 llm, 1, llm, zvertiid, 32, 'inst(X)', t_ops, t_wrt)181 CALL histdef(fileid, 'upwd', '', '', nbp_lon, jj_nb, nhoriid, & 182 nbp_lev, 1, nbp_lev, zvertiid, 32, 'inst(X)', t_ops, t_wrt) 193 183 ! dnwd 194 CALL histdef(fileid, 'dnwd', '', '', iim, jj_nb, nhoriid, &195 llm, 1, llm, zvertiid, 32, 'inst(X)', t_ops, t_wrt)184 CALL histdef(fileid, 'dnwd', '', '', nbp_lon, jj_nb, nhoriid, & 185 nbp_lev, 1, nbp_lev, zvertiid, 32, 'inst(X)', t_ops, t_wrt) 196 186 197 187 ! phi 198 DO k=1, llm188 DO k=1,nbp_lev 199 189 IF (k<10) THEN 200 190 WRITE(nvar,'(i1)') k … … 206 196 nvar='phi_lev'//trim(nvar) 207 197 208 CALL histdef(fileid, nvar, '', '', iim, jj_nb, nhoriid, &209 llm, 1, llm, zvertiid, 32, 'inst(X)', t_ops, t_wrt)198 CALL histdef(fileid, nvar, '', '', nbp_lon, jj_nb, nhoriid, & 199 nbp_lev, 1, nbp_lev, zvertiid, 32, 'inst(X)', t_ops, t_wrt) 210 200 END DO 211 201
Note: See TracChangeset
for help on using the changeset viewer.