Test NEMS workflow with v55 adcirc for ATL+PR 120 m mesh for the stroms about coastalapp

[saeed-moghimi-noaa]

Assigned to @ps2epsi

[saeed-moghimi-noaa]

Hi @mtcontre and @WPringle
Where we can get the ATL+PR 120 m mesh NAVD88 +10 m mesh?
Thanks,
-Saeed

[mtcontre]

Hi Saeed,
Can you please check this mesh:
https://drive.google.com/drive/folders/1WZm30VWkbsINHslfV7qFIaHhUqHUF7Fm?usp=sharing

It is EC+PR, but I'm not sure if it is clipped to 6 or 10 m. Otherwise, the last release that we used for GESTOFS is in the following link:
https://drive.google.com/drive/folders/1yQLhdzDDqWMftm2H3GQlKRQI8-86BfbG?usp=sharing
However, that mesh doesn't have the insert in PR.
Best!

[mtcontre]

Hi Saeed,
Can you please check this mesh:
https://drive.google.com/drive/folders/1WZm30VWkbsINHslfV7qFIaHhUqHUF7Fm?usp=sharing

It is EC+PR, but I'm not sure if it is clipped to 6 or 10 m. Otherwise, the last release that we used for GESTOFS is in the following link:
https://drive.google.com/drive/folders/1yQLhdzDDqWMftm2H3GQlKRQI8-86BfbG?usp=sharing
However, that mesh doesn't have the insert in PR.
Best!

So the first link has the EC mesh up to 10 m. It needs the insert in PR.
Best!

[saeed-moghimi-noaa]

Thanks All.

@mtcontre
Please share the latest 120 m with 10m above NAVD88 with William for inclusion of PR. Please make sure it is with +10m
land elevation.

@WPringle
Please send the final PAC and ATL+PR meshes to @JaimeCalzadaNOAA for testing.

[saeed-moghimi-noaa]

Hi Yuji @ps2epsi
You should be able to run the previous HWRF in ATL through cap.  The Cap will read the ATM data and interpolation should happen on the fly. You just need to create a new:

https://github.com/NOAA-EMC/NSEM-workflow/blob/master/parm/nems.configure.atm_ocn_wav_1loop.IN
file and edit the time loop part to:

runSeq::
  @$_coupling_interval_sec_
    ATM -> OCN   
    WAV -> OCN   :remapMethod=redist
    ATM -> WAV   
    OCN -> WAV   :remapMethod=redist
    ATM
    OCN
    WAV
  @
::

[ps2epsi]

Saeed - I have a question about remapMethod. If remapMethod=redist, the simulation was failed by the following errors (in PET???.ESMF_LogFile)

20210610 111632.582 ERROR PET000 ESMCI_Array.C:6046

ESMCI::Array::tRedistStore() Invalid argument - srcArray and dstArray must provide identical number of exclusive elements

20210610 111632.591 ERROR PET000 ESMCI_Array.C:5827

ESMCI::Array::redistStore() Invalid argument - Internal subroutine call returned Error

20210610 111632.591 ERROR PET000 ESMCI_ArrayBundle.C:991

