Changeset 1530 for trunk/LMDZ.TITAN/libf
- Timestamp:
- Apr 6, 2016, 8:51:34 AM (9 years ago)
- Location:
- trunk/LMDZ.TITAN/libf/phytitan
- Files:
-
- 3 deleted
- 34 edited
-
aaam_bud.F (modified) (13 diffs)
-
ajsec.F (modified) (4 diffs)
-
ballon.F (modified) (13 diffs)
-
calchim.F (modified) (3 diffs)
-
clmain.F (modified) (12 diffs)
-
cltrac.F (modified) (1 diff)
-
coefkzmin.F (modified) (1 diff)
-
common_mod.F90 (modified) (1 diff)
-
conf_dat3d.F (deleted)
-
diagphy.F (modified) (2 diffs)
-
drag_noro.F (modified) (2 diffs)
-
grid_noro.F (modified) (5 diffs)
-
gwprofil.F (modified) (1 diff)
-
gwstress.F (modified) (1 diff)
-
hgardfou.F (modified) (3 diffs)
-
ini_histday.h (modified) (19 diffs)
-
ini_histins.h (modified) (17 diffs)
-
ini_histmth.h (modified) (21 diffs)
-
inicondens.F (modified) (3 diffs)
-
iniinterp_horiz.F (deleted)
-
interp_horiz.F (deleted)
-
iophy.F90 (modified) (18 diffs)
-
orodrag.F (modified) (1 diff)
-
orosetup.F (modified) (1 diff)
-
phyetat0.F90 (modified) (1 diff)
-
phyredem.F90 (modified) (1 diff)
-
physiq.F (modified) (1 diff)
-
radlwsw.F (modified) (1 diff)
-
soil.F (modified) (1 diff)
-
sources.F (modified) (1 diff)
-
sugwd.F (modified) (1 diff)
-
tgmdat_mod.F90 (modified) (1 diff)
-
ustarhb.F (modified) (1 diff)
-
vdif_kcay.F (modified) (1 diff)
-
yamada.F (modified) (1 diff)
-
yamada4.F (modified) (1 diff)
-
zenang.F (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.TITAN/libf/phytitan/aaam_bud.F
r102 r1530 7 7 c 8 8 use dimphy 9 use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat, klon_glo 9 10 implicit none 10 11 c====================================================================== … … 52 53 c =================== 53 54 c 54 c iim--common-I: Number of longitude intervals55 c jjm--common-I: Number of latitude intervals55 c nbp_lon--common-I: Number of longitude intervals 56 c nbp_lat-1--common-I: Number of latitude intervals 56 57 c klon-common-I: Number of points seen by the physics 57 c iim*(jjm-1)+2 for instance58 c nbp_lon*(nbp_lat-1-1)+2 for instance 58 59 c klev-common-I: Number of vertical layers 59 60 c====================================================================== … … 78 79 c====================================================================== 79 80 80 #include "dimensions.h"81 81 c 82 82 c ARGUMENTS … … 106 106 C PUT AAM QUANTITIES AT ZERO: 107 107 C 108 if( iim+1.gt.801.or.jjm+1.gt.401)then108 if(nbp_lon+1.gt.801.or.nbp_lat.gt.401)then 109 109 print *,' Pb de dimension dans aaam_bud' 110 110 stop … … 114 114 hadley=1.e18 115 115 hadday=1.e18*1.e7 116 dlat=xpi/float(jjm) 117 dlon=2.*xpi/float(iim) 116 IF (klon_glo.EQ.1) THEN 117 dlat=xpi 118 ELSE 119 dlat=xpi/float(nbp_lat-1) 120 ENDIF 121 dlon=2.*xpi/float(nbp_lon) 118 122 119 123 do iax=1,3 … … 140 144 zlat(1)=plat(l)*xpi/180. 141 145 142 do i=1, iim+1146 do i=1,nbp_lon+1 143 147 144 148 zs(i,1)=phis(l)/rg … … 154 158 155 159 156 do j = 2, jjm160 do j = 2,nbp_lat-1 157 161 158 162 C Values at Greenwich (Periodicity) 159 163 160 zs( iim+1,j)=phis(l+1)/rg161 ps( iim+1,j)=p(l+1,1)162 ssou( iim+1,j)=dragu(l+1)+liftu(l+1)163 ssov( iim+1,j)=dragv(l+1)+liftv(l+1)164 blsu( iim+1,j)=clu(l+1)165 blsv( iim+1,j)=clv(l+1)166 zlon( iim+1)=-plon(l+1)*xpi/180.164 zs(nbp_lon+1,j)=phis(l+1)/rg 165 ps(nbp_lon+1,j)=p(l+1,1) 166 ssou(nbp_lon+1,j)=dragu(l+1)+liftu(l+1) 167 ssov(nbp_lon+1,j)=dragv(l+1)+liftv(l+1) 168 blsu(nbp_lon+1,j)=clu(l+1) 169 blsv(nbp_lon+1,j)=clv(l+1) 170 zlon(nbp_lon+1)=-plon(l+1)*xpi/180. 167 171 zlat(j)=plat(l+1)*xpi/180. 168 172 169 ub( iim+1,j)=0.170 vb( iim+1,j)=0.173 ub(nbp_lon+1,j)=0. 174 vb(nbp_lon+1,j)=0. 171 175 do k=1,nlev 172 ub(iim+1,j)=ub(iim+1,j)+u(l+1,k)*(p(l+1,k)-p(l+1,k+1))/rg 173 vb(iim+1,j)=vb(iim+1,j)+v(l+1,k)*(p(l+1,k)-p(l+1,k+1))/rg 176 ub(nbp_lon+1,j)=ub(nbp_lon+1,j)+u(l+1,k)* 177 & (p(l+1,k)-p(l+1,k+1))/rg 178 vb(nbp_lon+1,j)=vb(nbp_lon+1,j)+v(l+1,k)* 179 & (p(l+1,k)-p(l+1,k+1))/rg 174 180 enddo 175 181 176 182 177 do i=1, iim183 do i=1,nbp_lon 178 184 179 185 l=l+1 … … 201 207 202 208 l=l+1 203 ub(1, jjm+1)=0.204 vb(1, jjm+1)=0.209 ub(1,nbp_lat)=0. 210 vb(1,nbp_lat)=0. 205 211 do k=1,nlev 206 ub(1, jjm+1)=ub(1,jjm+1)+u(l,k)*(p(l,k)-p(l,k+1))/rg207 vb(1, jjm+1)=vb(1,jjm+1)+v(l,k)*(p(l,k)-p(l,k+1))/rg208 enddo 209 zlat( jjm+1)=plat(l)*xpi/180.210 211 do i=1, iim+1212 zs(i, jjm+1)=phis(l)/rg213 ps(i, jjm+1)=p(l,1)214 ssou(i, jjm+1)=dragu(l)+liftu(l)215 ssov(i, jjm+1)=dragv(l)+liftv(l)216 blsu(i, jjm+1)=clu(l)217 blsv(i, jjm+1)=clv(l)218 ub(i, jjm+1)=ub(1,jjm+1)219 vb(i, jjm+1)=vb(1,jjm+1)212 ub(1,nbp_lat)=ub(1,nbp_lat)+u(l,k)*(p(l,k)-p(l,k+1))/rg 213 vb(1,nbp_lat)=vb(1,nbp_lat)+v(l,k)*(p(l,k)-p(l,k+1))/rg 214 enddo 215 zlat(nbp_lat)=plat(l)*xpi/180. 216 217 do i=1,nbp_lon+1 218 zs(i,nbp_lat)=phis(l)/rg 219 ps(i,nbp_lat)=p(l,1) 220 ssou(i,nbp_lat)=dragu(l)+liftu(l) 221 ssov(i,nbp_lat)=dragv(l)+liftv(l) 222 blsu(i,nbp_lat)=clu(l) 223 blsv(i,nbp_lat)=clv(l) 224 ub(i,nbp_lat)=ub(1,nbp_lat) 225 vb(i,nbp_lat)=vb(1,nbp_lat) 220 226 enddo 221 227 … … 223 229 C MOMENT ANGULAIRE 224 230 C 225 DO j=1, jjm226 DO i=1, iim231 DO j=1,nbp_lat-1 232 DO i=1,nbp_lon 227 233 228 234 raam(1)=raam(1)-rea**3*dlon*dlat*0.5* … … 261 267 C 262 268 263 DO j=1, jjm264 DO i=1, iim269 DO j=1,nbp_lat-1 270 DO i=1,nbp_lon 265 271 tmou(1)=tmou(1)-rea**2*dlon*0.5*sin(zlon(i)) 266 272 c *(zs(i,j)-zs(i,j+1)) … … 272 278 ENDDO 273 279 274 DO j=2, jjm275 DO i=1, iim280 DO j=2,nbp_lat-1 281 DO i=1,nbp_lon 276 282 tmou(1)=tmou(1)+rea**2*dlat*0.5*sin(zlat(j)) 277 283 c *(zs(i+1,j)-zs(i,j)) … … 288 294 C 289 295 l=1 290 DO j=2, jjm291 DO i=1, iim296 DO j=2,nbp_lat-1 297 DO i=1,nbp_lon 292 298 l=l+1 293 299 tsso(1)=tsso(1)-rea**3*cos(zlat(j))*dlon*dlat* … … 335 341 c00 format(F12.5,5(1x,F12.5)) 336 342 337 write(iam+1,*)((zs(i,j),i=1, iim),j=1,jjm+1)338 write(iam+1,*)((ps(i,j),i=1, iim),j=1,jjm+1)339 write(iam+1,*)((ub(i,j),i=1, iim),j=1,jjm+1)340 write(iam+1,*)((vb(i,j),i=1, iim),j=1,jjm+1)341 write(iam+1,*)((ssou(i,j),i=1, iim),j=1,jjm+1)342 write(iam+1,*)((ssov(i,j),i=1, iim),j=1,jjm+1)343 write(iam+1,*)((blsu(i,j),i=1, iim),j=1,jjm+1)344 write(iam+1,*)((blsv(i,j),i=1, iim),j=1,jjm+1)343 write(iam+1,*)((zs(i,j),i=1,nbp_lon),j=1,nbp_lat) 344 write(iam+1,*)((ps(i,j),i=1,nbp_lon),j=1,nbp_lat) 345 write(iam+1,*)((ub(i,j),i=1,nbp_lon),j=1,nbp_lat) 346 write(iam+1,*)((vb(i,j),i=1,nbp_lon),j=1,nbp_lat) 347 write(iam+1,*)((ssou(i,j),i=1,nbp_lon),j=1,nbp_lat) 348 write(iam+1,*)((ssov(i,j),i=1,nbp_lon),j=1,nbp_lat) 349 write(iam+1,*)((blsu(i,j),i=1,nbp_lon),j=1,nbp_lat) 350 write(iam+1,*)((blsv(i,j),i=1,nbp_lon),j=1,nbp_lat) 345 351 346 352 RETURN -
trunk/LMDZ.TITAN/libf/phytitan/ajsec.F
r1048 r1530 7 7 8 8 use dimphy 9 use mod_grid_phy_lmdz, only: nbp_lev 9 10 use cpdet_mod, only: t2tpot, tpot2t 10 11 IMPLICIT none … … 27 28 c d_qfi-----output-R-Incrementation des traceurs 28 29 c====================================================================== 29 #include "dimensions.h"30 30 #include "YOMCST.h" 31 31 REAL paprs(klon,klev+1), pplay(klon,klev) … … 66 66 c------------------------------------- passage en temperature potentielle 67 67 ! ADAPTATION GCM POUR CP(T) 68 call t2tpot(klon* llm,tfi,zh,ppk)68 call t2tpot(klon*nbp_lev,tfi,zh,ppk) 69 69 c 70 70 DO k = limbas, limhau … … 168 168 c------------------------------------- et calcul du d_t 169 169 ! ADAPTATION GCM POUR CP(T) 170 call tpot2t(klon* llm,zh,zt,ppk)170 call tpot2t(klon*nbp_lev,zh,zt,ppk) 171 171 172 172 DO k = limbas, limhau -
trunk/LMDZ.TITAN/libf/phytitan/ballon.F
r102 r1530 3 3 4 4 use dimphy 5 use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat 5 6 implicit none 6 7 … … 94 95 c====================================================================== 95 96 96 #include "dimensions.h"97 97 #include "YOMCST.h" 98 98 c … … 110 110 INTEGER jj,ii,ll 111 111 112 REAL zlon(iim+1),zlat(jjm+1) 113 save zlon,zlat 112 REAL,SAVE,ALLOCATABLE :: zlon(:),zlat(:) 114 113 115 114 REAL time … … 141 140 print*,"BALLOONS ACTIVATED" 142 141 142 allocate(zlon(nbp_lon+1)) 143 allocate(zlat(nbp_lat)) 144 143 145 C Latitudes: 144 146 zlat(1)=plat(1)*RPI/180. 145 do j = 2, jjm146 k=(j-2)* iim+2147 do j = 2,nbp_lat-1 148 k=(j-2)*nbp_lon+2 147 149 zlat(j)=plat(k)*RPI/180. 148 150 enddo 149 zlat( jjm+1)=plat(klon)*RPI/180.151 zlat(nbp_lat)=plat(klon)*RPI/180. 150 152 151 153 C Longitudes: 152 do i = 1, iim154 do i = 1,nbp_lon 153 155 k=i+1 154 156 zlon(i)=plon(k)*RPI/180. 155 157 enddo 156 zlon( iim+1)=zlon(1)+2.*RPI158 zlon(nbp_lon+1)=zlon(1)+2.*RPI 157 159 158 160 c verif init lat de 90 à -90, lon de -180 à 180 … … 256 258 257 259 use dimphy 260 use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat 258 261 implicit none 259 262 … … 283 286 c====================================================================== 284 287 285 #include "dimensions.h"286 288 #include "YOMCST.h" 287 289 c … … 289 291 c 290 292 real map_u(klon),map_v(klon),map_a(klon) 291 real latit( jjm+1),longit(iim)293 real latit(nbp_lat),longit(nbp_lon) 292 294 real phi,lam,ubal,vbal,abal 293 295 c … … 296 298 INTEGER i,j,k 297 299 INTEGER jj,ii 298 REAL ujj( iim+1),vjj(iim+1),ajj(iim+1)300 REAL ujj(nbp_lon+1),vjj(nbp_lon+1),ajj(nbp_lon+1) 299 301 REAL factlat,factlon 300 302 … … 302 304 c------------------------------------------------- 303 305 jj=1 ! POLE NORD 304 do j=2, jjm306 do j=2,nbp_lat-1 305 307 if (phi.lt.latit(j)) jj=j 306 308 enddo … … 309 311 c pole nord 310 312 if (jj.eq.1) then 311 do i=1, iim313 do i=1,nbp_lon 312 314 ujj(i) = map_u(i+1)*factlat + map_u(1)*(1-factlat) 313 315 vjj(i) = map_v(i+1)*factlat + map_v(1)*(1-factlat) … … 315 317 enddo 316 318 c pole sud 317 elseif (jj.eq. jjm) then318 do i=1, iim319 k = (jj-2)* iim+1+i319 elseif (jj.eq.nbp_lat-1) then 320 do i=1,nbp_lon 321 k = (jj-2)*nbp_lon+1+i 320 322 ujj(i) = map_u(klon)*factlat + map_u(k)*(1-factlat) 321 323 vjj(i) = map_v(klon)*factlat + map_v(k)*(1-factlat) … … 324 326 c autres latitudes 325 327 else 326 do i=1, iim327 k = (jj-2)* iim+1+i328 ujj(i) = map_u(k+ iim)*factlat + map_u(k)*(1-factlat)329 vjj(i) = map_v(k+ iim)*factlat + map_v(k)*(1-factlat)330 ajj(i) = map_a(k+ iim)*factlat + map_a(k)*(1-factlat)328 do i=1,nbp_lon 329 k = (jj-2)*nbp_lon+1+i 330 ujj(i) = map_u(k+nbp_lon)*factlat + map_u(k)*(1-factlat) 331 vjj(i) = map_v(k+nbp_lon)*factlat + map_v(k)*(1-factlat) 332 ajj(i) = map_a(k+nbp_lon)*factlat + map_a(k)*(1-factlat) 331 333 enddo 332 334 endif 333 ujj( iim+1)=ujj(1)334 vjj( iim+1)=vjj(1)335 ajj( iim+1)=ajj(1)335 ujj(nbp_lon+1)=ujj(1) 336 vjj(nbp_lon+1)=vjj(1) 337 ajj(nbp_lon+1)=ajj(1) 336 338 337 339 c Interpolation in longitudes 338 340 c------------------------------------------------- 339 341 ii=1 ! lon=-180 340 do i=2, iim342 do i=2,nbp_lon 341 343 if (lam.gt.longit(i)) ii=i 342 344 enddo … … 370 372 c====================================================================== 371 373 372 #include "dimensions.h"373 374 #include "YOMCST.h" 374 375 c -
trunk/LMDZ.TITAN/libf/phytitan/calchim.F
r1379 r1530 19 19 USE comgeomphy, only: rlatd 20 20 USE moyzon_mod, only: tmoy,playmoy,zlaymoy,zlevmoy,klat 21 use mod_grid_phy_lmdz, only: nbp_lat 21 22 implicit none 22 #include "dimensions.h"23 23 #include "clesphys.h" 24 #include "paramet.h"25 24 #include "YOMCST.h" 26 25 … … 105 104 c ************************************ 106 105 107 allocate(krpd(15,ND+1,NLRT, jjp1),krate(NLEV,NR),md(NLEV,NC))106 allocate(krpd(15,ND+1,NLRT,nbp_lat),krate(NLEV,NR),md(NLEV,NC)) 108 107 109 108 c Verification du nombre de composes: coherence common_mod et nqmax-nmicro … … 265 264 c taux de photodissociations: 266 265 c -------------------------- 267 call disso(krpd, jjp1)266 call disso(krpd,nbp_lat) 268 267 269 268 c reactions chimiques: -
trunk/LMDZ.TITAN/libf/phytitan/clmain.F
r1443 r1530 35 35 USE interface_surf 36 36 use dimphy 37 use mod_grid_phy_lmdz, only: nbp_lev 37 38 use cpdet_mod, only: t2tpot 38 39 IMPLICIT none … … 66 67 cAA la premiere couche 67 68 c====================================================================== 68 #include "dimensions.h"69 69 c$$$ PB ajout pour soil 70 70 #include "dimsoil.h" … … 324 324 325 325 ! ADAPTATION GCM POUR CP(T) 326 call t2tpot(knon* llm,yt,yteta,ppk)326 call t2tpot(knon*nbp_lev,yt,yteta,ppk) 327 327 328 328 yzlev(1:knon,1)=0. … … 479 479 USE interface_surf 480 480 use dimphy 481 use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat, nbp_lev 481 482 use cpdet_mod, only: t2tpot,tpot2t,cpdet 482 483 … … 486 487 c Objet: diffusion verticale de "h" 487 488 c====================================================================== 488 #include "dimensions.h"489 489 #include "YOMCST.h" 490 490 #include "dimsoil.h" … … 594 594 595 595 c passage en enthalpie potentielle 596 call t2tpot(knon* llm,t,local_h,ppk)596 call t2tpot(knon*nbp_lev,t,local_h,ppk) 597 597 c print*,"tpot en entree de clqh=",local_h(klon/2,:) 598 598 … … 685 685 ! ADAPTATION GCM POUR CP(T) 686 686 CALL interfsurf_hq(itime, dtime, rmu0, 687 e klon, iim, jjm, knon,687 e klon, nbp_lon, nbp_lat-1, knon, 688 688 e rlon, rlat, cufi, cvfi, 689 689 e debut, lafin, soil_model, nsoilmx,tsoil, … … 733 733 ENDDO 734 734 ENDDO 735 call tpot2t(knon* llm,local_h,d_t,ppk)735 call tpot2t(knon*nbp_lev,local_h,d_t,ppk) 736 736 737 737 c print*,"tpot en sortie de clqh=",local_h(klon/2,:) … … 781 781 c flux_v---output-R- (diagnostic) flux du vent: (kg m/s)/(m**2 s) 782 782 c====================================================================== 783 #include "dimensions.h"784 783 #include "iniprint.h" 785 784 INTEGER knon … … 908 907 c pcfh-----output-R- coefficients a calculer (chaleur et humidite) 909 908 c====================================================================== 910 #include "dimensions.h"911 909 #include "YOMCST.h" 912 910 #include "iniprint.h" … … 1172 1170 1173 1171 use dimphy 1172 use mod_grid_phy_lmdz, only: nbp_lev 1174 1173 use cpdet_mod, only: cpdet 1175 1174 IMPLICIT none … … 1188 1187 c pcfh-----output-R- coefficients a calculer (chaleur et humidite) 1189 1188 c====================================================================== 1190 #include "dimensions.h"1191 1189 #include "YOMCST.h" 1192 1190 #include "iniprint.h" -
trunk/LMDZ.TITAN/libf/phytitan/cltrac.F
r102 r1530 27 27 c flux_tr--output-R- flux de tr 28 28 c====================================================================== 29 #include "dimensions.h"30 29 REAL dtime 31 30 REAL coef(klon,klev) -
trunk/LMDZ.TITAN/libf/phytitan/coefkzmin.F
r1056 r1530 9 9 IMPLICIT NONE 10 10 11 #include "dimensions.h"12 11 #include "YOMCST.h" 13 12 -
trunk/LMDZ.TITAN/libf/phytitan/common_mod.F90
r1525 r1530 90 90 use dimphy 91 91 IMPLICIT NONE 92 #include "dimensions.h"93 92 #include "microtab.h" 94 93 -
trunk/LMDZ.TITAN/libf/phytitan/diagphy.F
r1048 r1530 51 51 implicit none 52 52 53 #include "dimensions.h"54 53 #include "YOMCST.h" 55 54 C … … 212 211 IMPLICIT NONE 213 212 C 214 #include "dimensions.h"215 213 #include "YOMCST.h" 216 214 C -
trunk/LMDZ.TITAN/libf/phytitan/drag_noro.F
r1356 r1530 11 11 use dimphy 12 12 IMPLICIT none 13 14 #include "dimensions.h"15 #include "paramet.h"16 13 17 14 c====================================================================== … … 30 27 c paprs---input-R-Pressure in semi layers (Pa) 31 28 c pplay---input-R-Pressure model-layers (Pa) 32 c pgeop---input-R-Geopotential model layers (reference ground)29 c pgeop---input-R-Geopotential model layers (reference to ground) 33 30 c pn2-----input-R-Brunt-Vaisala freq.^2 at 1/2 layers 34 31 c t-------input-R-temperature (K) -
trunk/LMDZ.TITAN/libf/phytitan/grid_noro.F
r1356 r1530 52 52 C======================================================================= 53 53 54 use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat 54 55 IMPLICIT none 55 56 56 #include "dimensions.h"57 57 #include "YOMCST.h" 58 58 … … 79 79 C INTERMEDIATE FIELDS (CORRELATIONS OF OROGRAPHY GRADIENT) 80 80 81 REAL ztz( iim+1,jjm+1),zxtzx(iim+1,jjm+1)82 REAL zytzy( iim+1,jjm+1),zxtzy(iim+1,jjm+1)83 REAL weight( iim+1,jjm+1)81 REAL ztz(nbp_lon+1,nbp_lat),zxtzx(nbp_lon+1,nbp_lat) 82 REAL zytzy(nbp_lon+1,nbp_lat),zxtzy(nbp_lon+1,nbp_lat) 83 REAL weight(nbp_lon+1,nbp_lat) 84 84 85 85 C CORRELATIONS OF USN OROGRAPHY GRADIENTS … … 107 107 c 108 108 c 109 if( iim.ne.imar) STOP 'Problem dim. x'110 if( jjm.ne.jmar-1) STOP 'Problem dim. y'109 if(nbp_lon.ne.imar) STOP 'Problem dim. x' 110 if(nbp_lat-1.ne.jmar-1) STOP 'Problem dim. y' 111 111 IF (imar.GT.2200 .OR. jmar.GT.1100) THEN 112 112 PRINT*, 'imar or jmar too big', imar, jmar … … 114 114 ENDIF 115 115 116 IF(imar+1.ne. iim+1.or.jmar.ne.jjm+1)THEN116 IF(imar+1.ne.nbp_lon+1.or.jmar.ne.nbp_lat)THEN 117 117 print *,' imar or jmar bad dimensions:',imar,jmar 118 118 call abort … … 297 297 C FIRST FILTER, MOVING AVERAGE OVER 9 POINTS. 298 298 299 CALL MVA9(zmea, iim+1,jjm+1)300 CALL MVA9(zstd, iim+1,jjm+1)301 CALL MVA9(zpic, iim+1,jjm+1)302 CALL MVA9(zval, iim+1,jjm+1)303 CALL MVA9(zxtzx, iim+1,jjm+1)304 CALL MVA9(zxtzy, iim+1,jjm+1)305 CALL MVA9(zytzy, iim+1,jjm+1)299 CALL MVA9(zmea,nbp_lon+1,nbp_lat) 300 CALL MVA9(zstd,nbp_lon+1,nbp_lat) 301 CALL MVA9(zpic,nbp_lon+1,nbp_lat) 302 CALL MVA9(zval,nbp_lon+1,nbp_lat) 303 CALL MVA9(zxtzx,nbp_lon+1,nbp_lat) 304 CALL MVA9(zxtzy,nbp_lon+1,nbp_lat) 305 CALL MVA9(zytzy,nbp_lon+1,nbp_lat) 306 306 307 307 DO ii = 1, imar -
trunk/LMDZ.TITAN/libf/phytitan/gwprofil.F
r1056 r1530 38 38 use dimphy 39 39 IMPLICIT NONE 40 41 #include "dimensions.h"42 #include "paramet.h"43 40 44 41 #include "YOMCST.h" -
trunk/LMDZ.TITAN/libf/phytitan/gwstress.F
r1056 r1530 54 54 use dimphy 55 55 implicit none 56 57 #include "dimensions.h"58 #include "paramet.h"59 56 60 57 #include "YOMCST.h" -
trunk/LMDZ.TITAN/libf/phytitan/hgardfou.F
r102 r1530 3 3 ! 4 4 SUBROUTINE hgardfou (t,tsol,text) 5 6 use dimphy 7 IMPLICIT none 5 8 c====================================================================== 6 9 c Verifier la temperature 7 10 c====================================================================== 8 use dimphy9 IMPLICIT none10 #include "dimensions.h"11 11 #include "YOMCST.h" 12 12 REAL t(klon,klev), tsol(klon) … … 84 84 ok = .FALSE. 85 85 DO i = 1, jbad 86 PRINT *,'i,temperature sol=',jadrs(i),zt(jadrs(i))86 PRINT *,'i,temperature =',jadrs(i),zt(jadrs(i)) 87 87 ENDDO 88 88 ENDIF … … 101 101 ok = .FALSE. 102 102 DO i = 1, jbad 103 PRINT *,'i,temperature sol=',jadrs(i),zt(jadrs(i))103 PRINT *,'i,temperature =',jadrs(i),zt(jadrs(i)) 104 104 ENDDO 105 105 ENDIF -
trunk/LMDZ.TITAN/libf/phytitan/ini_histday.h
r1056 r1530 22 22 23 23 CALL histdef(nid_day, "phis", "Surface geop. height", "-", 24 . iim,jj_nb,nhori, 1,1,1, nvert, 32,24 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 25 25 . "once", zsto,zout) 26 26 27 27 CALL histdef(nid_day, "aire", "Grid area", "-", 28 . iim,jj_nb,nhori, 1,1,1, nvert, 32,28 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 29 29 . "once", zsto,zout) 30 30 31 31 ccccccc axe Ls 32 32 CALL histdef(nid_day, "ls", "Solar longitude", "degrees", 33 . iim,jj_nb,nhori, 1,1,1, nvert, 32,33 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 34 34 . "ave(X)", zsto,zout) 35 35 … … 37 37 38 38 CALL histdef(nid_day, "tsol", "Surface Temperature", "K", 39 . iim,jj_nb,nhori, 1,1,1, nvert, 32,39 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 40 40 . "ave(X)", zsto,zout) 41 41 42 42 CALL histdef(nid_day, "psol", "Surface Pressure", "Pa", 43 . iim,jj_nb,nhori, 1,1,1, nvert, 32,43 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 44 44 . "ave(X)", zsto,zout) 45 45 46 46 c CALL histdef(nid_day, "ue", "Zonal energy transport", "-", 47 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,47 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 48 48 c . "ave(X)", zsto,zout) 49 49 50 50 c CALL histdef(nid_day, "ve", "Merid energy transport", "-", 51 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,51 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 52 52 c . "ave(X)", zsto,zout) 53 53 … … 60 60 61 61 CALL histdef(nid_day, "temp", "Air temperature", "K", 62 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,62 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 63 63 . "ave(X)", zsto,zout) 64 64 65 65 CALL histdef(nid_day, "pres", "Air pressure", "Pa", 66 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,66 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 67 67 . "ave(X)", zsto,zout) 68 68 69 69 CALL histdef(nid_day, "geop", "Geopotential height", "m", 70 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,70 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 71 71 . "ave(X)", zsto,zout) 72 72 73 73 CALL histdef(nid_day, "vitu", "Zonal wind", "m/s", 74 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,74 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 75 75 . "ave(X)", zsto,zout) 76 76 77 77 CALL histdef(nid_day, "vitv", "Meridional wind", "m/s", 78 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,78 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 79 79 . "ave(X)", zsto,zout) 80 80 81 81 CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s", 82 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,82 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 83 83 . "ave(X)", zsto,zout) 84 84 85 85 CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2", 86 . iim,jj_nb,nhori, 1,1,1, nvert, 32,86 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 87 87 . "ave(X)", zsto1,zout) 88 88 89 89 CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2", 90 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,90 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 91 91 . "ave(X)", zsto,zout) 92 92 93 93 CALL histdef(nid_day, "dudyn", "Dynamics dU", "m/s2", 94 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,94 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 95 95 . "ave(X)", zsto,zout) 96 96 … … 101 101 c DO iq=1,nmicro 102 102 c CALL histdef(nid_day, tname(iq), ttext(iq), "n/m2", 103 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,103 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 104 104 c . "ave(X)", zsto,zout) 105 105 c ENDDO 106 106 CALL histdef(nid_day, "qaer","nb tot aer" , "n/m2", 107 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,107 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 108 108 . "ave(X)", zsto,zout) 109 109 110 110 if (clouds.eq.1) then 111 111 CALL histdef(nid_day, "qnoy","nb tot noy" , "n/m2", 112 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,112 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 113 113 . "ave(X)", zsto,zout) 114 114 CALL histdef(nid_day, "qgl1","V tot gl1" , "m3/m2", 115 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,115 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 116 116 . "ave(X)", zsto,zout) 117 117 CALL histdef(nid_day, "qgl2","V tot gl2" , "m3/m2", 118 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,118 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 119 119 . "ave(X)", zsto,zout) 120 120 CALL histdef(nid_day, "qgl3","V tot gl3" , "m3/m2", 121 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,121 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 122 122 . "ave(X)", zsto,zout) 123 123 c-------------- 124 124 c ----- SATURATION ESP NUAGES 125 125 CALL histdef(nid_day,"ch4sat", "saturation CH4", "--", 126 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,126 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 127 127 . "ave(X)", zsto,zout) 128 128 CALL histdef(nid_day,"c2h6sat", "saturation C2H6", "--", 129 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,129 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 130 130 . "ave(X)", zsto,zout) 131 131 CALL histdef(nid_day,"c2h2sat", "saturation C2H2", "--", 132 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,132 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 133 133 . "ave(X)", zsto,zout) 134 134 c -------------- 135 135 c ----- RESERVOIR DE SURFACE 136 136 CALL histdef(nid_day, "reserv", "Reservoir surface","m", 137 . iim,jj_nb,nhori, 1,1,1, nvert, 32,137 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 138 138 . "ave(X)", zsto,zout) 139 139 c -------------- 140 140 c ----- ECHANGE GAZ SURF/ATM (evaporation) 141 141 CALL histdef(nid_day, "evapch4", "Evaporation CH4","m", 142 . iim,jj_nb,nhori, 1,1,1, nvert, 32,142 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 143 143 . "ave(X)", zsto,zout) 144 144 c -------------- 145 145 c ----- PRECIPITATIONS (precipitations cumulatives) 146 146 CALL histdef(nid_day,"prech4","Precip CH4","m", 147 . iim,jj_nb,nhori, 1,1,1, nvert, 32,147 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 148 148 . "ave(X)", zsto,zout) 149 149 CALL histdef(nid_day,"prec2h6","Precip C2H6", 150 . "m", iim,jj_nb,nhori, 1,1,1, nvert, 32,150 . "m",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 151 151 . "ave(X)", zsto,zout) 152 152 CALL histdef(nid_day,"prec2h2","Precip C2H2", 153 . "m", iim,jj_nb,nhori, 1,1,1, nvert, 32,153 . "m",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 154 154 . "ave(X)", zsto,zout) 155 155 CALL histdef(nid_day,"prenoy","Precip NOY", 156 . "um/s", iim,jj_nb,nhori, 1,1,1, nvert, 32,156 . "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 157 157 . "ave(X)", zsto,zout) 158 158 CALL histdef(nid_day,"preaer","Precip AER", 159 . "um/s", iim,jj_nb,nhori, 1,1,1, nvert, 32,159 . "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 160 160 . "ave(X)", zsto,zout) 161 161 c -------------- 162 162 c ----- FLUX GLACE 163 163 CALL histdef(nid_day,"flxgl1", "flux gl CH4", 164 . "kg/m2/s", iim,jj_nb,nhori, klev,1,klev,nvert,32,164 . "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 165 165 . "ave(X)", zsto,zout) 166 166 CALL histdef(nid_day,"flxgl2", "flux gl C2H6", 167 . "kg/m2/s", iim,jj_nb,nhori, klev,1,klev,nvert,32,167 . "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 168 168 . "ave(X)", zsto,zout) 169 169 CALL histdef(nid_day,"flxgl3", "flux gl C2H2", 170 . "kg/m2/s", iim,jj_nb,nhori, klev,1,klev,nvert,32,170 . "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 171 171 . "ave(X)", zsto,zout) 172 172 c -------------- 173 173 c ----- RAYON DES GOUTTES 174 174 CALL histdef(nid_day,"rcldbar", "rayon moyen goutte", 175 . "m", iim,jj_nb,nhori, klev,1,klev,nvert, 32,175 . "m",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 176 176 . "ave(X)", zsto,zout) 177 177 endif … … 182 182 DO iq=nmicro+1,nqmax 183 183 CALL histdef(nid_day, tname(iq), ttext(iq), "ppm", 184 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,184 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 185 185 . "ave(X)", zsto,zout) 186 186 ENDDO … … 198 198 199 199 CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2", 200 . iim,jj_nb,nhori, 1,1,1, nvert, 32,200 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 201 201 . "ave(X)", zsto1,zout) 202 202 203 203 CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2", 204 . iim,jj_nb,nhori, 1,1,1, nvert, 32,204 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 205 205 . "ave(X)", zsto1,zout) 206 206 207 207 CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2", 208 . iim,jj_nb,nhori, 1,1,1, nvert, 32,208 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 209 209 . "ave(X)", zsto1,zout) 210 210 … … 212 212 213 213 CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2", 214 . iim,jj_nb,nhori, klev,1,klev,nvert,214 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 215 215 . 32, "ave(X)", zsto1,zout) 216 216 217 217 CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2", 218 . iim,jj_nb,nhori, klev,1,klev,nvert,218 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 219 219 . 32, "ave(X)", zsto1,zout) 220 220 … … 224 224 write(str2,'(i2.2)') k 225 225 CALL histdef(nid_day,"thv"//str2,"Haze Opa Vis", 226 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,226 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 227 227 . "ave(X)",zsto1,zout) 228 228 ENDDO … … 231 231 write(str2,'(i2.2)') k 232 232 CALL histdef(nid_day,"thi"//str2,"Haze Opa IR", 233 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,233 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 234 234 . "ave(X)",zsto1,zout) 235 235 ENDDO … … 240 240 write(str2,'(i2.2)') k 241 241 CALL histdef(nid_day,"khv"//str2,"Haze ext Vis ", 242 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,242 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 243 243 . "ave(X)",zsto1,zout) 244 244 ENDDO … … 247 247 write(str2,'(i2.2)') k 248 248 CALL histdef(nid_day,"khi"//str2,"Haze ext IR ", 249 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,249 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 250 250 . "ave(X)",zsto1,zout) 251 251 ENDDO … … 256 256 write(str2,'(i2.2)') k 257 257 CALL histdef(nid_day,"tgv"//str2,"Gas Opa Vis", 258 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,258 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 259 259 . "ave(X)",zsto1,zout) 260 260 ENDDO … … 263 263 write(str2,'(i2.2)') k 264 264 CALL histdef(nid_day,"tgi"//str2,"Gas Opa IR", 265 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,265 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 266 266 . "ave(X)",zsto1,zout) 267 267 ENDDO … … 272 272 write(str2,'(i2.2)') k 273 273 CALL histdef(nid_day,"kgv"//str2,"Gas ext Vis ", 274 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,274 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 275 275 . "ave(X)",zsto1,zout) 276 276 ENDDO … … 279 279 write(str2,'(i2.2)') k 280 280 CALL histdef(nid_day,"kgi"//str2,"Gas ext IR ", 281 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,281 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 282 282 . "ave(X)",zsto1,zout) 283 283 ENDDO … … 287 287 if (clouds.eq.1) then 288 288 CALL histdef(nid_day,"tcld","Cld Opa proxy", 289 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,289 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 290 290 . "ave(X)",zsto,zout) 291 291 … … 293 293 c ----- EXTINCTION NUAGES 294 294 CALL histdef(nid_day,"kcld","Cld Ext proxy", 295 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,295 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 296 296 . "ave(X)",zsto,zout) 297 297 endif … … 303 303 304 304 CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s", 305 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,305 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 306 306 . "ave(X)", zsto,zout) 307 307 308 308 CALL histdef(nid_day, "dtphy", "Physics dT", "K/s", 309 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,309 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 310 310 . "ave(X)", zsto,zout) 311 311 312 312 CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s", 313 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,313 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 314 314 . "ave(X)", zsto,zout) 315 315 316 316 CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s", 317 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,317 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 318 318 . "ave(X)", zsto,zout) 319 319 320 320 CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s", 321 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,321 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 322 322 . "ave(X)", zsto,zout) 323 323 324 324 CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s", 325 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,325 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 326 326 . "ave(X)", zsto,zout) 327 327 328 328 c CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s", 329 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,329 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 330 330 c . "ave(X)", zsto,zout) 331 331 … … 338 338 c call histdef(nid_day, "taux", 339 339 c $ "Zonal wind stress", "Pa", 340 c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,340 c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 341 341 c $ "ave(X)", zsto,zout) 342 342 343 343 c call histdef(nid_day, "tauy", 344 344 c $ "Meridional xind stress", "Pa", 345 c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,345 c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 346 346 c $ "ave(X)", zsto,zout) 347 347 348 348 c CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-", 349 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,349 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 350 350 c . "ave(X)", zsto,zout) 351 351 352 352 c CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-", 353 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,353 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 354 354 c . "ave(X)", zsto,zout) 355 355 -
trunk/LMDZ.TITAN/libf/phytitan/ini_histins.h
r1056 r1530 21 21 22 22 CALL histdef(nid_ins, "phis", "Surface geop. height", "-", 23 . iim,jj_nb,nhori, 1,1,1, nvert, 32,23 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 24 24 . "once", zsto,zout) 25 25 26 26 CALL histdef(nid_ins, "aire", "Grid area", "-", 27 . iim,jj_nb,nhori, 1,1,1, nvert, 32,27 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 28 28 . "once", zsto,zout) 29 29 30 30 ccccccc axe Ls 31 31 CALL histdef(nid_ins, "ls", "Solar longitude", "degrees", 32 . iim,jj_nb,nhori, 1,1,1, nvert, 32,32 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 33 33 . "inst(X)", zsto,zout) 34 34 … … 36 36 37 37 CALL histdef(nid_ins, "tsol", "Surface Temperature", "K", 38 . iim,jj_nb,nhori, 1,1,1, nvert, 32,38 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 39 39 . "inst(X)", zsto,zout) 40 40 41 41 CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa", 42 . iim,jj_nb,nhori, 1,1,1, nvert, 32,42 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 43 43 . "inst(X)", zsto,zout) 44 44 45 45 c CALL histdef(nid_ins, "ue", "Zonal energy transport", "-", 46 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,46 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 47 47 c . "inst(X)", zsto,zout) 48 48 49 49 c CALL histdef(nid_ins, "ve", "Merid energy transport", "-", 50 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,50 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 51 51 c . "inst(X)", zsto,zout) 52 52 … … 59 59 60 60 CALL histdef(nid_ins, "temp", "Air temperature", "K", 61 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,61 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 62 62 . "inst(X)", zsto,zout) 63 63 64 64 CALL histdef(nid_ins, "pres", "Air pressure", "Pa", 65 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,65 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 66 66 . "inst(X)", zsto,zout) 67 67 68 68 CALL histdef(nid_ins, "geop", "Geopotential height", "m", 69 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,69 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 70 70 . "inst(X)", zsto,zout) 71 71 72 72 CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s", 73 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,73 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 74 74 . "inst(X)", zsto,zout) 75 75 76 76 CALL histdef(nid_ins, "vitv", "Meridional wind", "m/s", 77 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,77 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 78 78 . "inst(X)", zsto,zout) 79 79 80 80 CALL histdef(nid_ins, "vitw", "Vertical wind", "Pa/s", 81 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,81 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 82 82 . "inst(X)", zsto,zout) 83 83 84 84 CALL histdef(nid_ins, "tops", "Solar rad. at TOA", "W/m2", 85 . iim,jj_nb,nhori, 1,1,1, nvert, 32,85 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 86 86 . "inst(X)", zsto,zout) 87 87 88 88 c CALL histdef(nid_ins, "duvdf", "Boundary-layer dU", "m/s2", 89 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,89 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 90 90 c . "inst(X)", zsto,zout) 91 91 92 92 c CALL histdef(nid_ins, "dudyn", "Dynamics dU", "m/s2", 93 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,93 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 94 94 c . "inst(X)", zsto,zout) 95 95 … … 105 105 DO iq=1,nmicro 106 106 CALL histdef(nid_ins, tname(iq), ttext(iq), "n/m2", 107 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,107 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 108 108 . "inst(X)", zsto,zout) 109 109 ENDDO … … 112 112 DO iq=nmicro+1,nqmax 113 113 CALL histdef(nid_ins, tname(iq), ttext(iq), "ppm", 114 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,114 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 115 115 . "inst(X)", zsto,zout) 116 116 ENDDO … … 123 123 124 124 CALL histdef(nid_ins, "topl", "IR rad. at TOA", "W/m2", 125 . iim,jj_nb,nhori, 1,1,1, nvert, 32,125 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 126 126 . "inst(X)", zsto,zout) 127 127 128 128 CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2", 129 . iim,jj_nb,nhori, 1,1,1, nvert, 32,129 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 130 130 . "inst(X)", zsto,zout) 131 131 132 132 CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2", 133 . iim,jj_nb,nhori, 1,1,1, nvert, 32,133 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 134 134 . "inst(X)", zsto,zout) 135 135 … … 137 137 138 138 CALL histdef(nid_ins, "SWnet", "Net SW flux","W/m2", 139 . iim,jj_nb,nhori, klev,1,klev,nvert,139 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 140 140 . 32, "inst(X)", zsto,zout) 141 141 142 142 CALL histdef(nid_ins, "LWnet", "Net LW flux","W/m2", 143 . iim,jj_nb,nhori, klev,1,klev,nvert,143 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 144 144 . 32, "inst(X)", zsto,zout) 145 145 … … 149 149 write(str2,'(i2.2)') k 150 150 CALL histdef(nid_ins,"thv"//str2,"Haze Opa Vis", 151 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,151 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 152 152 . "ins(X)",zsto,zout) 153 153 ENDDO … … 156 156 write(str2,'(i2.2)') k 157 157 CALL histdef(nid_ins,"thi"//str2,"Haze Opa IR", 158 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,158 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 159 159 . "ins(X)",zsto,zout) 160 160 ENDDO … … 165 165 write(str2,'(i2.2)') k 166 166 CALL histdef(nid_ins,"khv"//str2,"Haze ext Vis ", 167 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,167 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 168 168 . "ins(X)",zsto,zout) 169 169 ENDDO … … 172 172 write(str2,'(i2.2)') k 173 173 CALL histdef(nid_ins,"khi"//str2,"Haze ext IR ", 174 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,174 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 175 175 . "ins(X)",zsto,zout) 176 176 ENDDO … … 181 181 write(str2,'(i2.2)') k 182 182 CALL histdef(nid_ins,"tgv"//str2,"Haze Opa Vis", 183 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,183 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 184 184 . "ins(X)",zsto,zout) 185 185 ENDDO … … 188 188 write(str2,'(i2.2)') k 189 189 CALL histdef(nid_ins,"tgi"//str2,"Haze Opa IR", 190 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,190 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 191 191 . "ins(X)",zsto,zout) 192 192 ENDDO … … 197 197 write(str2,'(i2.2)') k 198 198 CALL histdef(nid_ins,"kgv"//str2,"Haze ext Vis ", 199 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,199 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 200 200 . "ins(X)",zsto,zout) 201 201 ENDDO … … 204 204 write(str2,'(i2.2)') k 205 205 CALL histdef(nid_ins,"kgi"//str2,"Haze ext IR ", 206 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,206 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 207 207 . "ins(X)",zsto,zout) 208 208 ENDDO … … 214 214 215 215 CALL histdef(nid_ins, "dtdyn", "Dynamics dT", "K/s", 216 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,216 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 217 217 . "inst(X)", zsto,zout) 218 218 219 219 CALL histdef(nid_ins, "dtphy", "Physics dT", "K/s", 220 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,220 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 221 221 . "inst(X)", zsto,zout) 222 222 223 223 CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s", 224 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,224 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 225 225 . "inst(X)", zsto,zout) 226 226 227 227 CALL histdef(nid_ins, "dtajs", "Dry adjust. dT", "K/s", 228 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,228 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 229 229 . "inst(X)", zsto,zout) 230 230 231 231 CALL histdef(nid_ins, "dtswr", "SW radiation dT", "K/s", 232 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,232 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 233 233 . "inst(X)", zsto,zout) 234 234 235 235 CALL histdef(nid_ins, "dtlwr", "LW radiation dT", "K/s", 236 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,236 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 237 237 . "inst(X)", zsto,zout) 238 238 239 239 c CALL histdef(nid_ins, "dtec", "Cinetic dissip dT", "K/s", 240 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,240 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 241 241 c . "inst(X)", zsto,zout) 242 242 243 243 c CALL histdef(nid_ins, "dvvdf", "Boundary-layer dV", "m/s2", 244 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,244 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 245 245 c . "inst(X)", zsto,zout) 246 246 … … 253 253 c call histdef(nid_ins, "taux", 254 254 c $ "Zonal wind stress", "Pa", 255 c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,255 c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 256 256 c $ "inst(X)", zsto,zout) 257 257 258 258 c call histdef(nid_ins, "tauy", 259 259 c $ "Meridional xind stress", "Pa", 260 c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,260 c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 261 261 c $ "inst(X)", zsto,zout) 262 262 263 263 c CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-", 264 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,264 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 265 265 c . "inst(X)", zsto,zout) 266 266 267 267 c CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-", 268 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,268 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 269 269 c . "inst(X)", zsto,zout) 270 270 -
trunk/LMDZ.TITAN/libf/phytitan/ini_histmth.h
r1356 r1530 22 22 23 23 CALL histdef(nid_mth, "phis", "Surface geop. height", "-", 24 . iim,jj_nb,nhori, 1,1,1, nvert, 32,24 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 25 25 . "once", zsto,zout) 26 26 27 27 CALL histdef(nid_mth, "aire", "Grid area", "-", 28 . iim,jj_nb,nhori, 1,1,1, nvert, 32,28 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 29 29 . "once", zsto,zout) 30 30 31 31 ccccccc axe Ls 32 32 CALL histdef(nid_mth, "ls", "Solar longitude", "degrees", 33 . iim,jj_nb,nhori, 1,1,1, nvert, 32,33 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 34 34 . "ave(X)", zsto,zout) 35 35 … … 37 37 38 38 CALL histdef(nid_mth, "tsol", "Surface Temperature", "K", 39 . iim,jj_nb,nhori, 1,1,1, nvert, 32,39 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 40 40 . "ave(X)", zsto,zout) 41 41 42 42 CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa", 43 . iim,jj_nb,nhori, 1,1,1, nvert, 32,43 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 44 44 . "ave(X)", zsto,zout) 45 45 46 46 c CALL histdef(nid_mth, "ue", "Zonal energy transport", "-", 47 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,47 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 48 48 c . "ave(X)", zsto,zout) 49 49 50 50 c CALL histdef(nid_mth, "ve", "Merid energy transport", "-", 51 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,51 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 52 52 c . "ave(X)", zsto,zout) 53 53 … … 60 60 61 61 CALL histdef(nid_mth, "temp", "Air temperature", "K", 62 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,62 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 63 63 . "ave(X)", zsto,zout) 64 64 65 65 CALL histdef(nid_mth, "pres", "Air pressure", "Pa", 66 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,66 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 67 67 . "ave(X)", zsto,zout) 68 68 69 69 CALL histdef(nid_mth, "geop", "Geopotential height", "m", 70 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,70 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 71 71 . "ave(X)", zsto,zout) 72 72 73 73 CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s", 74 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,74 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 75 75 . "ave(X)", zsto,zout) 76 76 77 77 CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s", 78 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,78 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 79 79 . "ave(X)", zsto,zout) 80 80 81 81 CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s", 82 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,82 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 83 83 . "ave(X)", zsto,zout) 84 84 85 85 c CALL histdef(nid_mth, "Kz", "vertical diffusion coef", "m2/s", 86 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,86 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 87 87 c . "ave(X)", zsto,zout) 88 88 89 89 CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2", 90 . iim,jj_nb,nhori, 1,1,1, nvert, 32,90 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 91 91 . "ave(X)", zsto1,zout) 92 92 93 93 CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2", 94 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,94 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 95 95 . "ave(X)", zsto,zout) 96 96 97 97 CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2", 98 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,98 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 99 99 . "ave(X)", zsto,zout) 100 100 … … 105 105 c DO iq=1,nmicro 106 106 c CALL histdef(nid_mth, tname(iq), ttext(iq), "n/m2", 107 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,107 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 108 108 c . "ave(X)", zsto,zout) 109 109 c ENDDO 110 110 CALL histdef(nid_mth, "qaer","nb tot aer" , "n/m2", 111 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,111 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 112 112 . "ave(X)", zsto,zout) 113 113 114 114 if (clouds.eq.1) then 115 115 CALL histdef(nid_mth, "qnoy","nb tot noy" , "n/m2", 116 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,116 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 117 117 . "ave(X)", zsto,zout) 118 118 CALL histdef(nid_mth, "qgl1","V tot gl1" , "m3/m2", 119 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,119 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 120 120 . "ave(X)", zsto,zout) 121 121 CALL histdef(nid_mth, "qgl2","V tot gl2" , "m3/m2", 122 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,122 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 123 123 . "ave(X)", zsto,zout) 124 124 CALL histdef(nid_mth, "qgl3","V tot gl3" , "m3/m2", 125 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,125 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 126 126 . "ave(X)", zsto,zout) 127 127 c-------------- 128 128 c ----- SATURATION ESP NUAGES 129 129 CALL histdef(nid_mth,"ch4sat", "saturation CH4", "--", 130 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,130 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 131 131 . "ave(X)", zsto,zout) 132 132 CALL histdef(nid_mth,"c2h6sat", "saturation C2H6", "--", 133 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,133 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 134 134 . "ave(X)", zsto,zout) 135 135 CALL histdef(nid_mth,"c2h2sat", "saturation C2H2", "--", 136 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,136 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 137 137 . "ave(X)", zsto,zout) 138 138 c -------------- 139 139 c ----- RESERVOIR DE SURFACE 140 140 CALL histdef(nid_mth, "reserv", "Reservoir surface","m", 141 . iim,jj_nb,nhori, 1,1,1, nvert, 32,141 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 142 142 . "ave(X)", zsto,zout) 143 143 c -------------- 144 144 c ----- ECHANGE GAZ SURF/ATM (evaporation) 145 145 CALL histdef(nid_mth, "evapch4", "Evaporation CH4","m", 146 . iim,jj_nb,nhori, 1,1,1, nvert, 32,146 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 147 147 . "ave(X)", zsto,zout) 148 148 c -------------- 149 149 c ----- PRECIPITATIONS (precipitations moyennes) 150 150 CALL histdef(nid_mth,"prech4","Precip CH4","um/s", 151 . iim,jj_nb,nhori, 1,1,1, nvert, 32,151 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 152 152 . "ave(X)", zsto,zout) 153 153 CALL histdef(nid_mth,"prec2h6","Precip C2H6", 154 . "um/s", iim,jj_nb,nhori, 1,1,1, nvert, 32,154 . "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 155 155 . "ave(X)", zsto,zout) 156 156 CALL histdef(nid_mth,"prec2h2","Precip C2H2", 157 . "um/s", iim,jj_nb,nhori, 1,1,1, nvert, 32,157 . "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 158 158 . "ave(X)", zsto,zout) 159 159 CALL histdef(nid_mth,"prenoy","Precip NOY", 160 . "um/s", iim,jj_nb,nhori, 1,1,1, nvert, 32,160 . "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 161 161 . "ave(X)", zsto,zout) 162 162 CALL histdef(nid_mth,"preaer","Precip AER", 163 . "um/s", iim,jj_nb,nhori, 1,1,1, nvert, 32,163 . "um/s",nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 164 164 . "ave(X)", zsto,zout) 165 165 c -------------- 166 166 c ----- FLUX GLACE 167 167 CALL histdef(nid_mth,"flxgl1", "flux gl CH4", 168 . "kg/m2/s", iim,jj_nb,nhori, klev,1,klev,nvert,32,168 . "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 169 169 . "ave(X)", zsto,zout) 170 170 CALL histdef(nid_mth,"flxgl2", "flux gl C2H6", 171 . "kg/m2/s", iim,jj_nb,nhori, klev,1,klev,nvert,32,171 . "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 172 172 . "ave(X)", zsto,zout) 173 173 CALL histdef(nid_mth,"flxgl3", "flux gl C2H2", 174 . "kg/m2/s", iim,jj_nb,nhori, klev,1,klev,nvert,32,174 . "kg/m2/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 175 175 . "ave(X)", zsto,zout) 176 176 c -------------- 177 177 c ----- Source/puits GLACE 178 178 CALL histdef(nid_mth,"solch4", "dQ gl CH4", 179 . "m3/m3", iim,jj_nb,nhori, klev,1,klev,nvert, 32,179 . "m3/m3",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 180 180 . "ave(X)", zsto,zout) 181 181 CALL histdef(nid_mth,"solc2h6", "dQ gl C2H6", 182 . "m3/m3", iim,jj_nb,nhori, klev,1,klev,nvert, 32,182 . "m3/m3",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 183 183 . "ave(X)", zsto,zout) 184 184 CALL histdef(nid_mth,"solc2h2", "dQ gl C2H2", 185 . "m3/m3", iim,jj_nb,nhori, klev,1,klev,nvert, 32,185 . "m3/m3",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 186 186 . "ave(X)", zsto,zout) 187 187 c -------------- 188 188 c ----- RAYON DES GOUTTES 189 189 CALL histdef(nid_mth,"rcldbar", "rayon moyen goutte", 190 . "m", iim,jj_nb,nhori, klev,1,klev,nvert, 32,190 . "m",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 191 191 . "ave(X)", zsto,zout) 192 192 endif … … 197 197 DO iq=nmicro+1,nqmax 198 198 CALL histdef(nid_mth, tname(iq), ttext(iq), "ppm", 199 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,199 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 200 200 . "ave(X)", zsto,zout) 201 201 ENDDO … … 203 203 c DO iq=nmicro+1,nqmax 204 204 c CALL histdef(nid_mth, "c_"//tname(iq), "c_"//ttext(iq), 205 c . "ppm/s", iim,jj_nb,nhori, klev,1,klev,nvert, 32,205 c . "ppm/s",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 206 206 c . "ave(X)", zsto,zout) 207 207 c ENDDO … … 219 219 220 220 CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2", 221 . iim,jj_nb,nhori, 1,1,1, nvert, 32,221 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 222 222 . "ave(X)", zsto1,zout) 223 223 224 224 CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2", 225 . iim,jj_nb,nhori, 1,1,1, nvert, 32,225 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 226 226 . "ave(X)", zsto1,zout) 227 227 228 228 CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2", 229 . iim,jj_nb,nhori, 1,1,1, nvert, 32,229 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 230 230 . "ave(X)", zsto1,zout) 231 231 … … 233 233 234 234 CALL histdef(nid_mth, "SWnet", "Net SW flux","W/m2", 235 . iim,jj_nb,nhori, klev,1,klev,nvert,235 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 236 236 . 32, "ave(X)", zsto1,zout) 237 237 238 238 c CALL histdef(nid_mth, "SWup", "upward SW flux","W/m2", 239 c . iim,jj_nb,nhori, klev,1,klev,nvert,239 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 240 240 c . 32, "ave(X)", zsto1,zout) 241 241 242 242 c CALL histdef(nid_mth, "SWdn", "downward SW flux","W/m2", 243 c . iim,jj_nb,nhori, klev,1,klev,nvert,243 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 244 244 c . 32, "ave(X)", zsto1,zout) 245 245 246 246 CALL histdef(nid_mth, "LWnet", "Net LW flux","W/m2", 247 . iim,jj_nb,nhori, klev,1,klev,nvert,247 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 248 248 . 32, "ave(X)", zsto1,zout) 249 249 250 250 c CALL histdef(nid_mth, "LWup", "upward LW flux","W/m2", 251 c . iim,jj_nb,nhori, klev,1,klev,nvert,251 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 252 252 c . 32, "ave(X)", zsto1,zout) 253 253 254 254 c CALL histdef(nid_mth, "LWdn", "downward LW flux","W/m2", 255 c . iim,jj_nb,nhori, klev,1,klev,nvert,255 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 256 256 c . 32, "ave(X)", zsto1,zout) 257 257 258 258 CALL histdef(nid_mth, "fluxvdf", "PBL net flux","W/m2", 259 . iim,jj_nb,nhori, klev,1,klev,nvert,259 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 260 260 . 32, "ave(X)", zsto,zout) 261 261 262 262 CALL histdef(nid_mth, "fluxdyn", "Dyn. net flux","W/m2", 263 . iim,jj_nb,nhori, klev,1,klev,nvert,263 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 264 264 . 32, "ave(X)", zsto,zout) 265 265 266 266 CALL histdef(nid_mth, "fluxajs", "Dry adj. net flux","W/m2", 267 . iim,jj_nb,nhori, klev,1,klev,nvert,267 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 268 268 . 32, "ave(X)", zsto,zout) 269 269 270 270 c CALL histdef(nid_mth, "fluxec", "Cin. net flux","W/m2", 271 c . iim,jj_nb,nhori, klev,1,klev,nvert,271 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 272 272 c . 32, "ave(X)", zsto,zout) 273 273 … … 277 277 write(str2,'(i2.2)') k 278 278 CALL histdef(nid_mth,"thv"//str2,"Haze Opa Vis", 279 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,279 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 280 280 . "ave(X)",zsto1,zout) 281 281 ENDDO … … 284 284 write(str2,'(i2.2)') k 285 285 CALL histdef(nid_mth,"thi"//str2,"Haze Opa IR", 286 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,286 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 287 287 . "ave(X)",zsto1,zout) 288 288 ENDDO … … 293 293 write(str2,'(i2.2)') k 294 294 CALL histdef(nid_mth,"khv"//str2,"Haze ext Vis ", 295 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,295 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 296 296 . "ave(X)",zsto1,zout) 297 297 ENDDO … … 300 300 write(str2,'(i2.2)') k 301 301 CALL histdef(nid_mth,"khi"//str2,"Haze ext IR ", 302 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,302 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 303 303 . "ave(X)",zsto1,zout) 304 304 ENDDO … … 309 309 write(str2,'(i2.2)') k 310 310 CALL histdef(nid_mth,"tgv"//str2,"Gas Opa Vis", 311 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,311 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 312 312 . "ave(X)",zsto1,zout) 313 313 ENDDO … … 316 316 write(str2,'(i2.2)') k 317 317 CALL histdef(nid_mth,"tgi"//str2,"Haze Opa IR", 318 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,318 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 319 319 . "ave(X)",zsto1,zout) 320 320 ENDDO … … 325 325 write(str2,'(i2.2)') k 326 326 CALL histdef(nid_mth,"kgv"//str2,"Gas ext Vis ", 327 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,327 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 328 328 . "ave(X)",zsto1,zout) 329 329 ENDDO … … 332 332 write(str2,'(i2.2)') k 333 333 CALL histdef(nid_mth,"kgi"//str2,"Gas ext IR ", 334 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,334 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 335 335 . "ave(X)",zsto1,zout) 336 336 ENDDO … … 340 340 if (clouds.eq.1) then 341 341 CALL histdef(nid_mth,"tcld","Cld Opa proxy", 342 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,342 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 343 343 . "ave(X)",zsto,zout) 344 344 … … 346 346 c ----- EXTINCTION NUAGES 347 347 CALL histdef(nid_mth,"kcld","Cld Ext proxy", 348 . "m-1", iim,jj_nb,nhori,klev,1,klev,nvert,32,348 . "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 349 349 . "ave(X)",zsto,zout) 350 350 endif … … 355 355 write(str2,'(i2.2)') k 356 356 CALL histdef(nid_mth,"occcld"//str2,"occ cld", 357 . "--", iim,jj_nb,nhori,klev,1,klev,nvert,32,357 . "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 358 358 . "ave(X)",zsto,zout) 359 359 enddo … … 365 365 366 366 CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s", 367 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,367 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 368 368 . "ave(X)", zsto,zout) 369 369 370 370 CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s", 371 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,371 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 372 372 . "ave(X)", zsto,zout) 373 373 374 374 CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s", 375 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,375 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 376 376 . "ave(X)", zsto,zout) 377 377 378 378 CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s", 379 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,379 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 380 380 . "ave(X)", zsto,zout) 381 381 382 382 CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s", 383 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,383 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 384 384 . "ave(X)", zsto,zout) 385 385 386 386 CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s", 387 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,387 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 388 388 . "ave(X)", zsto,zout) 389 389 390 390 c CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s", 391 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,391 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 392 392 c . "ave(X)", zsto,zout) 393 393 394 394 c CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2", 395 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,395 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 396 396 c . "ave(X)", zsto,zout) 397 397 … … 404 404 c call histdef(nid_mth, "taux", 405 405 c $ "Zonal wind stress", "Pa", 406 c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,406 c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 407 407 c $ "ave(X)", zsto,zout) 408 408 409 409 c call histdef(nid_mth, "tauy", 410 410 c $ "Meridional xind stress", "Pa", 411 c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,411 c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 412 412 c $ "ave(X)", zsto,zout) 413 413 414 414 c CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-", 415 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,415 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 416 416 c . "ave(X)", zsto,zout) 417 417 418 418 c CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-", 419 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,419 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 420 420 c . "ave(X)", zsto,zout) 421 421 -
trunk/LMDZ.TITAN/libf/phytitan/inicondens.F
r102 r1530 10 10 11 11 use dimphy 12 use mod_grid_phy_lmdz, only: nbp_lev 12 13 IMPLICIT NONE 13 #include "dimensions.h"14 14 #include "YOMCST.h" 15 15 … … 17 17 c ----------- 18 18 INTEGER ny 19 real press( llm),temp(llm) ! pressure in mbar !20 REAL yc( llm,ny)19 real press(nbp_lev),temp(nbp_lev) ! pressure in mbar ! 20 REAL yc(nbp_lev,ny) 21 21 character*10 nomy(ny+1) 22 22 … … 28 28 do ic=1,ny 29 29 print*, 'traceur CH(', ic, ')=', nomy(ic),'------------' 30 do l=1, llm30 do l=1,nbp_lev 31 31 32 32 c Par defaut, yc est a 1 c'est a dire qu'on ne condense pas -
trunk/LMDZ.TITAN/libf/phytitan/iophy.F90
r1056 r1530 24 24 25 25 subroutine init_iophy_new(rlat,rlon) 26 USE dimphy 26 USE dimphy, only: klon 27 27 USE mod_phys_lmdz_para 28 USE mod_grid_phy_lmdz 28 USE mod_grid_phy_lmdz, only: nbp_lon, nbp_lat, klon_glo 29 29 USE ioipsl 30 30 implicit none 31 include 'dimensions.h'32 31 real,dimension(klon),intent(in) :: rlon 33 32 real,dimension(klon),intent(in) :: rlat … … 51 50 52 51 !$OMP MASTER 53 ALLOCATE(io_lat(jjm+1-1/(iim*jjm))) 54 io_lat(1)=rlat_glo(1) 55 io_lat(jjm+1-1/(iim*jjm))=rlat_glo(klon_glo) 56 IF ((iim*jjm) > 1) then 57 DO i=2,jjm 58 io_lat(i)=rlat_glo(2+(i-2)*iim) 52 ALLOCATE(io_lat(nbp_lat)) 53 IF (klon_glo == 1) THEN 54 io_lat(1)=rlat_glo(1) 55 ELSE 56 io_lat(1)=rlat_glo(1) 57 io_lat(nbp_lat)=rlat_glo(klon_glo) 58 DO i=2,nbp_lat-1 59 io_lat(i)=rlat_glo(2+(i-2)*nbp_lon) 59 60 ENDDO 60 61 ENDIF 61 62 62 ALLOCATE(io_lon(iim)) 63 io_lon(:)=rlon_glo(2-1/(iim*jjm):iim+1-1/(iim*jjm)) 63 ALLOCATE(io_lon(nbp_lon)) 64 IF (klon_glo == 1) THEN 65 io_lon(1)=rlon_glo(1) 66 ELSE 67 io_lon(1:nbp_lon)=rlon_glo(2:nbp_lon+1) 68 ENDIF 64 69 65 70 ddid=(/ 1,2 /) 66 dsg=(/ iim, jjm+1-1/(iim*jjm)/)67 dsl=(/ iim, jj_nb /)71 dsg=(/ nbp_lon, nbp_lat /) 72 dsl=(/ nbp_lon, jj_nb /) 68 73 dpf=(/ 1,jj_begin /) 69 dpl=(/ iim, jj_end /)74 dpl=(/ nbp_lon, jj_end /) 70 75 dhs=(/ ii_begin-1,0 /) 71 if (mpi_rank==mpi_size-1) then76 IF (mpi_rank==mpi_size-1) THEN 72 77 dhe=(/0,0/) 73 else74 dhe=(/ iim-ii_end,0 /)75 endif78 ELSE 79 dhe=(/ nbp_lon-ii_end,0 /) 80 ENDIF 76 81 77 82 call flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe, & … … 85 90 USE dimphy 86 91 USE mod_phys_lmdz_para 87 use ioipsl 88 implicit none89 i nclude 'dimensions.h'90 real,dimension( iim),intent(in) :: lon91 real,dimension( jjm+1-1/(iim*jjm)),intent(in) :: lat92 use ioipsl, only: flio_dom_set 93 USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat 94 implicit none 95 real,dimension(nbp_lon),intent(in) :: lon 96 real,dimension(nbp_lat),intent(in) :: lat 92 97 93 98 INTEGER,DIMENSION(2) :: ddid … … 100 105 101 106 !$OMP MASTER 102 allocate(io_lat( jjm+1-1/(iim*jjm)))107 allocate(io_lat(nbp_lat)) 103 108 io_lat(:)=lat(:) 104 allocate(io_lon( iim))109 allocate(io_lon(nbp_lon)) 105 110 io_lon(:)=lon(:) 106 111 107 112 ddid=(/ 1,2 /) 108 dsg=(/ iim, jjm+1-1/(iim*jjm)/)109 dsl=(/ iim, jj_nb /)113 dsg=(/ nbp_lon, nbp_lat /) 114 dsl=(/ nbp_lon, jj_nb /) 110 115 dpf=(/ 1,jj_begin /) 111 dpl=(/ iim, jj_end /)116 dpl=(/ nbp_lon, jj_end /) 112 117 dhs=(/ ii_begin-1,0 /) 113 118 if (mpi_rank==mpi_size-1) then 114 119 dhe=(/0,0/) 115 120 else 116 dhe=(/ iim-ii_end,0 /)121 dhe=(/ nbp_lon-ii_end,0 /) 117 122 endif 118 123 … … 125 130 126 131 subroutine histbeg_phy(name,itau0,zjulian,dtime,nhori,nid_day) 127 USE dimphy128 USE mod_ phys_lmdz_para129 use ioipsl 132 USE mod_phys_lmdz_para, only: jj_begin, jj_end, jj_nb, is_sequential 133 USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat 134 use ioipsl, only: histbeg 130 135 use write_field 131 136 implicit none 132 include 'dimensions.h'133 137 134 138 character*(*), intent(IN) :: name … … 141 145 !$OMP MASTER 142 146 if (is_sequential) then 143 call histbeg(name, iim,io_lon, jj_nb,io_lat(jj_begin:jj_end), &144 1, iim,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day)147 call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), & 148 1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day) 145 149 else 146 call histbeg(name, iim,io_lon, jj_nb,io_lat(jj_begin:jj_end), &147 1, iim,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day,phys_domain_id)150 call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), & 151 1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day,phys_domain_id) 148 152 endif 149 153 !$OMP END MASTER … … 154 158 plon,plat,plon_bounds,plat_bounds, & 155 159 nname,itau0,zjulian,dtime,nnhori,nnid_day) 156 USE dimphy 160 USE dimphy, only: klon 157 161 USE mod_phys_lmdz_para 158 USE mod_grid_phy_lmdz 159 use ioipsl 160 use write_field 161 implicit none 162 include 'dimensions.h' 162 USE mod_grid_phy_lmdz, only: klon_glo, nbp_lon, nbp_lat 163 use ioipsl, only: histbeg 164 165 implicit none 163 166 164 167 real,dimension(klon),intent(in) :: rlon … … 186 189 REAL, allocatable, dimension(:) :: npplat, npplon 187 190 REAL, allocatable, dimension(:,:) :: npplat_bounds, npplon_bounds 188 INTEGER, PARAMETER :: jjmp1=jjm+1-1/jjm 189 REAL, dimension(iim,jjmp1) :: zx_lon, zx_lat 191 REAL, dimension(nbp_lon,nbp_lat) :: zx_lon, zx_lat 190 192 191 193 CALL gather(rlat,rlat_glo) … … 212 214 endif 213 215 ! 214 IF ( tabij(i).LE. iim) THEN216 IF ( tabij(i).LE.nbp_lon) THEN 215 217 plat_bounds(i,1)=rlat_glo(tabij(i)) 216 218 ELSE 217 plat_bounds(i,1)=rlat_glo(tabij(i)- iim)219 plat_bounds(i,1)=rlat_glo(tabij(i)-nbp_lon) 218 220 ENDIF 219 plat_bounds(i,2)=rlat_glo(tabij(i)+ iim)221 plat_bounds(i,2)=rlat_glo(tabij(i)+nbp_lon) 220 222 ! 221 223 ! print*,'CFMIP_iophy point i lon lon_bds',i,plon_bounds(i,1),rlon_glo(tabij(i)),plon_bounds(i,2) … … 233 235 ENDDO 234 236 235 CALL gr_fi_ecrit(1,klon, iim,jjmp1,rlon_glo,zx_lon)236 if (( iim*jjm).gt.1) then237 DO i = 1, iim237 CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlon_glo,zx_lon) 238 if ((nbp_lon*nbp_lat).gt.1) then 239 DO i = 1, nbp_lon 238 240 zx_lon(i,1) = rlon_glo(i+1) 239 zx_lon(i, jjmp1) = rlon_glo(i+1)241 zx_lon(i,nbp_lat) = rlon_glo(i+1) 240 242 ENDDO 241 243 endif 242 CALL gr_fi_ecrit(1,klon, iim,jjmp1,rlat_glo,zx_lat)244 CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlat_glo,zx_lat) 243 245 244 246 DO i=1,pim … … 249 251 250 252 if (ipt(i).EQ.1) then 251 plon_bounds(i,1)=zx_lon( iim,jpt(i))253 plon_bounds(i,1)=zx_lon(nbp_lon,jpt(i)) 252 254 plon_bounds(i,2)=360.+zx_lon(ipt(i)+1,jpt(i)) 253 255 endif 254 256 255 if (ipt(i).EQ. iim) then257 if (ipt(i).EQ.nbp_lon) then 256 258 plon_bounds(i,2)=360.+zx_lon(1,jpt(i)) 257 259 endif … … 265 267 endif 266 268 267 if (jpt(i).EQ. jjmp1) then268 plat_bounds(i,1)=zx_lat(ipt(i), jjmp1)+0.001269 plat_bounds(i,2)=zx_lat(ipt(i), jjmp1)-0.001269 if (jpt(i).EQ.nbp_lat) then 270 plat_bounds(i,1)=zx_lat(ipt(i),nbp_lat)+0.001 271 plat_bounds(i,2)=zx_lat(ipt(i),nbp_lat)-0.001 270 272 endif 271 273 ! … … 323 325 324 326 subroutine histwrite2d_phy(nid,lpoint,name,itau,field) 325 USE dimphy 327 USE dimphy, only: klon 326 328 USE mod_phys_lmdz_para 327 USE ioipsl 328 implicit none329 i nclude 'dimensions.h'329 USE ioipsl, only: histwrite 330 USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat 331 implicit none 330 332 331 333 integer,intent(in) :: nid … … 336 338 REAL,dimension(klon_mpi) :: buffer_omp 337 339 INTEGER, allocatable, dimension(:) :: index2d 338 REAL :: Field2d( iim,jj_nb)340 REAL :: Field2d(nbp_lon,jj_nb) 339 341 340 342 integer :: ip … … 347 349 CALL grid1Dto2D_mpi(buffer_omp,Field2d) 348 350 if(.NOT.lpoint) THEN 349 ALLOCATE(index2d( iim*jj_nb))350 ALLOCATE(fieldok( iim*jj_nb))351 CALL histwrite(nid,name,itau,Field2d, iim*jj_nb,index2d)351 ALLOCATE(index2d(nbp_lon*jj_nb)) 352 ALLOCATE(fieldok(nbp_lon*jj_nb)) 353 CALL histwrite(nid,name,itau,Field2d,nbp_lon*jj_nb,index2d) 352 354 else 353 355 ALLOCATE(fieldok(npstn)) … … 378 380 379 381 subroutine histwrite3d_phy(nid,lpoint,name,itau,field) 380 USE dimphy 382 USE dimphy, only: klon 381 383 USE mod_phys_lmdz_para 382 383 use ioipsl 384 implicit none 385 include 'dimensions.h' 384 USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat 385 use ioipsl, only: histwrite 386 implicit none 386 387 387 388 integer,intent(in) :: nid … … 391 392 real,dimension(:,:),intent(in) :: field ! --> field(klon,:) 392 393 REAL,dimension(klon_mpi,size(field,2)) :: buffer_omp 393 REAL :: Field3d( iim,jj_nb,size(field,2))394 REAL :: Field3d(nbp_lon,jj_nb,size(field,2)) 394 395 INTEGER :: ip, n, nlev 395 396 INTEGER, ALLOCATABLE, dimension(:) :: index3d … … 409 410 CALL grid1Dto2D_mpi(buffer_omp,field3d) 410 411 if(.NOT.lpoint) THEN 411 ALLOCATE(index3d( iim*jj_nb*nlev))412 ALLOCATE(fieldok( iim*jj_nb,nlev))413 CALL histwrite(nid,name,itau,Field3d, iim*jj_nb*nlev,index3d)412 ALLOCATE(index3d(nbp_lon*jj_nb*nlev)) 413 ALLOCATE(fieldok(nbp_lon*jj_nb,nlev)) 414 CALL histwrite(nid,name,itau,Field3d,nbp_lon*jj_nb*nlev,index3d) 414 415 else 415 416 nlev=size(field,2) -
trunk/LMDZ.TITAN/libf/phytitan/orodrag.F
r1056 r1530 10 10 use dimphy 11 11 IMPLICIT NONE 12 13 #include "dimensions.h"14 #include "paramet.h"15 12 16 13 c -
trunk/LMDZ.TITAN/libf/phytitan/orosetup.F
r1056 r1530 97 97 use dimphy 98 98 implicit none 99 100 #include "dimensions.h"101 #include "paramet.h"102 99 103 100 #include "YOMCST.h" -
trunk/LMDZ.TITAN/libf/phytitan/phyetat0.F90
r1525 r1530 21 21 ! Objet: Lecture de l'etat initial pour la physique 22 22 !====================================================================== 23 #include "dimensions.h"24 23 #include "netcdf.inc" 25 24 #include "dimsoil.h" -
trunk/LMDZ.TITAN/libf/phytitan/phyredem.F90
r1525 r1530 20 20 ! Objet: Lecture de l'etat initial pour la physique 21 21 !====================================================================== 22 #include "dimensions.h"23 22 #include "netcdf.inc" 24 23 #include "dimsoil.h" -
trunk/LMDZ.TITAN/libf/phytitan/physiq.F
r1525 r1530 73 73 USE time_phylmdz_mod, only: itau_phy,day_ref,annee_ref,nday 74 74 USE logic_mod, only: iflag_trac,moyzon_ch,moyzon_mu 75 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, nbp_lev 75 76 IMPLICIT none 76 77 c====================================================================== -
trunk/LMDZ.TITAN/libf/phytitan/radlwsw.F
r1356 r1530 41 41 USE write_field_phy 42 42 IMPLICIT none 43 #include "dimensions.h"44 43 #include "YOMCST.h" 45 44 #include "clesphys.h" -
trunk/LMDZ.TITAN/libf/phytitan/soil.F
r102 r1530 48 48 use dimphy 49 49 IMPLICIT NONE 50 #include "dimensions.h"51 50 #include "YOMCST.h" 52 51 #include "dimsoil.h" -
trunk/LMDZ.TITAN/libf/phytitan/sources.F
r808 r1530 23 23 use dimphy 24 24 IMPLICIT NONE 25 #include "dimensions.h"26 #include "microtab.h"27 25 #include "YOMCST.h" 28 26 c -
trunk/LMDZ.TITAN/libf/phytitan/sugwd.F
r1056 r1530 76 76 use dimphy 77 77 IMPLICIT NONE 78 79 #include "dimensions.h"80 #include "paramet.h"81 78 82 79 #include "YOEGWD.h" -
trunk/LMDZ.TITAN/libf/phytitan/tgmdat_mod.F90
r1461 r1530 1 1 MODULE TGMDAT_MOD 2 2 IMPLICIT NONE 3 include "dimensions.h"4 3 5 INTEGER,PARAMETER :: NLAYER=llm,NLEVEL=NLAYER+16 4 REAL,SAVE :: PI=3.14159265358979323846 7 5 ! RGAS IS THE UNIVERSAL GAS CONSTANT IN UNITS OF: M SEC-2 AMU K-1 KM -
trunk/LMDZ.TITAN/libf/phytitan/ustarhb.F
r495 r1530 18 18 use dimphy 19 19 IMPLICIT none 20 #include "dimensions.h"21 20 #include "YOMCST.h" 22 21 c -
trunk/LMDZ.TITAN/libf/phytitan/vdif_kcay.F
r102 r1530 8 8 use dimphy 9 9 IMPLICIT NONE 10 #include "dimensions.h"11 10 c....................................................................... 12 11 c -
trunk/LMDZ.TITAN/libf/phytitan/yamada.F
r102 r1530 8 8 use dimphy 9 9 IMPLICIT NONE 10 #include "dimensions.h"11 10 c....................................................................... 12 11 c -
trunk/LMDZ.TITAN/libf/phytitan/yamada4.F
r814 r1530 9 9 use dimphy 10 10 IMPLICIT NONE 11 #include "dimensions.h"12 11 #include "YOMCST.h" 13 12 c....................................................................... -
trunk/LMDZ.TITAN/libf/phytitan/zenang.F
r119 r1530 26 26 use dimphy 27 27 IMPLICIT none 28 #include "dimensions.h"29 28 #include "YOMCST.h" 30 29 #include "comorbit.h"
Note: See TracChangeset
for help on using the changeset viewer.
