Opened 5 years ago
#64 new defect
day fraction "fract" can only be 0 or 1 ?
| Reported by: | emillour | Owned by: | |
|---|---|---|---|
| Priority: | major | Milestone: | |
| Component: | MARS GCM | Version: | |
| Keywords: | Cc: |
Description
In the Mars physics (Generic also), variable "fract(:)" seems (computations in routine solang) to refer to the fraction of the cell in which it is daytime but in practice it is set to either 0 or 1:
DO jl=1,kgrid
IF (pmu0(jl).gt.0.) THEN
pfract(jl)=1.
c pmu0(jl)=sqrt(1224.*pmu0(jl)*pmu0(jl)+1.)/35.
ELSE
c pmu0(jl)=0.
pfract(jl)=0.
ENDIF
ENDDO
The variable is stored as "fract" in "comsaison_h".
Should it not be allowed to vary between 0 and 1 for meshes in partial day and night? (it seems to be the case in the Earth and Venus models)
Note: See
TracTickets for help on using
tickets.
