noaa-ocs-modeling/PaHM

Setup and run PaHM for `BEST` and `OFCL` tracks of florence

Closed this issue · 12 comments

FariborzDaneshvar-NOAA commented

PaHM was downloaded and compiled on NHC_Colab_2 cluster on PW.
directory: /home/Fariborz.Daneshvar/PaHM

The following commands were used for compilation:

  • module load intel/2022.1.2 netcdf cmake
  • ./build.sh --compiler intel --platform pworks

Then a copy of ./templates/pahm_control-ike.in was created in the PaHM directory and was named pahm_control_florence.in.

For each scenario, the following lines were updated in the pahm_control_florence.in: bestTrackFileName, meshFileName, outFileName, begDateTime, endDateTime, and modelType = 10

Finally, the following command was used to run each scenario:

  • ./pahm pahm_control_florence.in
FariborzDaneshvar-NOAA commented

@pvelissariou1 I am using pahm <control_file_name> command in the pahm directory, as you suggested, but getting pahm: command not found` error (here is a screenshot)

Image

pvelissariou1 commented
  • 👍1
FariborzDaneshvar-NOAA commented

Florence 2018, BEST track run

The following lines were updated in pahm_control_florence.in:

  • bestTrackFileName = /lustre/hurricanes/florence_2018_Fariborz_v4_15_ensembles/nhc_track/hurricane-track.dat
  • meshFileName = /lustre/hurricanes/florence_2018_Fariborz_v4_15_ensembles/mesh/mesh_w_bdry.grd
  • outFileName = pahm_windout_florence_BEST
  • begDateTime = 2018-08-30 06:00:00
  • endDateTime = 2018-09-18 12:00:00

Model output has signals of the storm track:

Image

FariborzDaneshvar-NOAA commented

Florence 2018, OFCL track run

The following lines were updated in pahm_control_florence.in:

  • bestTrackFileName = /lustre/hurricanes/florence_2018_Fariborz_v4_15_ensembles/nhc_track/hurricane-track.dat
  • meshFileName = /lustre/hurricanes/florence_2018_Fariborz_v3_39_ensembles/mesh/mesh_w_bdry.grd
  • outFileName = pahm_windout_florence_OFCL
  • begDateTime = 2018-08-30 12:00:00
  • endDateTime = 2018-09-12 18:00:00

This run failed with the following message:

(perturbation) [Fariborz.Daneshvar@nhccolab2-3 PaHM]$ ./pahm pahm_control_florence.in
The contol file : pahm_control_florence.in was found and will be opened for reading.
Finished processing the input fields from the control file ...
ERROR :: ProcessAsymmetricVortexData: 23 isotachs were nonzero.
ERROR :: ProcessAsymmetricVortexData: 23 isotachs were nonzero.
ERROR :: ProcessAsymmetricVortexData: 23 isotachs were nonzero.
ERROR :: ProcessAsymmetricVortexData: 23 isotachs were nonzero.
ERROR :: GetGAHMFields: The GAHM asymmetric data structure has more than 4 iSotachs in cycle 54.
ERROR :: Terminate: PaHM Terminating.
FariborzDaneshvar-NOAA commented

Comparison of BEST & OFCL track files created by the ondemand-storm-workflow runs:

Here are first few lines of two nhc_track files for comparison

Florence, 2018 - BEST:

AL, 06, 2018083006,   , BEST,   0, 128N,  169W,  20, 1008, LO,   0,    ,    0,    0,    0,    0, 1010,  150,  50,   0,   0,   L,   0,    ,   0,   0,     INVEST,,,,,,,,
AL, 06, 2018083012,   , BEST,   0, 128N,  179W,  25, 1007, LO,   0,    ,    0,    0,    0,    0, 1010,  150,  50,   0,   0,   L,   0,    , 270,   5,        SIX,,,,,,,,
AL, 06, 2018083018,   , BEST,   0, 128N,  190W,  25, 1007, LO,   0,    ,    0,    0,    0,    0, 1010,  150,  50,  35,   0,   L,   0,    , 270,   6,        SIX,,,,,,,,
AL, 06, 2018083100,   , BEST,   0, 131N,  202W,  30, 1006, LO,   0,    ,    0,    0,    0,    0, 1010,  150,  40,  40,   0,   L,   0,    , 284,   6,        SIX,,,,,,,,
AL, 06, 2018083106,   , BEST,   0, 134N,  214W,  30, 1006, LO,   0,    ,    0,    0,    0,    0, 1010,  150,  40,  40,   0,   L,   0,    , 284,   6,        SIX,,,,,,,,
AL, 06, 2018083112,   , BEST,   0, 136N,  226W,  30, 1006, LO,   0,    ,    0,    0,    0,    0, 1010,  150,  40,  40,   0,   L,   0,    , 280,   6,        SIX,,,,,,,,
AL, 06, 2018083118,   , BEST,   0, 138N,  238W,  30, 1006, TD,   0,    ,    0,    0,    0,    0, 1010,  150,  40,  40,   0,   L,   0,    , 280,   6,        SIX,,,,,,,,
AL, 06, 2018090100,   , BEST,   0, 140N,  249W,  30, 1006, TD,   0,    ,    0,    0,    0,    0, 1010,  150,  40,  40,   0,   L,   0,    , 281,   6,        SIX,,,,,,,,
AL, 06, 2018090106,   , BEST,   0, 143N,  261W,  35, 1005, TS,  34, NEQ,   30,   30,    0,   30, 1010,  150,  30,  45,   0,   L,   0,    , 285,   6,        SIX,,,,,,,,
AL, 06, 2018090112,   , BEST,   0, 148N,  272W,  35, 1004, TS,  34, NEQ,   30,   30,    0,   30, 1010,  150,  30,  45,   0,   L,   0,    , 295,   6,   FLORENCE,,,,,,,,

Florence, 2018 - OFCL:

AL, 06, 2018083012, 03, OFCL,   0, 128N,  179W,  25, 1007, LO,  34, NEQ,    0,    0,    0,    0, 1010,    0,  30,  35,   0,    ,   0, LAA,   0,   0,           ,,,,,,,,
AL, 06, 2018083018, 03, OFCL,   0, 128N,  190W,  25, 1007, LO,  34, NEQ,    0,    0,    0,    0, 1010,    0,  30,  35,   0,    ,   0, LAA, 270,   6,           ,,,,,,,,
AL, 06, 2018083100, 03, OFCL,   0, 131N,  202W,  30, 1007, LO,  34, NEQ,    0,    0,    0,    0, 1010,    0,  40,  40,   0,    ,   0, JLB, 284,   6,           ,,,,,,,,
AL, 06, 2018083106, 03, OFCL,   0, 134N,  210W,  30, 1006, LO,  34, NEQ,    0,    0,    0,    0, 1009,    0,  40,  40,   0,    ,   0, RJB, 291,   4,           ,,,,,,,,
AL, 06, 2018083112, 03, OFCL,   0, 135N,  221W,  30, 1006, LO,  34, NEQ,    0,    0,    0,    0, 1009,    0,  40,  40,   0,    ,   0, LAA, 275,   6,           ,,,,,,,,
AL, 06, 2018083118, 03, OFCL,   0, 138N,  240W,  30, 1004, TD,  34, NEQ,    0,    0,    0,    0, 1009,    0,  40,  40,   0,    ,   0, LAA, 279,  10,           ,,,,,,,,
AL, 06, 2018090100, 03, OFCL,   0, 140N,  249W,  30, 1004, TD,  34, NEQ,    0,    0,    0,    0, 1010,    0,  40,  40,   0,    ,   0, JLB, 283,   5,           ,,,,,,,,
AL, 06, 2018090106, 03, OFCL,   0, 143N,  261W,  35, 1003, TS,  34, NEQ,   30,   30,    0,   30, 1009,    0,  30,  45,   0,    ,   0, RJB, 285,   6,           ,,,,,,,,
AL, 06, 2018090112, 03, OFCL,   0, 146N,  273W,  35, 1003, TS,  34, NEQ,   30,   30,    0,   30, 1009,    0,  30,  45,   0,    ,   0, LAA, 285,   6,           ,,,,,,,,
AL, 06, 2018090118, 03, OFCL,   0, 154N,  283W,  40, 1002, TS,  34, NEQ,   40,   40,   20,   40, 1009,    0,  30,  50,   0,    ,   0, LAA, 310,   6,           ,,,,,,,,
pvelissariou1 commented
FariborzDaneshvar-NOAA commented

@SorooshMani-NOAA there are three storm track files in the nhc_track/ directory:
hurricane-track.dat, hurricane-track.dat_adj, and hurricane-track.dat_asymvort

The first one (hurricane-track.dat) was used in two runs mentioned above. Which one do you want me to use for new PaHM runs?

pvelissariou1 commented

@FariborzDaneshvar-NOAA You need to use hurricane-track.dat, the other two are generated from PaHM to be used by ADCIRC basically. hurricane-track.dat_adj is generated when using the Holland model in PaHM and hurricane-track.dat_asymvort is generated when using the GAHM model in PaHM.

  • 👍1
SorooshMani-NOAA commented

@FariborzDaneshvar-NOAA I'm talking about the one in setup/ensemble.dir/runs/track_files/original.22. The nhc_track/hurricane-track.dat is the one before perturbation. The other two are generated after PaHM processes the input, so they are actually PaHM outputs.

So this is where I would find the 3 files I talked about:

  • Actual Best track: <best-track-rundir>/nhc_track/hurricane-track.dat
  • Modified OFCL track: <ofcl-track-rundir>/nhc_track/hurricane-track.dat
  • Fake Best track: <ofcl-track-rundir>/setup/ensemble.dir/runs/track_files/original.22
FariborzDaneshvar-NOAA commented

@pvelissariou1 and @SorooshMani-NOAA, I ran PaHM for these four track files:

  1. Actual BEST track: hurricane-track_BEST.dat (path: /lustre/hurricanes/florence_2018_Fariborz_BEST_10_v2/nhc_track/hurricane-track.dat)
  2. BEST track used for SCHISM run: original_BEST.22 (path: /lustre/hurricanes/florence_2018_Fariborz_BEST_10_v2/setup/ensemble.dir/track_files/original.22)
  3. Modified OFCL track: hurricane-track_OFCL.dat (path: /lustre/hurricanes/florence_2018_Fariborz_OFCL_10_v2/nhc_track/hurricane-track.dat)
  4. Fake BEST track for SCHISM run: original_OFCL.22 (path: /lustre/hurricanes/florence_2018_Fariborz_OFCL_10_v2/setup/ensemble.dir/track_files/original.22)

Only the third run (for hurricane-track_OFCL.dat) failed with the following error message:

The contol file : pahm_control_florence_OFCL.in was found and will be opened for reading.
Finished processing the input fields from the control file ...
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source
libifcoremt.so.5   00002B2E0B6DB739  for__signal_handl     Unknown  Unknown
libpthread-2.17.s  00002B2E0EF75630  Unknown               Unknown  Unknown
libifcoremt.so.5   00002B2E0B74563B  for_len_trim          Unknown  Unknown
libifcoremt.so.5   00002B2E0B7455BF  for_trim              Unknown  Unknown
pahm               0000000000424CC9  Unknown               Unknown  Unknown
pahm               0000000000424847  Unknown               Unknown  Unknown
pahm               0000000000450215  Unknown               Unknown  Unknown
pahm               000000000048A18F  Unknown               Unknown  Unknown
pahm               000000000048A275  Unknown               Unknown  Unknown
pahm               0000000000402492  Unknown               Unknown  Unknown
libc-2.17.so       00002B2E0F1A4555  __libc_start_main     Unknown  Unknown
pahm               0000000000402399  Unknown               Unknown  Unknown
FariborzDaneshvar-NOAA commented

@pvelissariou1, I removed the last 20 lines of the hurricane-track_OFCL.dat with non-zero forecast hours (also copied below) and ran it again. This time it worked fine, Thanks!

AL, 06, 2018091218, 03, OFCL,   3, 309N,  725W, 105,  949, HU,  34, NEQ,  170,  140,  100,  140, 1013,    0,  15, 130,  25,    ,   0, SRS, 314,   7,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,   3, 309N,  725W, 105,  949, HU,  50, NEQ,  100,   80,   60,   80, 1013,    0,  15, 130,  25,    ,   0, SRS, 314,   7,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,   3, 309N,  725W, 105,  949, HU,  64, NEQ,   60,   60,   40,   60, 1013,    0,  15, 130,  25,    ,   0, SRS, 314,   7,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  12, 321N,  741W, 110,  948, HU,  34, NEQ,  170,  140,  100,  140, 1013,    0,  15, 135,  25,    ,   0, SRS, 312,   6,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  12, 321N,  741W, 110,  948, HU,  50, NEQ,  100,   80,   60,   80, 1013,    0,  15, 135,  25,    ,   0, SRS, 312,   6,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  12, 321N,  741W, 110,  948, HU,  64, NEQ,   60,   60,   40,   60, 1013,    0,  15, 135,  25,    ,   0, SRS, 312,   6,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  24, 334N,  759W, 110,  948, HU,  34, NEQ,  170,  140,  100,  140, 1013,    0,  15, 135,  25,    ,   0, SRS, 311,   5,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  24, 334N,  759W, 110,  948, HU,  50, NEQ,  100,   80,   60,   80, 1013,    0,  15, 135,  25,    ,   0, SRS, 311,   5,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  24, 334N,  759W, 110,  948, HU,  64, NEQ,   60,   60,   40,   60, 1013,    0,  15, 135,  25,    ,   0, SRS, 311,   5,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  36, 339N,  771W, 105,  954, HU,  34, NEQ,  170,  140,  100,  140, 1013,    0,  15, 130,  25,    ,   0, SRS, 297,   3,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  36, 339N,  771W, 105,  954, HU,  50, NEQ,  100,   80,   60,   80, 1013,    0,  15, 130,  25,    ,   0, SRS, 297,   3,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  36, 339N,  771W, 105,  954, HU,  64, NEQ,   60,   60,   40,   60, 1013,    0,  15, 130,  25,    ,   0, SRS, 297,   3,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  48, 340N,  779W, 100,  959, HU,  34, NEQ,  170,  140,  100,  120, 1013,    0,  15, 120,  25,    ,   0, SRS, 279,   2,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  48, 340N,  779W, 100,  959, HU,  50, NEQ,  100,   80,   60,   60, 1013,    0,  15, 120,  25,    ,   0, SRS, 279,   2,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  48, 340N,  779W, 100,  959, HU,  64, NEQ,   60,   60,   40,   40, 1013,    0,  15, 120,  25,    ,   0, SRS, 279,   2,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  72, 336N,  792W,  70,  987, HU,  34, NEQ,  170,  140,   70,   70, 1013,    0,  15,  85,  25,    ,   0, SRS, 250,   1,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  72, 336N,  792W,  70,  987, HU,  50, NEQ,   80,   80,   40,   40, 1013,    0,  15,  85,  25,    ,   0, SRS, 250,   1,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL,  96, 340N,  817W,  30, 1008, TD,  34, NEQ,    0,    0,    0,    0, 1013,    0,  15,  40,  25,    ,   0, SRS, 282,   3,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL, 120, 356N,  834W,  20, 1011, LO,  34, NEQ,    0,    0,    0,    0, 1013,    0,  15,  30,  25,    ,   0, SRS, 319,   3,           ,,,,,,,,
AL, 06, 2018091218, 03, OFCL, 144, 391N,  820W,  15, 1012, EX,  34, NEQ,    0,    0,    0,    0, 1013,    0,  15,  20,  25,    ,   0, SRS,  17,   5,           ,,,,,,,,
FariborzDaneshvar-NOAA commented

Update!

@pvelissariou1 Here are PaHM outputs for v_wind for all four PaHM runs:

PaHM_florence_v_wind_max