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)

Change History (0)

Note: See TracTickets for help on using tickets.