https://github.com/alphaparrot/ExoPlaSim
Tip revision: 9c1acdb9139eea51193af50fb94016a3f9e2543d authored by Adiv Paradise on 17 October 2017, 20:17:43 UTC
Fixed bug with fixedlon setting. There seem to be output bugs, and
Fixed bug with fixedlon setting. There seem to be output bugs, and
Tip revision: 9c1acdb
dump_ice.f90
! ===================
! SUBROUTINE ICE_DUMP
! ===================
subroutine ice_dump
use icemod
character (len=256) :: yfn
call step2cal(nstep,ntspd,iyea,imon,iday,ihou,imin)
write(yfn,'("ice_dump_",I4.4,4I2.2)') iyea,imon,iday,ihou,imin
open(72,file=yfn)
100 format(a16,i16)
write (72,100) 'nice',nice
write (72,100) 'newsurf',newsurf
write (72,100) 'nsnow',nsnow
write (72,100) 'ntskin',ntskin
write (72,100) 'ntspd',ntspd
write (72,100) 'noutput',noutput
write (72,100) 'nout',nout
write (72,100) 'nfluko',nfluko
write (72,100) 'ncpl_ice_ocean',ncpl_ice_ocean
write (72,100) 'nperpetual_ice',nperpetual_ice
write (72,100) 'nprint',nprint
write (72,100) 'nprhor',nprhor
write (72,100) 'nstep',nstep
write (72,100) 'nrestart',nrestart
write (72,100) 'naccuout',naccuout
200 format(a16,e20.10)
write (72,200) 'taunc',taunc
write (72,200) 'xmind',xmind
write (72,200) 'xdt',xdt
write (72,200) 'cicemin',cicemin
write (72,200) 'solar_day',solar_day
! call ice_dump_hor('xshfl',xshfl)
! call ice_dump_hor('xshdt',xshdt)
! call ice_dump_hor('xlhfl',xlhfl)
! call ice_dump_hor('xlhdt',xlhdt)
! call ice_dump_hor('xswfl',xswfl)
! call ice_dump_hor('xlwfl',xlwfl)
call ice_dump_hor('xts',xts)
call ice_dump_hor('xsst',xsst)
call ice_dump_hor('xmld',xmld)
call ice_dump_hor('xls',xls)
call ice_dump_hor('xiced',xiced)
call ice_dump_hor('xicec',xicec)
call ice_dump_hor('xsnow',xsnow)
! call ice_dump_hor('xsmelt',xsmelt)
! call ice_dump_hor('xsmflx',xsmflx)
! call ice_dump_hor('ximelt',ximelt)
! call ice_dump_hor('xsndch',xsndch)
call ice_dump_hor('xqmelt',xqmelt)
! call ice_dump_hor('xheat',xheat)
! call ice_dump_hor('xcflux',xcflux)
! call ice_dump_hor('xprs',xprs)
! call ice_dump_hor('xfresh',xfresh)
call ice_dump_hor('xroff',xroff)
! call ice_dump_hor('xtaux',xtaux)
! call ice_dump_hor('xtauy',xtauy)
! call ice_dump_hor('xust3',xust3)
call ice_dump_hor('xoflux',xoflux)
! call ice_dump_hor('xtsflux',xtsflux)
! call ice_dump_hor('xhnew',xhnew)
call ice_dump_hor('cheat',cheat)
call ice_dump_hor('cfresh',cfresh)
call ice_dump_hor('ctaux',ctaux)
call ice_dump_hor('ctauy',ctauy)
call ice_dump_hor('cust3',cust3)
call ice_dump_hor('xflxicea',xflxicea)
call ice_dump_hor('xheata',xheata)
call ice_dump_hor('xofluxa',xofluxa)
call ice_dump_hor('xqmelta',xqmelta)
call ice_dump_hor('xcfluxa',xcfluxa)
call ice_dump_hor('xsmelta',xsmelta)
call ice_dump_hor('ximelta',ximelta)
call ice_dump_hor('xtsfluxa',xtsfluxa)
call ice_dump_hor('xhnewa',xhnewa)
call ice_dump_hor14('xclsst',xclsst)
call ice_dump_hor14('xclicec',xclicec)
call ice_dump_hor14('xcliced',xcliced)
call ice_dump_hor14('xflxice',xflxice)
close(72)
end
subroutine ice_dump_hor(ya,pa)
use icemod
character (len=*) :: ya
real :: pa(NHOR)
write(72,'(/A)') trim(ya)
do j = 1 , NHOR , 4
write(72,'(i6,2x,4e18.10)') j,pa(j:j+3)
enddo
return
end
subroutine ice_dump_hor14(ya,pa)
use icemod
character (len=*) :: ya
real :: pa(NHOR,14)
do mon = 1 , 14
write(72,'(/I2,2X,A)') mon-1,trim(ya)
do j = 1 , NHOR , 4
write(72,'(i6,2x,4e18.10)') j,pa(j:j+3,mon)
enddo
enddo
return
end