Opened 5 years ago

Closed 4 years ago

#220 closed defect (wontfix)

Problem with negative precipitation

Reported by: maignan Owned by: somebody
Priority: major Milestone:
Component: Anthropogenic processes Version:
Keywords: Cc:

Description

It happens that some forcing files have some slightly negative precipitation. In this case the program stops in hydrol_soil_infilt with the message:

Error in the calculation of infilt tot -4.176114610371373E-008
 k, ji, jst, mc   1.17784270216726        2.16374183793683             2394
           1  0.277241865204730

FATAL ERROR FROM ROUTINE hydrol_soil_infilt
 --> We will STOP after hydrol_soil_infilt.
 -->
 -->

Fatal error from IOIPSL. STOP in ipslerr with code

We have to decide if we correct the forcing files or if we modify the code.
Such a modification could be made in intersurf_initialize_2d/intersurf_main_2d:

!@ PROPOSED MODIFICATION
    WHERE(zprecip_rain(:) .LT. 0.)
	zprecip_rain(:)=0.		
    ENDWHERE				
!@ END PROPOSED MODIFICATION
 !
    IF (check_INPUTS) THEN
...

Change History (2)

comment:1 Changed 5 years ago by jpolcher

For me this should clearly be solved in the forcing.
How can a forcing file provide negative rainfall ... this makes no physical sense.
We can put some test into the driver, but that is the most I would concede.

Jan

comment:2 Changed 4 years ago by maignan

  • Resolution set to wontfix
  • Status changed from new to closed
Note: See TracTickets for help on using tickets.