ESMCI::ArrayBundle::redistStore( Invalid argument - Internal subroutine call returned Error

20210610 111632.591 ERROR PET000 ESMCI_ArrayBundle_F.C:494

c_esmc_arraybundlerediststorenf( Invalid argument - Internal subroutine call returned Error If remapMethod=regrid, the simulation did not fail, but the warning message was shown:

20210610 121719.611 WARNING PET000 Specified option

'remapmethod=regrid' is not a vailid choice. Defaulting to BILINEAR for: 'inst_merid_wind_height10m'

20210610 123224.630 WARNING PET000 Specified option

'remapmethod=regrid' is not a vailid choice. Defaulting to BILINEAR for: 'inst_zonal_wind_height10m' Is this a good remapMethod? Do you know other remapMethod? Please let me know. Thanks, Yuji

…

On Tue, Jun 8, 2021 at 1:27 PM Saeed Moghimi ***@***.***> wrote: Hi Yuji You should be able to run the previous HWRF in ATL through cap. The Cap will read the ATM data and interpolation should happen on the fly. You just need to create a new: NSEM-workflow/parm/nems.configure.atm_ocn_wav_1loop.IN file and edit the time loop part to: runSeq:: @$_coupling_interval_sec_ ATM -> OCN WAV -> OCN :remapMethod=redist ATM -> WAV OCN -> WAV :remapMethod=redist ATM OCN WAV @ :: — You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub <#39 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AE6P45WLIJIYL5DCXAE4N7LTRZHIHANCNFSM456Y2NMQ> .

-- Yuji Funakoshi, Ph.D. Coast Survey Development Laboratory National Oceanic and Atmospheric Administration 1315 East West Highway - N/CS13 Silver Spring, MD 20910-3282 ph. [phone] ***@***.***

[saeed-moghimi-noaa]

Hi Yuji @ps2epsi

Please see the NUOPC manual here:
https://earthsystemmodeling.org/docs/nightly/develop/NUOPC_refdoc/node3.html#SECTION00033400000000000000

you can use redist only if the mesh on both side are the same. In this case HWRF is interpolated on the old HSOFS and OCN is on the 120m new ATL mesh. Therefore you need to change the remapmethod.

I suggest to go with bilinear or one of the nearest (I do not know the difference) please look into the manual.

As I mentioned in the issue above you need to create a new:

https://github.com/NOAA-EMC/NSEM-workflow/blob/master/parm/nems.configure.atm_ocn_wav_1loop.IN

and make these changes:

runSeq::
  @$_coupling_interval_sec_
    ATM -> OCN   :remapMethod=bilinear
    WAV -> OCN   :remapMethod=redist
    ATM -> WAV   :remapMethod=bilinear
    OCN -> WAV   :remapMethod=redist
    ATM
    OCN
    WAV
  @
::

Thanks,
-Saeed

[ps2epsi]

Saeed - Thank you for your advice. So far, the simulation is still running with "regrid", which uses the default. If it works without any problems, I will run the hindcast for other storms. Thanks, Yuji

…

On Thu, Jun 10, 2021 at 10:10 AM Saeed Moghimi ***@***.***> wrote: Hi Yuji @ps2epsi <https://github.com/ps2epsi> Please see the NUOPC manual here: https://earthsystemmodeling.org/docs/nightly/develop/NUOPC_refdoc/node3.html#SECTION00033400000000000000 [image: image] <https://user-images.githubusercontent.com/66274703/121531463-6b8f5900-c9cc-11eb-8f05-f690ae4168fa.png> you can use redist only if the mesh on both side are the same. In this case HWRF is interpolated on the old HSOFS and OCN is on the 120m new ATL mesh. Therefore you need to change the remapmethod. I suggest to go with bilinear or one of the nearest (I do not know the difference) please look into the manual. runSeq:: @$_coupling_interval_sec_ ATM -> OCN :remapMethod=bilinear WAV -> OCN :remapMethod=redist ATM -> WAV :remapMethod=bilinear OCN -> WAV :remapMethod=redist ATM OCN WAV @ :: Thanks, -Saeed — You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub <#39 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AE6P45RT3XJCH27IRKRWOT3TSDBV3ANCNFSM456Y2NMQ> .

-- Yuji Funakoshi, Ph.D. Coast Survey Development Laboratory National Oceanic and Atmospheric Administration 1315 East West Highway - N/CS13 Silver Spring, MD 20910-3282 ph. [phone] ***@***.***

[ps2epsi]

I've tested all of remapMethod, and here is my recommendation based on simulation runtime.

1. redist
2. nearest_stod
3. bilinear
4. patch and nearest_dtos
5. conserve
   "redist" is the fastest. "nearest_stod" is the best method if "redist" isn't available. "bilinear" works fine, but it's not good for high resolution model grid because runtime depends on grid number. "patch and nearest_dtos" is the same as "bilinear", but they are slower than "bilinear". "conserve" didn't work for this time.

[saeed-moghimi-noaa]

Harvey and Maria run finished for

ATM -> ADC
ATM + WW3 -> ADC