Introduction
User subroutine for the atmospheric model.
Atmospheric module
imode
corresponds to the number of calls of the usatdv function. Depending on the value of imode
, different operations are performed in the following example.
imode = 0
if (imode.eq.0) then
write(
nfecra,*)
'defining the dimensions of the 1D vertical arrays'
else
write(
nfecra,*)
'defining the coordinates and levels of the 1D vertical arrays'
endif
if (imode.eq.0) then
zvmax = 1975.d0
ztop = 11000.d0
zzmax = (int(zvmax)/1000)*1000.d0
do while(zzmax.le.(ztop-1000.d0))
zzmax = zzmax + 1000.d0
enddo
endif
imode = 1
ztop = 11000.d0
zzmax = (int(
zvert(ii))/1000)*1000.d0
do while(zzmax.le.(ztop-1000.d0))
zzmax = zzmax+1000.d0
ii = ii + 1
enddo
endif
enddo
endif
Data Entry for the atmospheric ground model
ifbt1d = 0
if (iappel.eq.1) then
call
getfbr(
'75',nlelt,lstelt)
do ilelt = 1, nlelt
ifbt1d = ifbt1d + 1
enddo
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
enddo
endif
if (iappel.eq.2) then
enddo
enddo
endif
deallocate(lstelt)