...
There is a small error in OpenIFS version 38r1 in the code that sets the SST for the 3KW1 case of Neale & Hoskins. This code incorrectly sets the latitudinal width to 30o instead of 60o. This error is corrected in later versions of OpenIFS.Expand the box below for details of the necessary code changes.
| Expand |
|---|
| title | Expand to see required code changes... |
|---|
|
The required code changes are small. In the file: src/ifs/phys_ec/gp_sstaqua.F90, find the lines 233:238: | Code Block |
|---|
! 3KW1
RCHI=3.0_JPRB
IF( ABS(PGEMU(IOFF+JROF-1)) <= SIN(RPHID) ) THEN
PSD_VF(JROF,YSD_VF%YSST%MP,IBL)=PSD_VF(JROF,YSD_VF%YSST%MP,IBL)+&
&RCHI*COS(PGELAM(IOFF+JROF-1)-RLAMBDA0)*&
&COS(0.5_JPRB*ZPI*(ASIN(PGEMU(IOFF+JROF-1))/RPHID))**2 |
The use of RPHID needs to be corrected to 60o: | Code Block |
|---|
! 3KW1
RCHI=3.0_JPRB
IF( ABS(PGEMU(IOFF+JROF-1)) <= SIN(2._JPRB*RPHID) ) THEN ! correct use of RPHID on this line
PSD_VF(JROF,YSD_VF%YSST%MP,IBL)=PSD_VF(JROF,YSD_VF%YSST%MP,IBL)+&
&RCHI*COS(PGELAM(IOFF+JROF-1)-RLAMBDA0)*&
&COS(0.5_JPRB*ZPI*(ASIN(PGEMU(IOFF+JROF-1))/(2._JPRB*RPHID)))**2 ! correct use of RPHID on this line |
|
...