run error when applying forcing file

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
javadi
Posts: 115
Joined: Mon Jan 27, 2020 5:55 pm
Location: UoT

run error when applying forcing file

#1 Unread post by 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
Attachments
IMG_20201031_203700.jpg
IMG_20201031_202951.jpg
IMG_20201031_202923.jpg

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: run error when applying forcing file

#2 Unread post by 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
Posts: 115
Joined: Mon Jan 27, 2020 5:55 pm
Location: UoT

Re: run error when applying forcing file

#3 Unread post by javadi »

HI,
when it comes to your question, Here is my file. You mean that the error comes from ....h file?
Attachments
upwelling.h
(1.65 KiB) Downloaded 160 times

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: run error when applying forcing file

#4 Unread post by 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
Posts: 115
Joined: Mon Jan 27, 2020 5:55 pm
Location: UoT

Re: run error when applying forcing file

#5 Unread post by 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

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: run error when applying forcing file

#6 Unread post by 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
Posts: 1173
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: run error when applying forcing file

#7 Unread post by 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

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: run error when applying forcing file

#8 Unread post by 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.
Attachments
h.png
h.png (76.8 KiB) Viewed 4233 times

Post Reply