error while adding NARR wind forcing to ROMS

[youngstor]

Hi,

I have been running ROMS with tidal forcing only (no wind forcing and surface flux is defined as '#define ANA_SMFLUX'). Now I intend to add wind surface forcing to the model on top of the tidal forcing. I have created separate wind forcing netcdf files (Uwind, Vwind) from NARR data. How should I modify the .in file to include the wind forcing?

My old .in file has something like this (it's a case with grid nesting):

Code: Select all

     NFFILES == 1

     FRCNAME == Projects/test_ca_refined/otps_wc.nc   \
                           Projects/test_ca_refined/otps_wc.nc   \

I simply replaced above with:

Code: Select all

   NFFILES == 2

    FRCNAME == Projects/test_ca_refined_real_wind/otps_wc.nc   \
                Projects/test_ca_refined_real_wind/roms_narr_Jan2009.nc   \
                Projects/test_ca_refined_real_wind/otps_wc.nc   \
                Projects/test_ca_refined_real_wind/roms_narr_Jan2009.nc   

I also replaced ''

Code: Select all

#define ANA_SMFLUX

' with '

Code: Select all

#define BULK_FLUXES

', and re-compiled the program.

However, when I ran the model, I got the error saying '

Code: Select all

GET_2DFLD   - unable to find requested variable: swrad
               in files:
               Projects/test_ca_refined_real_wind/otps_wc.nc
               Projects/test_ca_refined_real_wind/roms_narr_Jan2009.nc

'.

I looked it up and found 'swrad' stands for 'short wave radiation'. How to get around this error without having to create a separate netcdf forcing file for 'swrad'?

Btw, is it ok to set all the other surface forcing (e.g. Pair, Tair, Qair, rain, swrad, lwrad) to zero if I only consider wind forcing for the surface?

Thanks a lot,
Xiufeng

[kate]

You can set ANA_SRFLUX as well as the other ANA_ options. Setting ANA_PAIR to 1000.0 is better than setting it to zero. Zero air pressure does funny things. You might instead opt for creating a new bulk_flux.F that computes stresses as a function of the winds only. I've got such a thing with BULK_FLUX2D in the attached file.