I was trying to make a time-dependent Ocean file which included time-dependent variables like OPEN, SURF, CHLO, DMS and their related TFLAG, in order to use daily-resolution DMS grid data to replace the monthly averaged version. However, when I tried to use this time-dependent oceanfile in my CMAQ simulation, it appeared that I encountered the error message “*** ERROR ABORT in subroutine retrieve_ocean_d on PE 026 Could not read OPEN from OCEAN_1" in the log. In this OCEAN FILE, the SDATE was ‘2017198‘; 25 days DMS data were incorporated; for OPEN and SURF, I just copied the first day result 25 times. Detailed information about my Ocean File are listed as follows.
I also noticed a similar issue solved by @barronh . Could you give me some suggestions to run simulation with a time-dependent Ocean file?
My interpretation of the CMAQ code reading OCEAN files is that it is designed for time-independent. The centralized_io_mod reads data using the XTRACT3 command with JDATE=0 and JTIME=0. Obviously, the JDATE and JTIME are not equal to 0, but the XTRACT3 subroutine ignores JDATE and JTIME for time-independent files and returns.
Given the data sources being used to construct OCEAN files, day specific should be sufficient. Since CMAQ is typically run one day at a time, this means you can use time-independent files that are created for each day.
If this does not solve your problem, please be specific about what you need.
Thanks for your suggestion! I managed to change the JDATE and JTIME in centralized_io_mod.F to make Oceanfile time-dependent, it really worked well with CCTM simulation.
