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 soil model
To activate the model, the user has to set the and to specify the zone id on which the soil model is applied in cs_user_parameters.c, routine cs_user_model:
Then the user may change default coefficient values for soil categories in cs_user_parameters.f90, routine usipsu
Initialization of atmospheric
The user has to specify the percentage of each categories for all faces of the soil zone in cs_user_initialization.c:
f->
val[soil_id + f->
dim * elt_id] = 0.;
f->
val[4 + f->
dim * elt_id] = 100.;
}
}