run error when applying forcing file

[javadi]

Hi users,
I have faced a FORTRAN runtime error in running time. Firstly, I made the grid for my region. The model run without any problem when applying my real grid. Then, I created my forcing file using d_ecmwf2roms script. I have defined 13 files in upwelling.in file and active smflux in upwelling.h. Here is snapshot of files. The model compile without any error; however, when I run, I have this error.

Please help me with your interesting comments

[kate]

I see that you turned off ANA_SMFLUX. Did you turn on anything else to compensate? I use BULK_FLUXES with:

Code: Select all

# if defined BULK_FLUXES || defined CCSM_FLUXES
#  define LONGWAVE_OUT
#  undef DIURNAL_SRFLUX
#  define SOLAR_SOURCE
#  define EMINUSP
#  undef ALBEDO_CLOUD
#  define ALBEDO_CURVE  /* for water */
#  undef ICE_ALB_EC92   /* for ice */
#  define ALBEDO_CSIM   /* for ice */
#  undef ALBEDO_FILE    /* for both */
#  undef ALBEDO_DIRDIFF
#  define ICE_SHORTWAVE_R
#  undef LONGWAVE
# endif

You probably don't need all that albedo stuff without sea ice.

[javadi]

HI,
when it comes to your question, Here is my file. You mean that the error comes from ....h file?

[kate]

Maybe? I don't see BULK_FLUXES in your .h file, so without ANA_SMFLUX, the model will only read your sms_era file, not many of the rest of them.
There may be other errors in your roms.in file though. The line:

Code: Select all

allocate ( FRC(max_Ffiles,Ngrids) )

happens once in read_phypar.F, so I assume that's what's on line 2176 of your read_phypar.f90. It gets invoked when it reads the NFFILES line of roms.in. Do you have more than one such line? There needs to be one and only one such line.

[javadi]

Thank you so much Kate.
Yes there is one line.
To be honest, I do not understand exactly how to turn on flag. Can you correct my upwelling.h? I try to use smflux. Now i should use my own grid named roms_grd.nc and forcing files. How active them in upwelling.h?
please help me

[kate]

You can see my various .h and .in files here:
https://github.com/kshedstrom/Apps
Of course, they go with my ROMS branch, so you have to know what you are doing to pick and choose features. I've described things as best I can here:
https://github.com/kshedstrom/roms_manu ... manual.pdf

[jcwarner]

check your ocean.in file make sure you dont have nffiles twice
NFFILES ==
...
NFFILES ==

the FRC array is allocated in read_phypar
CASE ('NFFILES')
...
allocate ( FRC(max_Ffiles,Ngrids) )

so it seems to be doing that call twice.
-j

[hpd14thu]

From your .h file, I find you use ANA_SMFLUX, but you undefine ANA_SSFLUX and STFLUX.
Well, define ana_.. means you use anaytical functions in ROMS/fucntions to calculate this variable. However, it seems that you want to calculate them by Monin-Obukhov similarities method in bulk.F but not anaytical functions.
Maybe defining BULK_FLUXES and undefine ana_... like Kate shown above is better for your case.