[MITgcm-support] Using wind speed as forcing

[Wenrui Jiang]

Hi,

Thanks for reading this message, and happy Thanksgiving!

I am trying to use the wind speed file as forcing, instead of the wind stress file. but I got this error:

(PID.TID 0000.0001) *** ERROR *** EXF_CHECK: read-in wind-stress but not u,v_wind components
(PID.TID 0000.0001) *** ERROR *** EXF_CHECK_INTERP: Invalid value for: uwind_lon_inc =  1.23456700E+05
(PID.TID 0000.0001) *** ERROR *** EXF_CHECK_INTERP: Fix it or Turn off uwind-interp (interpMethod=0)
(PID.TID 0000.0001) *** ERROR *** EXF_CHECK_INTERP: Invalid value for: vwind_lon_inc =  1.23456700E+05
(PID.TID 0000.0001) *** ERROR *** EXF_CHECK_INTERP: Fix it or Turn off vwind-interp (interpMethod=0)
(PID.TID 0000.0001) ** WARNING ** EXF_CHECK: wind-stress position irrelevant
(PID.TID 0000.0001) *** ERROR *** EXF_CHECK: detected  3 fatal error(s)
(PID.TID 0000.0001) *** ERROR *** S/R ALL_PROC_DIE: ending the run

I guess it might be because of the ALLOW_ATM_WIND flag. So I #undef it in EXF_OPTIONS.h. After recompiling, it still gives me the same error. In the EXF_OPTIONS.h file, there is this description:
C    The following CPP options:
C       ALLOW_ATM_WIND / useAtmWind (useWind)
C       ALLOW_ATM_TEMP               (TEMP)
C       ALLOW_DOWNWARD_RADIATION     (DOWN)
C       ALLOW_BULKFORMULAE           (BULK)
C       EXF_READ_EVAP                (EVAP)
C       ALLOW_READ_TURBFLUXES        (TURB)
C
C    permit all ocean-model forcing configurations listed in the 2 tables below.
C    The first configuration (A1,B1) is the flux-forced, ocean model.
C    Configurations A2,B3 and A2,B4 use pkg/exf open-water bulk formulae
C    to compute, from atmospheric variables, the missing surface fluxes.
C    The forcing fields in the rightmost column are defined in EXF_FIELDS.h
C    (ocean-model surface forcing field are defined in model/inc/FFIELDS.h)
C
C    (A) Surface momentum flux: [model: fu,fv ; exf: ustress,vstress]
C
C    # |useWind|        actions
C   ---|-------|-------------------------------------------------------------
C   (1)| False | Read-in ustress,vstress (if needed in B, compute wind-speed)
C      |       |
C   (2)| True  | Read-in uwind,vwind ; compute wind stress ustress,vstress.
C   ---|-------|-------------------------------------------------------------

But I couldn't find useWind or useAtmWind in runtime options or CPP flags. Could you please tell me what I should do?


Wenrui Jiang

Ph.D. Student,
Department of Earth & Planetary Sciences,
Johns Hopkins University.

Email: wjiang33 at jh.edu

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20211125/8933a881/attachment.html>