Opened 4 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.