Changeset 1530 for trunk/LMDZ.VENUS/libf/phyvenus
- Timestamp:
- Apr 6, 2016, 8:51:34 AM (9 years ago)
- Location:
- trunk/LMDZ.VENUS/libf/phyvenus
- Files:
-
- 3 deleted
- 61 edited
-
aaam_bud.F (modified) (13 diffs)
-
ajsec.F (modified) (4 diffs)
-
ballon.F (modified) (13 diffs)
-
blendrad.F (modified) (1 diff)
-
clmain.F (modified) (12 diffs)
-
cltrac.F (modified) (1 diff)
-
coefkzmin.F (modified) (1 diff)
-
concentrations2.F (modified) (1 diff)
-
conduction.F (modified) (1 diff)
-
conf_dat3d.F (deleted)
-
diagphy.F (modified) (2 diffs)
-
drag_noro.F (modified) (1 diff)
-
euvheat.F90 (modified) (1 diff)
-
flott_gwd_ran.F90 (modified) (1 diff)
-
grid_noro.F (modified) (5 diffs)
-
gwprofil.F (modified) (1 diff)
-
gwstress.F (modified) (1 diff)
-
hgardfou.F (modified) (1 diff)
-
hrtherm.F (modified) (1 diff)
-
ini_histday.h (modified) (8 diffs)
-
ini_histins.h (modified) (22 diffs)
-
ini_histmth.h (modified) (8 diffs)
-
iniinterp_horiz.F (deleted)
-
interp_horiz.F (deleted)
-
iophy.F90 (modified) (19 diffs)
-
jthermcalc.F (modified) (4 diffs)
-
jthermcalc_e107.F (modified) (1 diff)
-
load_ksi.F (modified) (1 diff)
-
lw_venus_ve.F (modified) (1 diff)
-
lwi.F (modified) (1 diff)
-
moldiff_red.F90 (modified) (5 diffs)
-
moldiffcoeff_red.F (modified) (1 diff)
-
molvis.F (modified) (1 diff)
-
new_cloud_sedim.F (modified) (1 diff)
-
nirco2abs.F (modified) (1 diff)
-
nlte_tcool.F (modified) (1 diff)
-
nltecool.F (modified) (1 diff)
-
nlthermeq.F (modified) (1 diff)
-
orodrag.F (modified) (1 diff)
-
orosetup.F (modified) (1 diff)
-
param_read.F (modified) (1 diff)
-
param_read_e107.F (modified) (1 diff)
-
phyetat0.F90 (modified) (1 diff)
-
phyredem.F90 (modified) (1 diff)
-
physiq.F (modified) (2 diffs)
-
phytrac_chimie.F (modified) (1 diff)
-
phytrac_emiss.F (modified) (2 diffs)
-
phytrac_relax.F (modified) (1 diff)
-
radlwsw.F (modified) (1 diff)
-
radlwsw.NewtonCool (modified) (1 diff)
-
soil.F (modified) (1 diff)
-
sugwd.F (modified) (1 diff)
-
sw_venus_cl.F (modified) (1 diff)
-
sw_venus_cl_1Dglobave.F (modified) (1 diff)
-
sw_venus_dc.F (modified) (1 diff)
-
sw_venus_dc_1Dglobave.F (modified) (1 diff)
-
sw_venus_ve.F (modified) (1 diff)
-
sw_venus_ve_1Dglobave.F (modified) (1 diff)
-
transp.F (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.VENUS/libf/phyvenus/aaam_bud.F
r101 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 100 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.VENUS/libf/phyvenus/ajsec.F
r1301 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.VENUS/libf/phyvenus/ballon.F
r1442 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 … … 258 260 259 261 use dimphy 262 use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat 260 263 implicit none 261 264 … … 285 288 c====================================================================== 286 289 287 #include "dimensions.h"288 290 #include "YOMCST.h" 289 291 c … … 291 293 c 292 294 real map_u(klon),map_v(klon),map_a(klon) 293 real latit( jjm+1),longit(iim)295 real latit(nbp_lat),longit(nbp_lon) 294 296 real phi,lam,ubal,vbal,abal 295 297 c … … 298 300 INTEGER i,j,k 299 301 INTEGER jj,ii 300 REAL ujj( iim+1),vjj(iim+1),ajj(iim+1)302 REAL ujj(nbp_lon+1),vjj(nbp_lon+1),ajj(nbp_lon+1) 301 303 REAL factlat,factlon 302 304 … … 304 306 c------------------------------------------------- 305 307 jj=1 ! POLE NORD 306 do j=2, jjm308 do j=2,nbp_lat-1 307 309 if (phi.lt.latit(j)) jj=j 308 310 enddo … … 311 313 c pole nord 312 314 if (jj.eq.1) then 313 do i=1, iim315 do i=1,nbp_lon 314 316 ujj(i) = map_u(i+1)*factlat + map_u(1)*(1-factlat) 315 317 vjj(i) = map_v(i+1)*factlat + map_v(1)*(1-factlat) … … 317 319 enddo 318 320 c pole sud 319 elseif (jj.eq. jjm) then320 do i=1, iim321 k = (jj-2)* iim+1+i321 elseif (jj.eq.nbp_lat-1) then 322 do i=1,nbp_lon 323 k = (jj-2)*nbp_lon+1+i 322 324 ujj(i) = map_u(klon)*factlat + map_u(k)*(1-factlat) 323 325 vjj(i) = map_v(klon)*factlat + map_v(k)*(1-factlat) … … 326 328 c autres latitudes 327 329 else 328 do i=1, iim329 k = (jj-2)* iim+1+i330 ujj(i) = map_u(k+ iim)*factlat + map_u(k)*(1-factlat)331 vjj(i) = map_v(k+ iim)*factlat + map_v(k)*(1-factlat)332 ajj(i) = map_a(k+ iim)*factlat + map_a(k)*(1-factlat)330 do i=1,nbp_lon 331 k = (jj-2)*nbp_lon+1+i 332 ujj(i) = map_u(k+nbp_lon)*factlat + map_u(k)*(1-factlat) 333 vjj(i) = map_v(k+nbp_lon)*factlat + map_v(k)*(1-factlat) 334 ajj(i) = map_a(k+nbp_lon)*factlat + map_a(k)*(1-factlat) 333 335 enddo 334 336 endif 335 ujj( iim+1)=ujj(1)336 vjj( iim+1)=vjj(1)337 ajj( iim+1)=ajj(1)337 ujj(nbp_lon+1)=ujj(1) 338 vjj(nbp_lon+1)=vjj(1) 339 ajj(nbp_lon+1)=ajj(1) 338 340 339 341 c Interpolation in longitudes 340 342 c------------------------------------------------- 341 343 ii=1 ! lon=-180 342 do i=2, iim344 do i=2,nbp_lon 343 345 if (lam.gt.longit(i)) ii=i 344 346 enddo … … 372 374 c====================================================================== 373 375 374 #include "dimensions.h"375 376 #include "YOMCST.h" 376 377 c -
trunk/LMDZ.VENUS/libf/phyvenus/blendrad.F
r1310 r1530 15 15 use dimphy 16 16 implicit none 17 #include "dimensions.h"18 17 c#include "dimradmars.h" 19 18 #include "nlteparams.h" -
trunk/LMDZ.VENUS/libf/phyvenus/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.VENUS/libf/phyvenus/cltrac.F
r1442 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.VENUS/libf/phyvenus/coefkzmin.F
r101 r1530 9 9 IMPLICIT NONE 10 10 11 #include "dimensions.h"12 11 #include "YOMCST.h" 13 12 -
trunk/LMDZ.VENUS/libf/phyvenus/concentrations2.F
r1452 r1530 22 22 ! declarations 23 23 24 #include "dimensions.h"25 24 #include "YOMCST.h" 26 25 #include "clesphys.h" -
trunk/LMDZ.VENUS/libf/phyvenus/conduction.F
r1310 r1530 17 17 c declarations: 18 18 c----------------------------------------------------------------------- 19 20 !#include "dimensions.h"21 !#include "dimphys.h"22 !#include "comcstfi.h"23 !#include "surfdat.h"24 !#include "chimiedata.h"25 !#include "conc.h"26 19 27 20 c arguments: -
trunk/LMDZ.VENUS/libf/phyvenus/diagphy.F
r1017 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.VENUS/libf/phyvenus/drag_noro.F
r1301 r1530 11 11 use dimphy 12 12 IMPLICIT none 13 14 #include "dimensions.h"15 #include "paramet.h"16 13 17 14 c====================================================================== -
trunk/LMDZ.VENUS/libf/phyvenus/euvheat.F90
r1442 r1530 32 32 ! ------------------ 33 33 ! 34 !#include "dimensions.h"35 34 #include "YOMCST.h" 36 35 #include "clesphys.h" -
trunk/LMDZ.VENUS/libf/phyvenus/flott_gwd_ran.F90
r826 r1530 13 13 use dimphy 14 14 implicit none 15 16 #include "dimensions.h"17 #include "paramet.h"18 15 19 16 #include "YOEGWD.h" -
trunk/LMDZ.VENUS/libf/phyvenus/grid_noro.F
r1301 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.VENUS/libf/phyvenus/gwprofil.F
r808 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.VENUS/libf/phyvenus/gwstress.F
r101 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.VENUS/libf/phyvenus/hgardfou.F
r101 r1530 9 9 c Verifier la temperature 10 10 c====================================================================== 11 #include "dimensions.h"12 11 #include "YOMCST.h" 13 12 REAL t(klon,klev), tsol(klon) -
trunk/LMDZ.VENUS/libf/phyvenus/hrtherm.F
r1442 r1530 15 15 16 16 17 #include "dimensions.h"18 17 #include "param.h" 19 18 #include "param_v4.h" -
trunk/LMDZ.VENUS/libf/phyvenus/ini_histday.h
r902 r1530 23 23 c 24 24 CALL histdef(nid_day, "phis", "Surface geop. height", "-", 25 . iim,jj_nb,nhori, 1,1,1, nvert, 32,25 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 26 26 . "once", zsto,zout) 27 27 c 28 28 CALL histdef(nid_day, "aire", "Grid area", "-", 29 . iim,jj_nb,nhori, 1,1,1, nvert, 32,29 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 30 30 . "once", zsto,zout) 31 31 c 32 32 CALL histdef(nid_day, "tsol", "Surface Temperature", "K", 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 c 36 36 CALL histdef(nid_day, "psol", "Surface Pressure", "Pa", 37 . iim,jj_nb,nhori, 1,1,1, nvert, 32,37 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 38 38 . "ave(X)", zsto,zout) 39 39 c 40 40 c CALL histdef(nid_day, "ue", "Zonal energy transport", "-", 41 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,41 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 42 42 c . "ave(X)", zsto,zout) 43 43 c 44 44 c CALL histdef(nid_day, "ve", "Merid energy transport", "-", 45 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,45 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 46 46 c . "ave(X)", zsto,zout) 47 47 c 48 48 c CALL histdef(nid_day, "cdragh", "Drag coef on T", "-", 49 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,49 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 50 50 c . "ave(X)", zsto,zout) 51 51 c 52 52 c CALL histdef(nid_day, "cdragm", "Drag coef on U", "-", 53 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,53 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 54 54 c . "ave(X)", zsto,zout) 55 55 c … … 62 62 c 63 63 CALL histdef(nid_day, "temp", "Air temperature", "K", 64 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,64 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 65 65 . "ave(X)", zsto,zout) 66 66 c 67 67 CALL histdef(nid_day, "pres", "Air pressure", "Pa", 68 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,68 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 69 69 . "ave(X)", zsto,zout) 70 70 c 71 71 CALL histdef(nid_day, "geop", "Geopotential height", "m", 72 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,72 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 73 73 . "ave(X)", zsto,zout) 74 74 c 75 75 CALL histdef(nid_day, "vitu", "Zonal wind", "m/s", 76 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,76 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 77 77 . "ave(X)", zsto,zout) 78 78 c 79 79 CALL histdef(nid_day, "vitv", "Meridional wind", "m/s", 80 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,80 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 81 81 . "ave(X)", zsto,zout) 82 82 c 83 83 CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s", 84 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,84 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 85 85 . "ave(X)", zsto,zout) 86 86 c 87 87 CALL histdef(nid_day, "dudyn", "Dynamics dU", "m/s2", 88 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,88 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 89 89 . "ave(X)", zsto,zout) 90 90 c 91 91 CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2", 92 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,92 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 93 93 . "ave(X)", zsto,zout) 94 94 c 95 95 c CALL histdef(nid_day, "mang", "Angular momentum", "kg m2/s", 96 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,96 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 97 97 c . "ave(X)", zsto,zout) 98 98 c 99 99 c CALL histdef(nid_day, "Kz", "vertical diffusion coef", "m2/s", 100 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,100 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 101 101 c . "ave(X)", zsto,zout) 102 102 c … … 107 107 WRITE(str2,'(i2.2)') iq 108 108 CALL histdef(nid_day, tname(iq), ttext(iq), "ppm", 109 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,109 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 110 110 . "ave(X)", zsto,zout) 111 111 ELSE … … 117 117 c 118 118 CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2", 119 . iim,jj_nb,nhori, 1,1,1, nvert, 32,119 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 120 120 . "ave(X)", zsto,zout) 121 121 c … … 130 130 c 131 131 CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2", 132 . iim,jj_nb,nhori, 1,1,1, nvert, 32,132 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 133 133 . "ave(X)", zsto,zout) 134 134 c 135 135 CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2", 136 . iim,jj_nb,nhori, 1,1,1, nvert, 32,136 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 137 137 . "ave(X)", zsto,zout) 138 138 c 139 139 CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2", 140 . iim,jj_nb,nhori, 1,1,1, nvert, 32,140 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 141 141 . "ave(X)", zsto,zout) 142 142 c … … 144 144 c 145 145 CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2", 146 . iim,jj_nb,nhori, klev,1,klev,nvert,146 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 147 147 . 32, "ave(X)", zsto,zout) 148 148 c 149 149 CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2", 150 . iim,jj_nb,nhori, klev,1,klev,nvert,150 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 151 151 . 32, "ave(X)", zsto,zout) 152 152 c 153 153 CALL histdef(nid_day, "fluxvdf", "PBL net flux","W/m2", 154 . iim,jj_nb,nhori, klev,1,klev,nvert,154 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 155 155 . 32, "ave(X)", zsto,zout) 156 156 c 157 157 CALL histdef(nid_day, "fluxdyn", "Dyn. net flux","W/m2", 158 . iim,jj_nb,nhori, klev,1,klev,nvert,158 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 159 159 . 32, "ave(X)", zsto,zout) 160 160 c 161 161 CALL histdef(nid_day, "fluxajs", "Dry adj. net flux","W/m2", 162 . iim,jj_nb,nhori, klev,1,klev,nvert,162 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 163 163 . 32, "ave(X)", zsto,zout) 164 164 c 165 165 c CALL histdef(nid_day, "fluxec", "Cin. net flux","W/m2", 166 c . iim,jj_nb,nhori, klev,1,klev,nvert,166 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 167 167 c . 32, "ave(X)", zsto,zout) 168 168 c … … 173 173 c 174 174 CALL histdef(nid_day, "dtdyn", "Dynamics dT", "K/s", 175 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,175 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 176 176 . "ave(X)", zsto,zout) 177 177 c 178 178 c CALL histdef(nid_day, "dtphy", "Physics dT", "K/s", 179 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,179 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 180 180 c . "ave(X)", zsto,zout) 181 181 c 182 182 CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s", 183 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,183 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 184 184 . "ave(X)", zsto,zout) 185 185 c 186 186 CALL histdef(nid_day, "dtajs", "Dry adjust. dT", "K/s", 187 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,187 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 188 188 . "ave(X)", zsto,zout) 189 189 c 190 190 CALL histdef(nid_day, "dtswr", "SW radiation dT", "K/s", 191 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,191 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 192 192 . "ave(X)", zsto,zout) 193 193 c 194 194 CALL histdef(nid_day, "dtlwr", "LW radiation dT", "K/s", 195 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,195 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 196 196 . "ave(X)", zsto,zout) 197 197 c 198 198 c CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s", 199 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,199 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 200 200 c . "ave(X)", zsto,zout) 201 201 c 202 202 CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2", 203 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,203 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 204 204 . "ave(X)", zsto,zout) 205 205 c 206 206 CALL histdef(nid_mth, "dugwo", "GW oro dU", "m/s2", 207 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,207 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 208 208 . "ave(X)", zsto,zout) 209 209 c 210 210 CALL histdef(nid_mth, "dugwno", "GW non-oro dU", "m/s2", 211 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,211 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 212 212 . "ave(X)", zsto,zout) 213 213 c 214 214 c CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2", 215 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,215 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 216 216 c . "ave(X)", zsto,zout) 217 217 c … … 223 223 c call histdef(nid_day, "taux", 224 224 c $ "Zonal wind stress", "Pa", 225 c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,225 c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 226 226 c $ "ave(X)", zsto,zout) 227 227 c 228 228 c call histdef(nid_day, "tauy", 229 229 c $ "Meridional xind stress", "Pa", 230 c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,230 c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 231 231 c $ "ave(X)", zsto,zout) 232 232 c 233 233 c CALL histdef(nid_day, "cdrm", "Momentum drag coef.", "-", 234 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,234 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 235 235 c . "ave(X)", zsto,zout) 236 236 c 237 237 c CALL histdef(nid_day, "cdrh", "Heat drag coef.", "-", 238 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,238 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 239 239 c . "ave(X)", zsto,zout) 240 240 c -
trunk/LMDZ.VENUS/libf/phyvenus/ini_histins.h
r1468 r1530 23 23 c 24 24 CALL histdef(nid_ins, "phis", "Surface geop. height", "-", 25 . iim,jj_nb,nhori, 1,1,1, nvert, 32,25 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 26 26 . "once", zsto,zout) 27 27 c 28 28 CALL histdef(nid_ins, "aire", "Grid area", "-", 29 . iim,jj_nb,nhori, 1,1,1, nvert, 32,29 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 30 30 . "once", zsto,zout) 31 31 c 32 32 CALL histdef(nid_ins, "tsol", "Surface Temperature", "K", 33 . iim,jj_nb,nhori, 1,1,1, nvert, 32,33 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 34 34 . "ins(X)", zsto,zout) 35 35 c 36 36 CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa", 37 . iim,jj_nb,nhori, 1,1,1, nvert, 32,37 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 38 38 . "ins(X)", zsto,zout) 39 39 c 40 40 c CALL histdef(nid_ins, "ue", "Zonal energy transport", "-", 41 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,41 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 42 42 c . "ins(X)", zsto,zout) 43 43 c 44 44 c CALL histdef(nid_ins, "ve", "Merid energy transport", "-", 45 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,45 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 46 46 c . "ins(X)", zsto,zout) 47 47 c 48 48 c CALL histdef(nid_ins, "cdragh", "Drag coef on T", "-", 49 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,49 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 50 50 c . "ins(X)", zsto,zout) 51 51 c 52 52 c CALL histdef(nid_ins, "cdragm", "Drag coef on U", "-", 53 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,53 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 54 54 c . "ins(X)", zsto,zout) 55 55 c … … 62 62 c 63 63 CALL histdef(nid_ins, "temp", "Air temperature", "K", 64 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,64 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 65 65 . "ins(X)", zsto,zout) 66 66 c 67 67 CALL histdef(nid_ins, "pres", "Air pressure", "Pa", 68 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,68 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 69 69 . "ins(X)", zsto,zout) 70 70 c 71 71 CALL histdef(nid_ins, "geop", "Geopotential height", "m", 72 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,72 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 73 73 . "ins(X)", zsto,zout) 74 74 c 75 75 CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s", 76 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,76 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 77 77 . "ins(X)", zsto,zout) 78 78 c 79 79 CALL histdef(nid_ins, "vitv", "Meridional wind", "m/s", 80 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,80 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 81 81 . "ins(X)", zsto,zout) 82 82 c 83 83 CALL histdef(nid_ins, "vitw", "Vertical wind", "Pa/s", 84 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,84 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 85 85 . "ins(X)", zsto,zout) 86 86 c 87 87 CALL histdef(nid_ins, "dudyn", "Dynamics dU", "m/s2", 88 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,88 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 89 89 . "ins(X)", zsto,zout) 90 90 c 91 91 CALL histdef(nid_ins, "duvdf", "Boundary-layer dU", "m/s2", 92 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,92 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 93 93 . "ins(X)", zsto,zout) 94 94 c 95 95 c CALL histdef(nid_ins, "mang", "Angular momentum", "kg m2/s", 96 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,96 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 97 97 c . "ins(X)", zsto,zout) 98 98 c 99 99 c CALL histdef(nid_ins, "Kz", "vertical diffusion coef", "m2/s", 100 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,100 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 101 101 c . "ins(X)", zsto,zout) 102 102 c 103 103 CALL histdef(nid_ins, "mmean", "Mean molecular mass", "g/mol", 104 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,104 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 105 105 . "ins(X)", zsto,zout) 106 106 107 107 CALL histdef(nid_ins, "rho", "Air density [mass/Vol]", "kg/m3", 108 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,108 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 109 109 . "ins(X)", zsto,zout) 110 110 … … 116 116 WRITE(str2,'(i2.2)') iq 117 117 CALL histdef(nid_ins, tname(iq), ttext(iq), "mol/mol", 118 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,118 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 119 119 . "ins(X)", zsto,zout) 120 120 ELSE … … 127 127 IF (callthermos .and. ok_chem) THEN 128 128 CALL histdef(nid_ins, "d_qmoldif CO2", "Dif molec" , "kg/kg", 129 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,129 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 130 130 . "ins(X)", zsto,zout) 131 131 CALL histdef(nid_ins, "d_qmoldif O3p", "Dif molec" , "kg/kg", 132 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,132 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 133 133 . "ins(X)", zsto,zout) 134 134 CALL histdef(nid_ins, "d_qmoldif N2", "Dif molec" , "kg/kg", 135 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,135 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 136 136 . "ins(X)", zsto,zout) 137 137 ENDIF 138 138 c 139 139 CALL histdef(nid_ins, "tops", "Solar rad. at TOA", "W/m2", 140 . iim,jj_nb,nhori, 1,1,1, nvert, 32,140 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 141 141 . "ins(X)", zsto,zout) 142 142 c … … 147 147 c 148 148 CALL histdef(nid_ins, "NBRTOTm1", "Nbr total droplet", 149 . "#/cm3", iim,jj_nb,nhori, klev,1,klev, nvert, 32,149 . "#/cm3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 150 150 . "ins(X)", zsto,zout) 151 151 c … … 155 155 c . for log normal distribution" , 156 156 c . "fraction", 157 c . iim,jj_nb,nhori, klev,1,klev, nvert, 32,157 c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32, 158 158 c . "ins(X)", zsto,zout) 159 159 c … … 163 163 c . for log normal distribution", 164 164 c . "fraction", 165 c . iim,jj_nb,nhori, klev,1,klev, nvert, 32,165 c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32, 166 166 c . "ins(X)", zsto,zout) 167 167 c … … 171 171 c 172 172 CALL histdef(nid_ins, "NBRTOTm2", "Nbr total droplet", 173 . "#/cm3", iim,jj_nb,nhori, klev,1,klev, nvert, 32,173 . "#/cm3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 174 174 . "ins(X)", zsto,zout) 175 175 c … … 179 179 c . for log normal distribution" , 180 180 c . "fraction", 181 c . iim,jj_nb,nhori, klev,1,klev, nvert, 32,181 c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32, 182 182 c . "ins(X)", zsto,zout) 183 183 c … … 187 187 c . for log normal distribution", 188 188 c . "fraction", 189 c . iim,jj_nb,nhori, klev,1,klev, nvert, 32,189 c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32, 190 190 c . "ins(X)", zsto,zout) 191 191 c … … 195 195 c 196 196 CALL histdef(nid_ins, "NBRTOTm3", "Nbr total droplet", "#/cm3", 197 . iim,jj_nb,nhori, klev,1,klev, nvert, 32,197 . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32, 198 198 . "ins(X)", zsto,zout) 199 199 c … … 203 203 c . for log normal distribution" , 204 204 c . "fraction", 205 c . iim,jj_nb,nhori, klev,1,klev, nvert, 32,205 c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32, 206 206 c . "ins(X)", zsto,zout) 207 207 c … … 211 211 c . for log normal distribution", 212 212 c . "fraction", 213 c . iim,jj_nb,nhori, klev,1,klev, nvert, 32,213 c . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32, 214 214 c . "ins(X)", zsto,zout) 215 215 c … … 222 222 CALL histdef(nid_ins, "WH2SO4", "Weight fraction H2SO4", 223 223 . "fraction", 224 . iim,jj_nb,nhori, klev,1,klev, nvert, 32,224 . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32, 225 225 . "ins(X)", zsto,zout) 226 226 c … … 229 229 CALL histdef(nid_ins, "rho_droplet", "density cloud droplet", 230 230 . "kg.m-3", 231 . iim,jj_nb,nhori, klev,1,klev, nvert, 32,231 . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32, 232 232 . "ins(X)", zsto,zout) 233 233 c … … 239 239 CALL histdef(nid_ins, "d_tr_sed_H2SO4", "var mmr from sedim", 240 240 . "kg/kg", 241 . iim,jj_nb,nhori, klev,1,klev, nvert, 32,241 . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32, 242 242 . "ins(X)", zsto,zout) 243 243 c … … 246 246 CALL histdef(nid_ins, "d_tr_sed_H2O", "var mmr from sedim", 247 247 . "kg/kg", 248 . iim,jj_nb,nhori, klev,1,klev, nvert, 32,248 . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32, 249 249 . "ins(X)", zsto,zout) 250 250 c … … 253 253 CALL histdef(nid_ins, "F_sedim", "tendency from sedim", 254 254 . "kg.m-2.s-1", 255 . iim,jj_nb,nhori, klev,1,klev, nvert, 32,255 . nbp_lon,jj_nb,nhori, klev,1,klev, nvert, 32, 256 256 . "ins(X)", zsto,zout) 257 257 c … … 267 267 c 268 268 CALL histdef(nid_ins, "topl", "IR rad. at TOA", "W/m2", 269 . iim,jj_nb,nhori, 1,1,1, nvert, 32,269 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 270 270 . "ins(X)", zsto,zout) 271 271 c 272 272 CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2", 273 . iim,jj_nb,nhori, 1,1,1, nvert, 32,273 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 274 274 . "ins(X)", zsto,zout) 275 275 c 276 276 CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2", 277 . iim,jj_nb,nhori, 1,1,1, nvert, 32,277 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 278 278 . "ins(X)", zsto,zout) 279 279 c … … 281 281 c 282 282 CALL histdef(nid_ins, "SWnet", "Net SW flux","W/m2", 283 . iim,jj_nb,nhori, klev,1,klev,nvert,283 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 284 284 . 32, "ins(X)", zsto,zout) 285 285 c 286 286 CALL histdef(nid_ins, "LWnet", "Net LW flux","W/m2", 287 . iim,jj_nb,nhori, klev,1,klev,nvert,287 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 288 288 . 32, "ins(X)", zsto,zout) 289 289 c 290 290 CALL histdef(nid_ins, "fluxvdf", "PBL net flux","W/m2", 291 . iim,jj_nb,nhori, klev,1,klev,nvert,291 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 292 292 . 32, "ins(X)", zsto,zout) 293 293 c 294 294 CALL histdef(nid_ins, "fluxdyn", "Dyn. net flux","W/m2", 295 . iim,jj_nb,nhori, klev,1,klev,nvert,295 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 296 296 . 32, "ins(X)", zsto,zout) 297 297 c 298 298 CALL histdef(nid_ins, "fluxajs", "Dry adj. net flux","W/m2", 299 . iim,jj_nb,nhori, klev,1,klev,nvert,299 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 300 300 . 32, "ins(X)", zsto,zout) 301 301 c 302 302 c CALL histdef(nid_ins, "fluxec", "Cin. net flux","W/m2", 303 c . iim,jj_nb,nhori, klev,1,klev,nvert,303 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 304 304 c . 32, "ins(X)", zsto,zout) 305 305 c … … 310 310 c 311 311 CALL histdef(nid_ins, "dtdyn", "Dynamics dT", "K/s", 312 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,312 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 313 313 . "ins(X)", zsto,zout) 314 314 c 315 315 c CALL histdef(nid_ins, "dtphy", "Physics dT", "K/s", 316 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,316 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 317 317 c . "ins(X)", zsto,zout) 318 318 c 319 319 CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s", 320 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,320 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 321 321 . "ins(X)", zsto,zout) 322 322 c 323 323 CALL histdef(nid_ins, "dtajs", "Dry adjust. dT", "K/s", 324 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,324 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 325 325 . "ins(X)", zsto,zout) 326 326 c 327 327 CALL histdef(nid_ins, "dtswr", "SW radiation dT", "K/s", 328 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,328 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 329 329 . "ins(X)", zsto,zout) 330 330 c 331 331 CALL histdef(nid_ins, "dtlwr", "LW radiation dT", "K/s", 332 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,332 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 333 333 . "ins(X)", zsto,zout) 334 334 c 335 335 c CALL histdef(nid_ins, "dtec", "Cinetic dissip dT", "K/s", 336 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,336 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 337 337 c . "ins(X)", zsto,zout) 338 338 c 339 339 CALL histdef(nid_ins, "duajs", "Dry convection dU", "m/s2", 340 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,340 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 341 341 . "ins(X)", zsto,zout) 342 342 c 343 343 CALL histdef(nid_ins, "dugwo", "GW oro dU", "m/s2", 344 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,344 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 345 345 . "ins(X)", zsto,zout) 346 346 c 347 347 CALL histdef(nid_ins, "dugwno", "GW non-oro dU", "m/s2", 348 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,348 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 349 349 . "ins(X)", zsto,zout) 350 350 c 351 351 c CALL histdef(nid_ins, "dvvdf", "Boundary-layer dV", "m/s2", 352 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,352 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 353 353 c . "ins(X)", zsto,zout) 354 354 c … … 360 360 c call histdef(nid_ins, "taux", 361 361 c $ "Zonal wind stress", "Pa", 362 c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,362 c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 363 363 c $ "ins(X)", zsto,zout) 364 364 c 365 365 c call histdef(nid_ins, "tauy", 366 366 c $ "Meridional xind stress", "Pa", 367 c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,367 c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 368 368 c $ "ins(X)", zsto,zout) 369 369 c 370 370 c CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-", 371 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,371 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 372 372 c . "ins(X)", zsto,zout) 373 373 c 374 374 c CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-", 375 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,375 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 376 376 c . "ins(X)", zsto,zout) 377 377 c -
trunk/LMDZ.VENUS/libf/phyvenus/ini_histmth.h
r1518 r1530 23 23 c 24 24 CALL histdef(nid_mth, "phis", "Surface geop. height", "-", 25 . iim,jj_nb,nhori, 1,1,1, nvert, 32,25 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 26 26 . "once", zsto,zout) 27 27 c 28 28 CALL histdef(nid_mth, "aire", "Grid area", "-", 29 . iim,jj_nb,nhori, 1,1,1, nvert, 32,29 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 30 30 . "once", zsto,zout) 31 31 c 32 32 CALL histdef(nid_mth, "tsol", "Surface Temperature", "K", 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 c 36 36 CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa", 37 . iim,jj_nb,nhori, 1,1,1, nvert, 32,37 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 38 38 . "ave(X)", zsto,zout) 39 39 c 40 40 c CALL histdef(nid_mth, "ue", "Zonal energy transport", "-", 41 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,41 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 42 42 c . "ave(X)", zsto,zout) 43 43 c 44 44 c CALL histdef(nid_mth, "ve", "Merid energy transport", "-", 45 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,45 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 46 46 c . "ave(X)", zsto,zout) 47 47 c 48 48 c CALL histdef(nid_mth, "cdragh", "Drag coef on T", "-", 49 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,49 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 50 50 c . "ave(X)", zsto,zout) 51 51 c 52 52 c CALL histdef(nid_mth, "cdragm", "Drag coef on U", "-", 53 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,53 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 54 54 c . "ave(X)", zsto,zout) 55 55 c … … 62 62 c 63 63 CALL histdef(nid_mth, "temp", "Air temperature", "K", 64 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,64 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 65 65 . "ave(X)", zsto,zout) 66 66 c 67 67 CALL histdef(nid_mth, "pres", "Air pressure", "Pa", 68 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,68 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 69 69 . "ave(X)", zsto,zout) 70 70 c 71 71 CALL histdef(nid_mth, "geop", "Geopotential height", "m", 72 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,72 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 73 73 . "ave(X)", zsto,zout) 74 74 c 75 75 CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s", 76 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,76 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 77 77 . "ave(X)", zsto,zout) 78 78 c 79 79 CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s", 80 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,80 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 81 81 . "ave(X)", zsto,zout) 82 82 c 83 83 CALL histdef(nid_mth, "vitw", "Vertical wind", "Pa/s", 84 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,84 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 85 85 . "ave(X)", zsto,zout) 86 86 c 87 87 CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2", 88 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,88 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 89 89 . "ave(X)", zsto,zout) 90 90 c 91 91 CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2", 92 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,92 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 93 93 . "ave(X)", zsto,zout) 94 94 c 95 95 c CALL histdef(nid_mth, "mang", "Angular momentum", "kg m2/s", 96 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,96 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 97 97 c . "ave(X)", zsto,zout) 98 98 c 99 99 c CALL histdef(nid_mth, "Kz", "vertical diffusion coef", "m2/s", 100 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,100 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 101 101 c . "ave(X)", zsto,zout) 102 102 c 103 103 CALL histdef(nid_mth, "mmean", "Mean molecular mass", "g/mol", 104 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,104 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 105 105 . "ave(X)", zsto,zout) 106 106 107 107 c CALL histdef(nid_mth, "rho", "Air density [mass/Vol]", "kg/m3", 108 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,108 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 109 109 c . "ave(X)", zsto,zout) 110 110 … … 116 116 WRITE(str2,'(i2.2)') iq 117 117 CALL histdef(nid_mth, tname(iq), ttext(iq), "mol/mol", 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 ELSE … … 127 127 IF (callthermos .and. ok_chem) THEN 128 128 CALL histdef(nid_mth, "d_qmoldifCO2", "Dif molec" , "kg/kg", 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_mth, "d_qmoldifO3p", "Dif molec" , "kg/kg", 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 CALL histdef(nid_mth, "d_qmoldifN2", "Dif molec" , "kg/kg", 135 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,135 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 136 136 . "ave(X)", zsto,zout) 137 137 ENDIF 138 138 c 139 139 CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2", 140 . iim,jj_nb,nhori, 1,1,1, nvert, 32,140 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 141 141 . "ave(X)", zsto,zout) 142 142 c … … 151 151 c 152 152 CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2", 153 . iim,jj_nb,nhori, 1,1,1, nvert, 32,153 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 154 154 . "ave(X)", zsto,zout) 155 155 c 156 156 CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2", 157 . iim,jj_nb,nhori, 1,1,1, nvert, 32,157 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 158 158 . "ave(X)", zsto,zout) 159 159 c 160 160 CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2", 161 . iim,jj_nb,nhori, 1,1,1, nvert, 32,161 . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 162 162 . "ave(X)", zsto,zout) 163 163 c … … 165 165 c 166 166 CALL histdef(nid_mth, "SWnet", "Net SW flux","W/m2", 167 . iim,jj_nb,nhori, klev,1,klev,nvert,167 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 168 168 . 32, "ave(X)", zsto,zout) 169 169 c 170 170 CALL histdef(nid_mth, "LWnet", "Net LW flux","W/m2", 171 . iim,jj_nb,nhori, klev,1,klev,nvert,171 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 172 172 . 32, "ave(X)", zsto,zout) 173 173 c 174 174 c CALL histdef(nid_mth, "fluxvdf", "PBL net flux","W/m2", 175 c . iim,jj_nb,nhori, klev,1,klev,nvert,175 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 176 176 c . 32, "ave(X)", zsto,zout) 177 177 c 178 178 c CALL histdef(nid_mth, "fluxdyn", "Dyn. net flux","W/m2", 179 c . iim,jj_nb,nhori, klev,1,klev,nvert,179 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 180 180 c . 32, "ave(X)", zsto,zout) 181 181 c 182 182 c CALL histdef(nid_mth, "fluxajs", "Dry adj. net flux","W/m2", 183 c . iim,jj_nb,nhori, klev,1,klev,nvert,183 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 184 184 c . 32, "ave(X)", zsto,zout) 185 185 c 186 186 c CALL histdef(nid_mth, "fluxec", "Cin. net flux","W/m2", 187 c . iim,jj_nb,nhori, klev,1,klev,nvert,187 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 188 188 c . 32, "ave(X)", zsto,zout) 189 189 c … … 194 194 c 195 195 CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s", 196 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,196 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 197 197 . "ave(X)", zsto,zout) 198 198 c 199 199 c CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s", 200 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,200 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 201 201 c . "ave(X)", zsto,zout) 202 202 c 203 203 CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s", 204 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,204 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 205 205 . "ave(X)", zsto,zout) 206 206 c 207 207 CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s", 208 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,208 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 209 209 . "ave(X)", zsto,zout) 210 210 c 211 211 CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s", 212 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,212 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 213 213 . "ave(X)", zsto,zout) 214 214 c 215 215 CALL histdef(nid_mth,"dtswrNLTE", "SWNLTE radiation dT", 216 . "K/s", iim,jj_nb,nhori, klev,1,klev,nvert,32,216 . "K/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 217 217 . "ave(X)", zsto,zout) 218 218 c 219 219 CALL histdef(nid_mth,"dtswrDCrisp","SWDCrisp radiation dT", 220 . "K/s", iim,jj_nb,nhori, klev,1,klev,nvert,32,220 . "K/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 221 221 . "ave(X)", zsto,zout) 222 222 c 223 223 CALL histdef(nid_mth, "dtlwr", "LW radiation 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 . "ave(X)", zsto,zout) 226 226 c 227 227 CALL histdef(nid_mth,"dtlwrNLTE", "LWNLTE radiation dT", 228 . "K/s", iim,jj_nb,nhori, klev,1,klev,nvert,32,228 . "K/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 229 229 . "ave(X)", zsto,zout) 230 230 c 231 231 CALL histdef(nid_mth,"dtlwrLTE","LW_LTE radiation dT", 232 . "K/s", iim,jj_nb,nhori, klev,1,klev,nvert,32,232 . "K/s",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 233 233 . "ave(X)", zsto,zout) 234 234 c 235 235 CALL histdef(nid_mth, "dteuv", "UV 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 . "ave(X)", zsto,zout) 238 238 CALL histdef(nid_mth, "dtcond", "Therm conduction", "K/s", 239 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,239 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 240 240 . "ave(X)", zsto,zout) 241 241 CALL histdef(nid_mth, "dumolvis", "molec viscosity (u)" 242 . ,"m/s2", iim,jj_nb,nhori, klev,1,klev,nvert,32,242 . ,"m/s2",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 243 243 . "ave(X)", zsto,zout) 244 244 CALL histdef(nid_mth, "dvmolvis", "molec viscosity (v)", 245 . "m/s2", iim,jj_nb,nhori, klev,1,klev,nvert,32,245 . "m/s2",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 246 246 . "ave(X)", zsto,zout) 247 247 248 248 c 249 249 c CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s", 250 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,250 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 251 251 c . "ave(X)", zsto,zout) 252 252 c 253 253 CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2", 254 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,254 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 255 255 . "ave(X)", zsto,zout) 256 256 c 257 257 CALL histdef(nid_mth, "dugwo", "GW oro dU", "m/s2", 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)", zsto,zout) 260 260 c 261 261 CALL histdef(nid_mth, "dugwno", "GW non-oro dU", "m/s2", 262 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,262 . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 263 263 . "ave(X)", zsto,zout) 264 264 265 265 c CALL histdef(nid_mth, "co2_vmr", "volume mixture ratio co2", 266 c . "cm-3", iim,jj_nb,nhori, klev,1,klev,nvert, 32,266 c . "cm-3",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 267 267 c . "ave(X)", zsto,zout) 268 268 269 269 c CALL histdef(nid_mth, "o_vmr", "volume mixture ratio o", 270 c . "cm-3", iim,jj_nb,nhori, klev,1,klev,nvert,32,270 c . "cm-3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 271 271 c . "ave(X)", zsto,zout) 272 272 c CALL histdef(nid_mth, "co_vmr", "volume mixture ratio co", 273 c . "cm-3", iim,jj_nb,nhori, klev,1,klev,nvert, 32,273 c . "cm-3",nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 274 274 c . "ave(X)", zsto,zout) 275 275 276 276 c CALL histdef(nid_mth, "n2_vmr", "volume mixture ratio n2", 277 c . "cm-3", iim,jj_nb,nhori, klev,1,klev,nvert,32,277 c . "cm-3",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32, 278 278 c . "ave(X)", zsto,zout) 279 279 c 280 280 c CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2", 281 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32,281 c . nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32, 282 282 c . "ave(X)", zsto,zout) 283 283 c … … 289 289 c call histdef(nid_mth, "taux", 290 290 c $ "Zonal wind stress", "Pa", 291 c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,291 c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 292 292 c $ "ave(X)", zsto,zout) 293 293 c 294 294 c call histdef(nid_mth, "tauy", 295 295 c $ "Meridional xind stress", "Pa", 296 c $ iim,jj_nb,nhori, 1,1,1, nvert, 32,296 c $ nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 297 297 c $ "ave(X)", zsto,zout) 298 298 c 299 299 c CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-", 300 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,300 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 301 301 c . "ave(X)", zsto,zout) 302 302 c 303 303 c CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-", 304 c . iim,jj_nb,nhori, 1,1,1, nvert, 32,304 c . nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32, 305 305 c . "ave(X)", zsto,zout) 306 306 c -
trunk/LMDZ.VENUS/libf/phyvenus/iophy.F90
r892 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 ! POUR VENUS, LA CARTE EST A L ENVERS !! … … 68 73 69 74 ddid=(/ 1,2 /) 70 dsg=(/ iim, jjm+1-1/(iim*jjm)/)71 dsl=(/ iim, jj_nb /)75 dsg=(/ nbp_lon, nbp_lat /) 76 dsl=(/ nbp_lon, jj_nb /) 72 77 dpf=(/ 1,jj_begin /) 73 dpl=(/ iim, jj_end /)78 dpl=(/ nbp_lon, jj_end /) 74 79 dhs=(/ ii_begin-1,0 /) 75 if (mpi_rank==mpi_size-1) then80 IF (mpi_rank==mpi_size-1) THEN 76 81 dhe=(/0,0/) 77 else78 dhe=(/ iim-ii_end,0 /)79 endif82 ELSE 83 dhe=(/ nbp_lon-ii_end,0 /) 84 ENDIF 80 85 81 86 call flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe, & … … 89 94 USE dimphy 90 95 USE mod_phys_lmdz_para 91 use ioipsl 92 implicit none93 i nclude 'dimensions.h'94 real,dimension( iim),intent(in) :: lon95 real,dimension( jjm+1-1/(iim*jjm)),intent(in) :: lat96 use ioipsl, only: flio_dom_set 97 USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat 98 implicit none 99 real,dimension(nbp_lon),intent(in) :: lon 100 real,dimension(nbp_lat),intent(in) :: lat 96 101 97 102 INTEGER,DIMENSION(2) :: ddid … … 104 109 105 110 !$OMP MASTER 106 allocate(io_lat( jjm+1-1/(iim*jjm)))111 allocate(io_lat(nbp_lat)) 107 112 io_lat(:)=lat(:) 108 allocate(io_lon( iim))113 allocate(io_lon(nbp_lon)) 109 114 io_lon(:)=lon(:) 110 115 111 116 ddid=(/ 1,2 /) 112 dsg=(/ iim, jjm+1-1/(iim*jjm)/)113 dsl=(/ iim, jj_nb /)117 dsg=(/ nbp_lon, nbp_lat /) 118 dsl=(/ nbp_lon, jj_nb /) 114 119 dpf=(/ 1,jj_begin /) 115 dpl=(/ iim, jj_end /)120 dpl=(/ nbp_lon, jj_end /) 116 121 dhs=(/ ii_begin-1,0 /) 117 122 if (mpi_rank==mpi_size-1) then 118 123 dhe=(/0,0/) 119 124 else 120 dhe=(/ iim-ii_end,0 /)125 dhe=(/ nbp_lon-ii_end,0 /) 121 126 endif 122 127 … … 129 134 130 135 subroutine histbeg_phy(name,itau0,zjulian,dtime,nhori,nid_day) 131 USE dimphy132 USE mod_ phys_lmdz_para133 use ioipsl 136 USE mod_phys_lmdz_para, only: jj_begin, jj_end, jj_nb, is_sequential 137 USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat 138 use ioipsl, only: histbeg 134 139 use write_field 135 140 implicit none 136 include 'dimensions.h'137 141 138 142 character*(*), intent(IN) :: name … … 145 149 !$OMP MASTER 146 150 if (is_sequential) then 147 call histbeg(name, iim,io_lon, jj_nb,io_lat(jj_begin:jj_end), &148 1, iim,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day)151 call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), & 152 1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day) 149 153 else 150 call histbeg(name, iim,io_lon, jj_nb,io_lat(jj_begin:jj_end), &151 1, iim,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day,phys_domain_id)154 call histbeg(name,nbp_lon,io_lon, jj_nb,io_lat(jj_begin:jj_end), & 155 1,nbp_lon,1,jj_nb,itau0, zjulian, dtime, nhori, nid_day,phys_domain_id) 152 156 endif 153 157 !$OMP END MASTER … … 158 162 plon,plat,plon_bounds,plat_bounds, & 159 163 nname,itau0,zjulian,dtime,nnhori,nnid_day) 160 USE dimphy 164 USE dimphy, only: klon 161 165 USE mod_phys_lmdz_para 162 USE mod_grid_phy_lmdz 163 use ioipsl 164 use write_field 165 implicit none 166 include 'dimensions.h' 166 USE mod_grid_phy_lmdz, only: klon_glo, nbp_lon, nbp_lat 167 use ioipsl, only: histbeg 168 169 implicit none 167 170 168 171 real,dimension(klon),intent(in) :: rlon … … 190 193 REAL, allocatable, dimension(:) :: npplat, npplon 191 194 REAL, allocatable, dimension(:,:) :: npplat_bounds, npplon_bounds 192 INTEGER, PARAMETER :: jjmp1=jjm+1-1/jjm 193 REAL, dimension(iim,jjmp1) :: zx_lon, zx_lat 195 REAL, dimension(nbp_lon,nbp_lat) :: zx_lon, zx_lat 194 196 195 197 CALL gather(rlat,rlat_glo) … … 216 218 endif 217 219 ! 218 IF ( tabij(i).LE. iim) THEN220 IF ( tabij(i).LE.nbp_lon) THEN 219 221 plat_bounds(i,1)=rlat_glo(tabij(i)) 220 222 ELSE 221 plat_bounds(i,1)=rlat_glo(tabij(i)- iim)223 plat_bounds(i,1)=rlat_glo(tabij(i)-nbp_lon) 222 224 ENDIF 223 plat_bounds(i,2)=rlat_glo(tabij(i)+ iim)225 plat_bounds(i,2)=rlat_glo(tabij(i)+nbp_lon) 224 226 ! 225 227 ! print*,'CFMIP_iophy point i lon lon_bds',i,plon_bounds(i,1),rlon_glo(tabij(i)),plon_bounds(i,2) … … 237 239 ENDDO 238 240 239 CALL gr_fi_ecrit(1,klon, iim,jjmp1,rlon_glo,zx_lon)240 if (( iim*jjm).gt.1) then241 DO i = 1, iim241 CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlon_glo,zx_lon) 242 if ((nbp_lon*nbp_lat).gt.1) then 243 DO i = 1, nbp_lon 242 244 zx_lon(i,1) = rlon_glo(i+1) 243 zx_lon(i, jjmp1) = rlon_glo(i+1)245 zx_lon(i,nbp_lat) = rlon_glo(i+1) 244 246 ENDDO 245 247 endif 246 CALL gr_fi_ecrit(1,klon, iim,jjmp1,rlat_glo,zx_lat)248 CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlat_glo,zx_lat) 247 249 248 250 DO i=1,pim … … 253 255 254 256 if (ipt(i).EQ.1) then 255 plon_bounds(i,1)=zx_lon( iim,jpt(i))257 plon_bounds(i,1)=zx_lon(nbp_lon,jpt(i)) 256 258 plon_bounds(i,2)=360.+zx_lon(ipt(i)+1,jpt(i)) 257 259 endif 258 260 259 if (ipt(i).EQ. iim) then261 if (ipt(i).EQ.nbp_lon) then 260 262 plon_bounds(i,2)=360.+zx_lon(1,jpt(i)) 261 263 endif … … 269 271 endif 270 272 271 if (jpt(i).EQ. jjmp1) then272 plat_bounds(i,1)=zx_lat(ipt(i), jjmp1)+0.001273 plat_bounds(i,2)=zx_lat(ipt(i), jjmp1)-0.001273 if (jpt(i).EQ.nbp_lat) then 274 plat_bounds(i,1)=zx_lat(ipt(i),nbp_lat)+0.001 275 plat_bounds(i,2)=zx_lat(ipt(i),nbp_lat)-0.001 274 276 endif 275 277 ! … … 327 329 328 330 subroutine histwrite2d_phy(nid,lpoint,name,itau,field) 329 USE dimphy 331 USE dimphy, only: klon 330 332 USE mod_phys_lmdz_para 331 USE ioipsl 332 implicit none333 i nclude 'dimensions.h'333 USE ioipsl, only: histwrite 334 USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat 335 implicit none 334 336 335 337 integer,intent(in) :: nid … … 340 342 REAL,dimension(klon_mpi) :: buffer_omp 341 343 INTEGER, allocatable, dimension(:) :: index2d 342 REAL :: Field2d( iim,jj_nb)344 REAL :: Field2d(nbp_lon,jj_nb) 343 345 344 346 integer :: ip … … 351 353 CALL grid1Dto2D_mpi(buffer_omp,Field2d) 352 354 if(.NOT.lpoint) THEN 353 ALLOCATE(index2d( iim*jj_nb))354 ALLOCATE(fieldok( iim*jj_nb))355 CALL histwrite(nid,name,itau,Field2d, iim*jj_nb,index2d)355 ALLOCATE(index2d(nbp_lon*jj_nb)) 356 ALLOCATE(fieldok(nbp_lon*jj_nb)) 357 CALL histwrite(nid,name,itau,Field2d,nbp_lon*jj_nb,index2d) 356 358 else 357 359 ALLOCATE(fieldok(npstn)) … … 382 384 383 385 subroutine histwrite3d_phy(nid,lpoint,name,itau,field) 384 USE dimphy 386 USE dimphy, only: klon 385 387 USE mod_phys_lmdz_para 386 387 use ioipsl 388 implicit none 389 include 'dimensions.h' 388 USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat 389 use ioipsl, only: histwrite 390 implicit none 390 391 391 392 integer,intent(in) :: nid … … 395 396 real,dimension(:,:),intent(in) :: field ! --> field(klon,:) 396 397 REAL,dimension(klon_mpi,size(field,2)) :: buffer_omp 397 REAL :: Field3d( iim,jj_nb,size(field,2))398 REAL :: Field3d(nbp_lon,jj_nb,size(field,2)) 398 399 INTEGER :: ip, n, nlev 399 400 INTEGER, ALLOCATABLE, dimension(:) :: index3d … … 413 414 CALL grid1Dto2D_mpi(buffer_omp,field3d) 414 415 if(.NOT.lpoint) THEN 415 ALLOCATE(index3d( iim*jj_nb*nlev))416 ALLOCATE(fieldok( iim*jj_nb,nlev))417 CALL histwrite(nid,name,itau,Field3d, iim*jj_nb*nlev,index3d)416 ALLOCATE(index3d(nbp_lon*jj_nb*nlev)) 417 ALLOCATE(fieldok(nbp_lon*jj_nb,nlev)) 418 CALL histwrite(nid,name,itau,Field3d,nbp_lon*jj_nb*nlev,index3d) 418 419 else 419 420 nlev=size(field,2) -
trunk/LMDZ.VENUS/libf/phyvenus/jthermcalc.F
r1442 r1530 18 18 19 19 c common variables and constants 20 include "dimensions.h"21 20 include "param.h" 22 21 include "param_v4.h" … … 1062 1061 1063 1062 c common variables and constants 1064 #include "dimensions.h"1065 1063 c#include "tracer.h" 1066 1064 #include "param.h" … … 1408 1406 1409 1407 c common variables and constants 1410 #include "dimensions.h"1411 1408 #include "param.h" 1412 1409 #include "param_v4.h" … … 1702 1699 1703 1700 ! common variables and constants 1704 include "dimensions.h"1705 1701 include "param.h" 1706 1702 include 'param_v4.h' -
trunk/LMDZ.VENUS/libf/phyvenus/jthermcalc_e107.F
r1310 r1530 16 16 17 17 c common variables and constants 18 #include "dimensions.h"19 18 #include "param.h" 20 19 #include "param_v4.h" -
trunk/LMDZ.VENUS/libf/phyvenus/load_ksi.F
r1301 r1530 4 4 IMPLICIT none 5 5 6 #include "dimensions.h"7 6 #include "YOMCST.h" 8 7 #include "comcstVE.h" -
trunk/LMDZ.VENUS/libf/phyvenus/lw_venus_ve.F
r1301 r1530 9 9 IMPLICIT none 10 10 11 #include "dimensions.h"12 11 #include "YOMCST.h" 13 12 C -
trunk/LMDZ.VENUS/libf/phyvenus/lwi.F
r1017 r1530 6 6 7 7 8 #include "dimensions.h"9 8 #include "YOMCST.h" 10 9 #include "timerad.h" -
trunk/LMDZ.VENUS/libf/phyvenus/moldiff_red.F90
r1442 r1530 8 8 implicit none 9 9 10 !#include "dimensions.h"11 10 !#include "dimphys.h" 12 11 #include "comcstfi.h" … … 965 964 use infotrac 966 965 IMPLICIT NONE 967 !#include "dimensions.h"968 966 969 967 INTEGER :: nl,nq,l, nn … … 1051 1049 use infotrac 1052 1050 IMPLICIT NONE 1053 !#include "dimensions.h"1054 1051 1055 1052 INTEGER :: nl,nq,il,l,i,iq,nlx,iz,ig … … 1394 1391 use infotrac 1395 1392 IMPLICIT NONE 1396 !#include "dimensions.h"1397 1393 INTEGER :: nl,nq,nlx,il,nn,iP,ig,compteur 1398 1394 INTEGER,DIMENSION(1) :: indP … … 1497 1493 ! use infotrac 1498 1494 IMPLICIT NONE 1499 !#include "dimensions.h"1500 1495 INTEGER :: nl,nq,nlx,il,nn,iP,ig,compteur 1501 1496 INTEGER,DIMENSION(1) :: indP -
trunk/LMDZ.VENUS/libf/phyvenus/moldiffcoeff_red.F
r1442 r1530 16 16 c 17 17 c======================================================================= 18 !#include "dimensions.h"19 !#include "dimphys.h"20 !#include "callkeys.h"21 !#include "comdiurn.h"22 !#include "chimiedata.h"23 !#include "tracer.h"24 !#include "conc.h"25 18 #include "diffusion.h" 26 19 -
trunk/LMDZ.VENUS/libf/phyvenus/molvis.F
r1310 r1530 20 20 c declarations: 21 21 c----------------------------------------------------------------------- 22 23 !#include "dimensions.h"24 !#include "dimphys.h"25 !#include "comcstfi.h"26 !#include "surfdat.h"27 !#include "chimiedata.h"28 !#include "conc.h"29 22 30 23 c arguments: -
trunk/LMDZ.VENUS/libf/phyvenus/new_cloud_sedim.F
r1525 r1530 19 19 c ------------- 20 20 #include "YOMCST.h" 21 c#include "dimensions.h"22 21 c#include "dimphys.h" 23 22 c#include "comcstfi.h" -
trunk/LMDZ.VENUS/libf/phyvenus/nirco2abs.F
r1442 r1530 54 54 c 55 55 56 !#include "dimensions.h"57 56 #include "YOMCST.h" 58 57 #include "clesphys.h" -
trunk/LMDZ.VENUS/libf/phyvenus/nlte_tcool.F
r1442 r1530 33 33 implicit none 34 34 35 include "dimensions.h"36 35 include "nlte_paramdef.h" 37 36 include "nlte_commons.h" -
trunk/LMDZ.VENUS/libf/phyvenus/nltecool.F
r1310 r1530 39 39 40 40 #include "nltedata.h" ! (Equivalent to the reading of the "nlte_escape.dat" file) 41 #include "dimensions.h"42 43 !#include "dimphys.h"44 !#include "chimiedata.h"45 !#include "conc.h" !Added to have "dynamic composition" in the scheme46 !#include "tracer.h" !"47 !#include "callkeys.h"48 41 49 42 c Input and output variables -
trunk/LMDZ.VENUS/libf/phyvenus/nlthermeq.F
r1310 r1530 10 10 use dimphy 11 11 implicit none 12 #include "dimensions.h"13 12 c#include "dimradmars.h" 14 13 #include "nlteparams.h" -
trunk/LMDZ.VENUS/libf/phyvenus/orodrag.F
r101 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.VENUS/libf/phyvenus/orosetup.F
r101 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.VENUS/libf/phyvenus/param_read.F
r1452 r1530 7 7 8 8 c common variables and constants 9 include "dimensions.h"10 9 include 'param.h' 11 10 include 'param_v4.h' -
trunk/LMDZ.VENUS/libf/phyvenus/param_read_e107.F
r1310 r1530 7 7 8 8 c common variables and constants 9 include "dimensions.h"10 9 include 'param.h' 11 10 include 'param_v4.h' -
trunk/LMDZ.VENUS/libf/phyvenus/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.VENUS/libf/phyvenus/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.VENUS/libf/phyvenus/physiq.F
r1525 r1530 70 70 ! use ieee_arithmetic 71 71 use time_phylmdz_mod, only: annee_ref, day_ref, itau_phy 72 use mod_grid_phy_lmdz, only: nbp_lon 72 73 use logic_mod, only: iflag_trac 73 74 IMPLICIT none … … 80 81 #define histins 81 82 c====================================================================== 82 #include "dimensions.h"83 integer jjmp184 parameter (jjmp1=jjm+1-1/jjm)85 83 #include "dimsoil.h" 86 84 #include "clesphys.h" -
trunk/LMDZ.VENUS/libf/phyvenus/phytrac_chimie.F
r1525 r1530 28 28 IMPLICIT none 29 29 30 c#include "dimensions.h"31 30 #include "clesphys.h" 32 31 c#include "temps.h" -
trunk/LMDZ.VENUS/libf/phyvenus/phytrac_emiss.F
r1525 r1530 32 32 USE comgeomphy 33 33 USE chemparam_mod,only:M_tr 34 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 34 35 IMPLICIT none 35 36 #include "YOMCST.h" 36 #include "dimensions.h"37 37 #include "clesphys.h" 38 #include "paramet.h"39 38 c====================================================================== 40 39 … … 124 123 125 124 ig_volcan(ilat,ilon)= 0 126 deltalat = 180./jjm 127 deltalon = 360./iim 125 if ((nbp_lon*nbp_lat)==1) then ! running a 1D simulation 126 deltalat=180. 127 deltalon=360. 128 else 129 deltalat = 180./(nbp_lat-1) 130 deltalon = 360./nbp_lon 131 endif 128 132 129 133 do i=1,nlon -
trunk/LMDZ.VENUS/libf/phyvenus/phytrac_relax.F
r1525 r1530 31 31 IMPLICIT none 32 32 #include "YOMCST.h" 33 #include "dimensions.h"34 33 #include "clesphys.h" 35 #include "paramet.h"36 34 c====================================================================== 37 35 -
trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.F
r1442 r1530 32 32 use write_field_phy 33 33 IMPLICIT none 34 #include "dimensions.h"35 34 #include "YOMCST.h" 36 35 #include "clesphys.h" -
trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.NewtonCool
r1310 r1530 29 29 USE write_field_phy 30 30 IMPLICIT none 31 #include "dimensions.h"32 31 #include "YOMCST.h" 33 32 #include "clesphys.h" -
trunk/LMDZ.VENUS/libf/phyvenus/soil.F
r101 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.VENUS/libf/phyvenus/sugwd.F
r101 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.VENUS/libf/phyvenus/sw_venus_cl.F
r1442 r1530 8 8 IMPLICIT none 9 9 10 #include "dimensions.h"11 10 #include "YOMCST.h" 12 11 C -
trunk/LMDZ.VENUS/libf/phyvenus/sw_venus_cl_1Dglobave.F
r1442 r1530 8 8 IMPLICIT none 9 9 10 #include "dimensions.h"11 10 #include "YOMCST.h" 12 11 C -
trunk/LMDZ.VENUS/libf/phyvenus/sw_venus_dc.F
r1442 r1530 8 8 IMPLICIT none 9 9 10 #include "dimensions.h"11 10 #include "YOMCST.h" 12 11 C -
trunk/LMDZ.VENUS/libf/phyvenus/sw_venus_dc_1Dglobave.F
r1442 r1530 8 8 IMPLICIT none 9 9 10 #include "dimensions.h"11 10 #include "YOMCST.h" 12 11 C -
trunk/LMDZ.VENUS/libf/phyvenus/sw_venus_ve.F
r1442 r1530 8 8 IMPLICIT none 9 9 10 #include "dimensions.h"11 10 #include "YOMCST.h" 12 11 C -
trunk/LMDZ.VENUS/libf/phyvenus/sw_venus_ve_1Dglobave.F
r1442 r1530 8 8 IMPLICIT none 9 9 10 #include "dimensions.h"11 10 #include "YOMCST.h" 12 11 C -
trunk/LMDZ.VENUS/libf/phyvenus/transp.F
r101 r1530 13 13 use dimphy 14 14 IMPLICIT none 15 #include "dimensions.h"16 15 #include "YOMCST.h" 17 16 c -
trunk/LMDZ.VENUS/libf/phyvenus/ustarhb.F
r892 r1530 18 18 use dimphy 19 19 IMPLICIT none 20 #include "dimensions.h"21 20 #include "YOMCST.h" 22 21 c -
trunk/LMDZ.VENUS/libf/phyvenus/vdif_kcay.F
r101 r1530 8 8 use dimphy 9 9 IMPLICIT NONE 10 #include "dimensions.h"11 10 c....................................................................... 12 11 c -
trunk/LMDZ.VENUS/libf/phyvenus/yamada.F
r101 r1530 8 8 use dimphy 9 9 IMPLICIT NONE 10 #include "dimensions.h"11 10 c....................................................................... 12 11 c -
trunk/LMDZ.VENUS/libf/phyvenus/yamada4.F
r101 r1530 9 9 use dimphy 10 10 IMPLICIT NONE 11 #include "dimensions.h"12 11 c....................................................................... 13 12 c -
trunk/LMDZ.VENUS/libf/phyvenus/zenang.F
r1301 r1530 26 26 use dimphy 27 27 IMPLICIT none 28 #include "dimensions.h"29 28 #include "YOMCST.h" 30 29 c================================================================
Note: See TracChangeset
for help on using the changeset viewer.
