- Timestamp:
- Jul 20, 2024, 12:00:23 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/etat0dyn_netcdf.F90
r4050 r5087 122 122 CALL startget_dyn3d('u' ,rlonu,rlatu,pls,y ,uvent,rlonv,rlatv) 123 123 CALL startget_dyn3d('v' ,rlonv,rlatv,pls(:,:jjm,:),y(:,:jjm,:),vvent, & 124 &rlonu,rlatu(:jjm))124 rlonu,rlatu(:jjm)) 125 125 CALL startget_dyn3d('t' ,rlonv,rlatu,pls,y ,t3d ,rlonu,rlatv) 126 126 tpot(:,:,:)=t3d(:,:,:) … … 228 228 !------------------------------------------------------------------------------- 229 229 iml=assert_eq([SIZE(lon_in),SIZE(pls,1),SIZE(workvar,1),SIZE(champ,1), & 230 &SIZE(lon_in2)], TRIM(modname)//" iml")230 SIZE(lon_in2)], TRIM(modname)//" iml") 231 231 jml=assert_eq( SIZE(lat_in),SIZE(pls,2),SIZE(workvar,2),SIZE(champ,2), & 232 &TRIM(modname)//" jml")232 TRIM(modname)//" jml") 233 233 lml=assert_eq( SIZE(pls,3),SIZE(workvar,3),SIZE(champ,3), & 234 &TRIM(modname)//" lml")234 TRIM(modname)//" lml") 235 235 jml2=SIZE(lat_in2) 236 236 … … 306 306 !------------------------------------------------------------------------------- 307 307 iml=assert_eq(SIZE(lon_in),SIZE(zs,1),SIZE(psol,1),SIZE(lon_in2), & 308 &TRIM(modname)//" iml")308 TRIM(modname)//" iml") 309 309 jml=assert_eq(SIZE(lat_in),SIZE(zs,2),SIZE(psol,2),TRIM(modname)//" jml") 310 310 jml2=SIZE(lat_in2) … … 333 333 !--- PSOL IS COMPUTED IN PASCALS 334 334 psol(:iml-1,:) = ps(:iml-1,:)*(1.0+(z(:iml-1,:)-zs(:iml-1,:))/287.0 & 335 &/ts(:iml-1,:))335 /ts(:iml-1,:)) 336 336 psol(iml,:)=psol(1,:) 337 337 DEALLOCATE(z,ps,ts)
Note: See TracChangeset
for help on using the changeset viewer.