Changeset 1756 for trunk/MESOSCALE/LMD_MM_MARS/SRC
- Timestamp:
- Jul 25, 2017, 11:38:07 AM (7 years ago)
- Location:
- trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/Venus
- Files:
-
- 1 added
- 2 edited
-
README (added)
-
create_readmeteo.F90 (modified) (13 diffs)
-
readmeteo.F90 (modified) (35 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/Venus/create_readmeteo.F90
r1733 r1756 12 12 !------------------------------------------------------------------------! 13 13 14 INTEGER, PARAMETER :: MONTHS_PER_YEAR = 12 15 INTEGER, PARAMETER :: mday(MONTHS_PER_YEAR) & 16 = (/61,66,66,65,60,54,50,46,47,47,51,56/) 17 14 INTEGER, PARAMETER :: MONTHS_PER_YEAR = 100 15 INTEGER, PARAMETER :: mday = 24 18 16 INTEGER :: start_day,init,i,month,day 19 17 INTEGER :: ierr,nid,nvarid … … 22 20 INTEGER :: no_please 23 21 INTEGER :: timedim,timelen 24 22 INTEGER :: inc 25 23 REAL, DIMENSION(100) :: param 26 24 REAL, DIMENSION(:), ALLOCATABLE :: time … … 32 30 interval = 0 33 31 34 32 !mday(:)=24 35 33 ! 36 34 ! Open input NETCDF file … … 38 36 write(*,*) "Scanning netcdf file ..." 39 37 ierr=NF_OPEN ("input_diagfi.nc",NF_NOWRITE,nid) 40 ierr=NF__OPEN ("input_diagfi.nc",NF_NOWRITE,4096,nid) 41 !ierr=nc_open("input_diagfi.nc",NC_NOWRITE,nid) 38 !ierr=NF90_OPEN("input_diagfi.nc",NF_NOWRITE,nid) 42 39 IF (ierr.NE.NF_NOERR) THEN 43 40 write(*,*)'**** Please create a symbolic link called input_diagfi.nc' … … 45 42 ENDIF 46 43 47 ierr=NF_INQ_DIMID(nid," Time",timedim)44 ierr=NF_INQ_DIMID(nid,"time_counter",timedim) 48 45 IF (ierr .NE. NF_NOERR) THEN 49 46 ierr=NF_INQ_DIMID(nid,"time",timedim) … … 56 53 ! 57 54 ALLOCATE(time(timelen)) 58 ierr = NF_INQ_VARID (nid, " Time",nvarid)55 ierr = NF_INQ_VARID (nid, "time_counter",nvarid) 59 56 IF (ierr .NE. NF_NOERR) THEN 60 57 ierr = NF_INQ_VARID (nid, "time",nvarid) … … 73 70 IF (ierr .NE. NF_NOERR) THEN 74 71 PRINT *, "Error: Readmeteo <ps> not found" 75 stop72 !stop 76 73 ENDIF 77 74 #ifdef NC_DOUBLE … … 92 89 ! beware, param(4) is the day reference of start and startfi 93 90 ! ...have to add time(1) to get the real starting date in diagfi 94 start_day=floor(param(4)+time(1)) 95 start_hour=nint((param(4)-floor(param(4))+time(1))*24) ! starting hour 96 start_hour=MOD(start_hour,24) 97 IF (interval .eq. 0) interval=nint(time(1)*24) ! interval between each time subscript 98 99 91 !start_day=floor(param(4)+time(1)) 92 !start_hour=nint((param(4)-floor(param(4))+time(1))*13) ! starting hour 93 !start_hour=MOD(start_hour,24) 94 !IF (interval .eq. 0) interval=nint(time(1)*24) ! interval between each time subscript 95 96 start_day=0 97 start_hour=0 98 interval=864 ! in hours, 864 HH = 1 MM = 1/100 Vd 100 99 PRINT *,'*****************' 101 100 PRINT *,'GCM data file starts at sol ',start_day,'and hour',start_hour 102 101 PRINT *,'GCM data interval is ',interval,'hours' 103 102 104 CALL wrf_day(start_day,month,day) 105 103 month=0 104 day=0 105 !interval=864 106 106 ! 107 107 ! User defined parameters … … 116 116 if (no_please == 0) stop 117 117 if (no_please == 1) then 118 my= 2024118 my=1000 119 119 n=timelen 120 start= 1120 start=0 121 121 interval_subs=1 122 122 else … … 124 124 write(*,*) "Starting Martian year ? ex: 24,25,26..." 125 125 read(*,*) my 126 my= 2000+my126 my=1000+my 127 127 write(*,*) "-- WRF data file information --" 128 128 write(*,*) "How many files do you want to create ? at least 2, max is",timelen … … 177 177 END DO 178 178 ! WRF time reference 179 hour= start_hour+(start-1)*interval179 hour=interval 180 180 inc_hour=interval*interval_subs 181 IF (hour >= 24) day=day+INT(hour/24) 182 hour=MOD(hour,24) 183 IF (day > mday(month)) THEN 184 day=day-mday(month) 185 month=month+1 181 IF (hour >= 36) day=day+INT(hour/36) 182 hour=MOD(hour,36) 183 IF (day >= mday) THEN 184 inc=INT(day/mday) 185 day=day-inc*mday 186 month=month+inc 186 187 END IF 187 IF (month > MONTHS_PER_YEAR) THEN188 IF (month >= MONTHS_PER_YEAR) THEN 188 189 my=my+1 189 month= 1190 month=0 190 191 END IF 191 192 DO i=1,n … … 207 208 END IF 208 209 write(6,fmt='(A1)') 'y' 209 IF (hour+inc_hour >= 24) day=day+INT((hour+inc_hour)/24) 210 hour=MOD(hour+inc_hour,24) 211 IF (day > mday(month)) THEN 212 day=day-mday(month) 213 month=month+1 210 IF (hour+inc_hour >= 36) day=day+INT((hour+inc_hour)/36) 211 hour=MOD(hour+inc_hour,36) 212 IF (day >= mday) THEN 213 inc=INT(day/mday) 214 day=day-inc*mday 215 month=month+inc 214 216 END IF 215 IF (month > MONTHS_PER_YEAR) THEN217 IF (month >= MONTHS_PER_YEAR) THEN 216 218 my=my+1 217 month= 1219 month=0 218 220 END IF 219 221 END DO … … 221 223 222 224 END 223 224 225 !--------------------------------------------------226 !--------------------------------------------------227 228 SUBROUTINE wrf_day(gcm_day,wrf_month,day)229 230 IMPLICIT NONE231 232 INTEGER, INTENT(INOUT) :: gcm_day233 INTEGER, INTENT(OUT) :: wrf_month,day234 235 INTEGER :: init,i236 INTEGER, PARAMETER :: MONTHS_PER_YEAR = 12237 INTEGER, PARAMETER :: mday(MONTHS_PER_YEAR) &238 = (/61,66,66,65,60,54,50,46,47,47,51,56/)239 240 !241 ! Find WRF month and day242 !243 244 IF (gcm_day >= 669) THEN !! gcm_day commence au jour 0245 PRINT *,'out of bounds ! martian year is 669 sols !'246 gcm_day=MOD(gcm_day,669)247 !STOP248 ENDIF249 250 251 252 init=gcm_day+1 !!+1 sinon on decale tout253 DO i=1,MONTHS_PER_YEAR254 wrf_month=i255 init=init-mday(i)256 IF (init <= 0) EXIT257 END DO258 259 PRINT *,'corresponding WRF month is ',wrf_month260 day=init+mday(wrf_month)261 PRINT *,'corresponding WRF day is ',day262 PRINT *,'*****************'263 264 END -
trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/Venus/readmeteo.F90
r1733 r1756 41 41 42 42 REAL :: ptop 43 REAL, PARAMETER :: grav= 3.7243 REAL, PARAMETER :: grav=8.89 44 44 LOGICAL, PARAMETER :: blank=.false. 45 45 … … 73 73 74 74 !! Intermediate data arrays 75 integer :: k,l,m,n,p 75 integer :: k,l,m,n,p,i,j 76 76 real, dimension(:), allocatable :: lat,lon,time,alt,aps,bps,levels,vertdsoil 77 77 real, dimension(:,:), allocatable :: vide,ones,ghtsfile … … 81 81 real, dimension(:,:,:), allocatable :: emissfile,co2icefile 82 82 real, dimension(:,:,:), allocatable :: tnsfile,unsfile,vnsfile 83 real, dimension(:,:,:,:), allocatable :: tfile,ufile,vfile,rfile,hfile 83 real, dimension(:,:,:,:), allocatable :: tfile,ufile,vfile,rfile,hfile,ghfile 84 real, dimension(:,:,:,:), allocatable :: pfile 84 85 real, dimension(:,:,:,:), allocatable :: eta_gcm 85 86 !real, dimension(:,:,:,:), allocatable :: tfileorig,ufileorig,vfileorig … … 173 174 SELECT CASE(ident) 174 175 CASE('LMD') 175 ierr=NF_INQ_DIMID(nid," Time",timedim)176 ierr=NF_INQ_DIMID(nid,"time_counter",timedim) 176 177 CASE('OXF') 177 178 ierr=NF_INQ_DIMID(nid,"time",timedim) … … 182 183 SELECT CASE(ident) 183 184 CASE('LMD') 184 ierr=NF_INQ_DIMID(nid,"lat itude",latdim)185 ierr=NF_INQ_DIMID(nid,"lat",latdim) 185 186 CASE('OXF') 186 187 ierr=NF_INQ_DIMID(nid,"lat",latdim) … … 191 192 SELECT CASE(ident) 192 193 CASE('LMD') 193 ierr=NF_INQ_DIMID(nid,"lon gitude",londim)194 ierr=NF_INQ_DIMID(nid,"lon",londim) 194 195 CASE('OXF') 195 196 ierr=NF_INQ_DIMID(nid,"lon",londim) … … 200 201 SELECT CASE(ident) 201 202 CASE('LMD') 202 ierr=NF_INQ_DIMID(nid," altitude",altdim)203 ierr=NF_INQ_DIMID(nid,"presnivs",altdim) 203 204 CASE('OXF') 204 205 ierr=NF_INQ_DIMID(nid,"sigma",altdim) … … 214 215 allocate(eta_gcm(lonlen,latlen,altlen,timelen)) 215 216 allocate(tfile(lonlen,latlen,altlen,timelen)) 217 allocate(pfile(lonlen,latlen,altlen,timelen)) 216 218 allocate(tsoilfile(lonlen,latlen,altlen,timelen)) 217 219 allocate(dsoilfile(lonlen,latlen,altlen,timelen)) … … 244 246 allocate(gwparam(lonlen,latlen,5)) 245 247 allocate(ghtsfile(lonlen,latlen)) !! no scan axis 248 allocate(ghfile(lonlen,latlen,altlen,timelen)) 246 249 allocate(vide(lonlen,latlen)) 247 250 allocate(ones(lonlen,latlen)) … … 270 273 271 274 tfile(:,:,:,:)=0 275 pfile(:,:,:,:)=0 272 276 tsoilfile(:,:,:,:)=0 273 277 isoilfile(:,:,:,:)=0 … … 300 304 gwparam(:,:,:)=0 301 305 ghtsfile(:,:)=0 306 ghfile(:,:,:,:)=0 302 307 vide(:,:)=0 303 308 ones(:,:)=0 … … 313 318 SELECT CASE(ident) 314 319 CASE('LMD') 315 ierr = NF_INQ_VARID (nid, " Time",nvarid)320 ierr = NF_INQ_VARID (nid, "time_counter",nvarid) 316 321 CASE('OXF') 317 322 ierr = NF_INQ_VARID (nid, "time",nvarid) … … 331 336 SELECT CASE(ident) 332 337 CASE('LMD') 333 ierr = NF_INQ_VARID (nid, "lat itude",nvarid)338 ierr = NF_INQ_VARID (nid, "lat",nvarid) 334 339 CASE('OXF') 335 340 ierr = NF_INQ_VARID (nid, "lat",nvarid) … … 349 354 SELECT CASE(ident) 350 355 CASE('LMD') 351 ierr = NF_INQ_VARID (nid, "lon gitude",nvarid)356 ierr = NF_INQ_VARID (nid, "lon",nvarid) 352 357 CASE('OXF') 353 358 ierr = NF_INQ_VARID (nid, "lon",nvarid) … … 364 369 print *,lon(1),' ... to ... ',lon(lonlen),' ... step: ',lon(lonlen)-lon(lonlen-1) 365 370 366 SELECT CASE(ident)367 CASE('LMD')368 print *, 'Hybrid coordinates'369 ierr = NF_INQ_VARID (nid, "aps", nvarid)370 IF (ierr .NE. NF_NOERR) THEN371 PRINT *, "Error: Readmeteo <aps> not found"372 stop373 ENDIF374 #ifdef NC_DOUBLE375 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, aps)376 #else377 ierr = NF_GET_VAR_REAL(nid, nvarid, aps)378 #endif379 ierr = NF_INQ_VARID (nid, "bps", nvarid)380 IF (ierr .NE. NF_NOERR) THEN381 PRINT *, "Error: Readmeteo <bps> not found"382 stop383 ENDIF384 #ifdef NC_DOUBLE385 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, bps)386 #else387 ierr = NF_GET_VAR_REAL(nid, nvarid, bps)388 #endif389 print *,aps(1),' ... to ... ',aps(altlen)390 print *,bps(1),' ... to ... ',bps(altlen)391 CASE('OXF')392 ierr = NF_INQ_VARID (nid, "sigma", nvarid)393 394 PRINT *, "Error: Readmeteo <sigma> not found"395 stop396 ENDIF397 #ifdef NC_DOUBLE398 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, bps)399 #else400 ierr = NF_GET_VAR_REAL(nid, nvarid, bps)401 #endif402 print *,bps(1),' ... to ... ',bps(altlen)403 END SELECT371 !SELECT CASE(ident) 372 !CASE('LMD') 373 !print *, 'Hybrid coordinates' 374 ! ierr = NF_INQ_VARID (nid, "aps", nvarid) 375 ! IF (ierr .NE. NF_NOERR) THEN 376 ! PRINT *, "Error: Readmeteo <aps> not found" 377 ! stop 378 ! ENDIF 379 !#ifdef NC_DOUBLE 380 ! ierr = NF_GET_VAR_DOUBLE(nid, nvarid, aps) 381 !#else 382 ! ierr = NF_GET_VAR_REAL(nid, nvarid, aps) 383 !#endif 384 ! ierr = NF_INQ_VARID (nid, "bps", nvarid) 385 ! IF (ierr .NE. NF_NOERR) THEN 386 ! PRINT *, "Error: Readmeteo <bps> not found" 387 ! stop 388 ! ENDIF 389 !#ifdef NC_DOUBLE 390 ! ierr = NF_GET_VAR_DOUBLE(nid, nvarid, bps) 391 !#else 392 ! ierr = NF_GET_VAR_REAL(nid, nvarid, bps) 393 !#endif 394 ! print *,aps(1),' ... to ... ',aps(altlen) 395 ! print *,bps(1),' ... to ... ',bps(altlen) 396 !CASE('OXF') 397 ! ierr = NF_INQ_VARID (nid, "sigma", nvarid) 398 ! IF (ierr .NE. NF_NOERR) THEN 399 ! PRINT *, "Error: Readmeteo <sigma> not found" 400 ! stop 401 ! ENDIF 402 !#ifdef NC_DOUBLE 403 ! ierr = NF_GET_VAR_DOUBLE(nid, nvarid, bps) 404 !#else 405 ! ierr = NF_GET_VAR_REAL(nid, nvarid, bps) 406 !#endif 407 ! print *,bps(1),' ... to ... ',bps(altlen) 408 !END SELECT 404 409 405 410 … … 443 448 444 449 print *,'Surface Pressure' 445 ierr = NF_INQ_VARID (nid,"ps ",nvarid)450 ierr = NF_INQ_VARID (nid,"psol",nvarid) 446 451 IF (ierr .NE. NF_NOERR) THEN 447 452 PRINT *, "Error: Readmeteo <ps> not found" … … 455 460 456 461 print *,'Ground Temperature' 457 ierr = NF_INQ_VARID (nid,"ts urf",nvarid)462 ierr = NF_INQ_VARID (nid,"tsol",nvarid) 458 463 IF (ierr .NE. NF_NOERR) THEN 459 464 PRINT *, "Error: Readmeteo <tsurf> not found" … … 518 523 ! vnsfile=vfileorig(:,:,1,:) 519 524 520 SELECT CASE(ident) 521 CASE('LMD') 522 print *,'Geopotential height at the ground'523 ierr = NF_INQ_VARID (nid," phisinit",nvarid)525 print *,'>>> Read 3D meteorological fields ! - This may take some time ...' 526 527 print *,'Geopotential height' 528 ierr = NF_INQ_VARID (nid,"geop",nvarid) 524 529 IF (ierr .NE. NF_NOERR) THEN 525 PRINT *, "Error: Readmeteo < phisinit> not found"530 PRINT *, "Error: Readmeteo <geop> not found" 526 531 stop 527 532 ENDIF 528 533 #ifdef NC_DOUBLE 529 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, gh tsfile)534 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, ghfile) 530 535 #else 531 ierr = NF_GET_VAR_REAL(nid, nvarid, gh tsfile)536 ierr = NF_GET_VAR_REAL(nid, nvarid, ghfile) 532 537 #endif 533 !**** from geopotential to geopotential height 534 ghtsfile=ghtsfile/grav 535 !**** 536 CASE('OXF') 537 ! 538 ! geopotential height ~ altimetry 539 ! 540 print *,'Geopotential height at the ground from file mountain_new.nc' 541 ierr=NF_OPEN("mountain_new.nc",NF_NOWRITE,nid3) 542 if (ierr.ne.NF_NOERR) then 543 write(*,*) "Error: Could not open that file either" 544 stop "Might as well stop here" 545 endif 546 ierr=NF_INQ_VARID(nid3,"orography",nvarid) 547 ierr=NF_GET_VAR_REAL(nid3,nvarid,ghtsfile) 548 if (ierr.ne.NF_NOERR) then 549 stop "Error: Failed reading phisinit" 550 endif 551 ierr=NF_CLOSE(nid3) 552 END SELECT 553 554 555 print *,'>>> Read 3D meteorological fields ! - This may take some time ...' 538 ghtsfile=ghfile(:,:,1,1)/grav !surface geop 539 540 print *,'Pressure' 541 ierr = NF_INQ_VARID (nid,"pres",nvarid) 542 IF (ierr .NE. NF_NOERR) THEN 543 ierr = NF_INQ_VARID (nid,"p",nvarid) 544 IF (ierr .NE. NF_NOERR) THEN 545 PRINT *, "Error: Readmeteo <p> not found" 546 stop 547 ENDIF 548 ENDIF 549 #ifdef NC_DOUBLE 550 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, pfile) 551 #else 552 ierr = NF_GET_VAR_REAL(nid, nvarid, pfile) 553 #endif 556 554 557 555 print *,'Temperature' … … 572 570 573 571 print *,'Zonal wind' 574 ierr = NF_INQ_VARID (nid," u",nvarid)572 ierr = NF_INQ_VARID (nid,"vitu",nvarid) 575 573 IF (ierr .NE. NF_NOERR) THEN 576 574 PRINT *, "Error: Readmeteo <u> not found" … … 585 583 586 584 print *,'Meridional wind' 587 ierr = NF_INQ_VARID (nid,"v ",nvarid)585 ierr = NF_INQ_VARID (nid,"vitv",nvarid) 588 586 IF (ierr .NE. NF_NOERR) THEN 589 587 PRINT *, "Error: Readmeteo <v> not found" … … 595 593 ierr = NF_GET_VAR_REAL(nid, nvarid, vfile) 596 594 #endif 597 vnsfile=ufile(:,:,1,:) 598 595 vnsfile=vfile(:,:,1,:) 599 596 600 597 !!------------------------ … … 762 759 763 760 !END SELECT 761 762 print*,'VENUS : rotation backward -> inversion of lat and long' 763 DO i=1,latlen 764 DO j=1,lonlen 765 psfile(j,i,:)=psfile(lonlen+1-j,latlen+1-i,:) 766 tsfile(i,j,:)=tsfile(lonlen+1-j,latlen+1-i,:) 767 ghtsfile(j,i)=ghtsfile(lonlen+1-j,latlen+1-i) 768 ghfile(j,i,:,:)=ghfile(lonlen+1-j,latlen+1-i,:,:) 769 pfile(j,i,:,:)=pfile(lonlen+1-j,latlen+1-i,:,:) 770 tfile(j,i,:,:)=tfile(lonlen+1-j,latlen+1-i,:,:) 771 ufile(j,i,:,:)=ufile(lonlen+1-j,latlen+1-i,:,:) 772 vfile(j,i,:,:)=vfile(lonlen+1-j,latlen+1-i,:,:) 773 tsoilfile(j,i,:,:)=tsoilfile(lonlen+1-j,latlen+1-i,:,:) 774 ENDDO 775 ENDDO 764 776 765 777 ierr=NF_CLOSE(nid) … … 876 888 write(1) SLAB 877 889 !print *,'The field '//DESC//' was written to '//output 878 879 890 !------------------------! 880 891 ! >>> Write a variable ! … … 894 905 write(1) SLAB 895 906 !print *,'The field '//DESC//' was written to '//output 896 897 907 !------------------------! 898 908 ! >>> Write a variable ! … … 903 913 DESC='Meridional wind' 904 914 XLVL=200100. 905 SLAB= vnsfile(:,:,time_out(l))915 SLAB=-1*vnsfile(:,:,time_out(l)) !VENUS ratoting backwards : v=-1*v 906 916 ! And now put everything in the destination file 907 917 ! ... Header … … 912 922 write(1) SLAB 913 923 !print *,'The field '//DESC//' was written to '//output 914 915 924 !------------------------! 916 925 ! >>> Write a variable ! … … 1366 1375 DO k = 1,altlen 1367 1376 XLVL=levels(k) 1368 1377 SLAB=tfile(:,:,k,time_out(l)) 1369 1378 ! And now put everything in the destination file 1370 1379 ! ... Header … … 1386 1395 DO k = 1,altlen 1387 1396 XLVL=levels(k) 1388 1397 SLAB=ufile(:,:,k,time_out(l)) 1389 1398 ! And now put everything in the destination file 1390 1399 ! ... Header … … 1394 1403 ! ... Data 1395 1404 write(1) SLAB 1405 !write(1) ufile(:,:,k,time_out(l)) 1396 1406 END DO 1397 1407 !print *,'The field '//DESC//' was written to '//output … … 1406 1416 DO k = 1,altlen 1407 1417 XLVL=levels(k) 1408 SLAB=vfile(:,:,k,time_out(l)) 1418 SLAB=-1*vfile(:,:,k,time_out(l)) ! VENUS rotation backward : v=-1*v 1409 1419 ! And now put everything in the destination file 1410 1420 ! ... Header … … 1430 1440 SELECT CASE(ident) 1431 1441 CASE('LMD') 1432 SLAB= aps(k)+bps(k)*psfile(:,:,time_out(l))1442 SLAB=pfile(:,:,k,time_out(l)) 1433 1443 CASE('OXF') 1434 SLAB=bps(k)*psfile(:,:,time_out(l))1444 !SLAB=bps(k)*psfile(:,:,time_out(l)) 1435 1445 END SELECT 1436 1446 ! And now put everything in the destination file … … 1453 1463 DO k = 1,altlen 1454 1464 XLVL=levels(k) 1455 !!*******1456 !! PROVISOIRE: normalement, il faudrait la hauteur geopotentielle1457 !!*******1458 !!however, not used by initialize_real1459 1465 SELECT CASE(ident) 1460 1466 CASE('LMD') 1461 SLAB= 10000.*alog(610./(aps(k)+bps(k)*psfile(:,:,time_out(l))))1467 SLAB=(ghfile(:,:,k,time_out(l)))/grav 1462 1468 CASE('OXF') 1463 SLAB=10000.*alog(610./(bps(k)*psfile(:,:,time_out(l))))1469 !SLAB=10000.*alog(610./(bps(k)*psfile(:,:,time_out(l)))) 1464 1470 END SELECT 1465 1471 ! And now put everything in the destination file … … 1563 1569 XLVL=levels(k) 1564 1570 SLAB=isoilfile(:,:,k,time_out(l)) 1571 1565 1572 ! And now put everything in the destination file 1566 1573 ! ... Header … … 1710 1717 deallocate(eta_gcm) 1711 1718 deallocate(tfile) 1719 deallocate(pfile) 1712 1720 deallocate(tsoilfile) 1713 1721 deallocate(isoilfile) … … 1737 1745 deallocate(co2icefile) 1738 1746 deallocate(ghtsfile) !! no scan axis 1747 deallocate(ghfile) 1739 1748 deallocate(vide) 1740 1749 deallocate(ones)
Note: See TracChangeset
for help on using the changeset viewer.
