Grid configurations consistent with new run_model.sh at FCOO

CarstenHansen · Oct 18, 2024 · 561058b · 561058b
1 parent 44a466f
commit 561058b
Show file tree

Hide file tree

Showing 15 changed files with 243 additions and 11 deletions.
diff --git a/fcoo/ARCT18/initdir/run_model.conf b/fcoo/ARCT18/initdir/run_model.conf
@@ -18,7 +18,9 @@
 # SIMLEN=198000
 # Grid data and model output data directories
   MOD_DEF=$RUNDIRS/mod_defs
-  NESTDIR=$DATA/nesting/$mod_ID
+# ARCT18 has no 'NESTFROM', so NESTDIR is blank
+  NESTDIR=
+  NESTTO=$DATA/nesting
   HOTDIR=$DATA/hotstart/$mod_ID
   NETCDF=$DATA/netcdf/$mod_ID    
 # Ice and wind file names in ww3_prnc.nml:                    

diff --git a/fcoo/GRL9/initdir/points.GRL9 b/fcoo/GRL9/initdir/points.GRL9
@@ -0,0 +1,24 @@
+$
+$ Added 20210101:
+   -23.110   65.253  'Flateyjard' $ Flateyjardufl
+$ 9 Iceland stations
+   -20.345   63.286  'Surtseyjar' $ Surtseyjardufl
+   -22.461   63.813  'Grindaviku' $ Grindavikurdufl
+   -22.877   64.052  'Gardskagad' $ Gardskagadufl
+   -15.189   64.197  'Hornafjard' $ Hornafjardardufl
+   -13.627   65.648  'Kogurdufl ' $ Kogurdufl
+   -24.778   65.698  'Blakksnes ' $ Blakksnes
+   -21.138   65.770  'Drangsnes ' $ Drangsnes
+   -18.192   66.292  'Grimseyjar' $ Grimseyjarsund
+   -23.367   66.438  'Straumnesd' $ Straumnesdufl
+$ Virtual Navigation stations
+   -37.0     60.3    'DanmaStr01'
+   -31.23    60.31   'DanmaStr02'
+   -44.0     59.5    'CapeFare_1'
+$ 12 nmi off Greenland Harbours 
+   -50.2     59.85   'Paamiut_1 '
+   -52.5     63.95   'Nuuk_1    '
+   -54.5     66.9    'Sisimiut_1'
+   -54.3     68.6    'Aasiaat_1 '
+   -37.6     65.25   'Tasiilaq_1'
+$
diff --git a/fcoo/GRL9/initdir/run_model.conf b/fcoo/GRL9/initdir/run_model.conf
@@ -0,0 +1,28 @@
+# mod_ID=ARCT18 model configuration to be parsed in run_model.sh
+#
+# Model set up directories
+  WW3DIR=/home/$USER/Models/WW3_git/WW3
+  MODDIR=$WW3DIR/fcoo/GRL9
+# input directories relative to MODDIR
+  INPDIR=initdir
+  GRIDDIR=setup
+# Model run base directories
+  RUNDIRS=/home/$USER/Rundirs/WW3RUN
+  # DATA=/scratch/${USER}_today/WW3
+  DATA=/home/$USER/Datadirs/WW3RUN
+# Model run datetimes
+  HOTDELTA=21600                         
+  FIELDDELTA=3600                          
+  POINTDELTA=1200                          
+  BOUNDDELTA=1200
+# SIMLEN=198000
+# Grid data and model output data directories
+  MOD_DEF=$RUNDIRS/mod_defs
+  NESTDIR=$DATA/nesting/$mod_ID
+  HOTDIR=$DATA/hotstart/$mod_ID
+  NETCDF=$DATA/netcdf/$mod_ID    
+# Ice and wind file names in ww3_prnc.nml:                    
+  prnc_ice=icec.nc
+  prnc_wind=wind.nc
+# Nest-to inner models
+  NEST1=NUUP_hires
diff --git a/fcoo/GRL9/initdir/ww3_ounf_GRL9.inp.tpl b/fcoo/GRL9/initdir/ww3_ounf_GRL9.inp.tpl
@@ -4,7 +4,7 @@ $--------------------------------------------------------------------- $
 $ First output time (yyyymmdd hhmmss), increment of output (s), 
 $ and (max) number of output times.
 $
-    __SIMSTART__  3600  721
+    __SIMSTART__   __FIELDDELTA__  4470
 $  19850101 000000  3600.  1000
 $
 $ Fields requested --------------------------------------------------- $

diff --git a/fcoo/GRL9/initdir/ww3_ounp_GRL9.inp.tpl b/fcoo/GRL9/initdir/ww3_ounp_GRL9.inp.tpl
@@ -5,7 +5,7 @@ $ First output time (yyyymmdd hhmmss), increment of output (s),
 $ and (max) number of output times.
 $
 $  19850101 000000  3600.  1000
-    __SIMSTART__  3600.  721
+    __SIMSTART__  __POINTDELTA__  4470
 $
 $ Points requested --------------------------------------------------- $
 $

diff --git a/fcoo/GRL9/initdir/ww3_shel_GRL9.nml.tpl b/fcoo/GRL9/initdir/ww3_shel_GRL9.nml.tpl
@@ -297,16 +297,16 @@
 ! -------------------------------------------------------------------- !
 &OUTPUT_DATE_NML
   DATE%FIELD%START         =  '__SIMSTART__'
-  DATE%FIELD%STRIDE        =  3600
+  DATE%FIELD%STRIDE        =  '__FIELDDELTA__'
   DATE%FIELD%STOP          =  '__SIMSTOP__'
   DATE%POINT%START         =  '__SIMSTART__'
-  DATE%POINT%STRIDE        =  1200
+  DATE%POINT%STRIDE        =  '__POINTDELTA__'
   DATE%POINT%STOP          =  '__SIMSTOP__'
   DATE%RESTART%START       =  '__SIMSTART__'
   DATE%RESTART%STRIDE      =  '__HOTDELTA__'
   DATE%RESTART%STOP        =  '__SIMSTOP__'
   DATE%BOUNDARY%START      =  '__SIMSTART__'
-  DATE%BOUNDARY%STRIDE     =  1200
+  DATE%BOUNDARY%STRIDE     =  '__BOUNDDELTA__'
   DATE%BOUNDARY%STOP       =  '__SIMSTOP__'
 /
 !  DATE%TRACK%START         =  '__SIMSTART__'

diff --git a/fcoo/GRL9/setup/namelists_GRL9_IC5_VE3_IS2.nml b/fcoo/GRL9/setup/namelists_GRL9_IC5_VE3_IS2.nml
@@ -9,5 +9,7 @@
        SDSSTRAIN =0., SDSSTRAIN2 = 0., FXFM3 = 40., FXPM3 = 20., SDSFACMTF = 400.,
        SDSCUM=0., SDSC5 =0. /
 &ROTD PLAT = 10.0, PLON = -125.0, UNROT = .FALSE. /
+! Boundary spectra to rotated grids: No. 1 in ww3_grid.nml: NUUP_hires, res 5400m
+&ROTB BPLAT(1) =  22.2959, BPLON(1) = 160.86456, /
 &PRO3 WDTHCG=2.0, WDTHTH=2.0 /
 END OF NAMELISTS
diff --git a/fcoo/GRL9/setup/nc2dat b/fcoo/GRL9/setup/nc2dat
@@ -0,0 +1 @@
+../../bin/nc2dat.py
diff --git a/fcoo/GRL9/setup/ww3_grid_GRL9.nml b/fcoo/GRL9/setup/ww3_grid_GRL9.nml
@@ -382,6 +382,20 @@ INBND_POINT(11) =     94 186 T  !
 !     OUTBND_LINE(I)%NP     = 0               ! number of points
 ! OR
 !     OUTBND_LINE(I)        = 0. 0. 0. 0. 0   ! included lines
+&OUTBND_COUNT_NML
+  OUTBND_COUNT%N_LINE   =  3
+/  
+&OUTBND_LINE_NML
+! nest1.ww3 to NUUP_hires grid (5400mX5400m)
+! The north pole of the rotated grid NUUP_5400 is specified in namelists.nml:
+!    &ROTB BPLAT(1) = 22.2959, BPLON(1) = 160.86456, /
+! The OUTBND_LINEs are specified in destination grid rotlon rotlat coordinates:
+  OUTBND_LINE(1)  =    -13.77658    -0.87531    -0.04833 0.00000   15
+  OUTBND_LINE(2)  =    -14.45320    -0.82698     0.00000 0.04833   38
+  OUTBND_LINE(3)  =    -14.40487     1.00956     0.04833 0.00000   23
+! nest2.ww3 etc. to other grids do not exist, yet.
+/
+! -------------------------------------------------------------------- !
 ! -------------------------------------------------------------------- !
 ! WAVEWATCH III - end of namelist                                      !
 ! -------------------------------------------------------------------- !
diff --git a/fcoo/NEA9/initdir/points.NEA9 b/fcoo/NEA9/initdir/points.NEA9
@@ -0,0 +1,79 @@
+$
+$ Added 20200812:
+    -6.100    50.103  '62107     ' $ 62107
+    -9.999    55.000  '62093     ' $ 62093
+   -11.401    59.100  '64045     ' $ 64045
+    -5.430    53.484  '62091     ' $ 62091
+    -9.268    53.228  'SmartBayWa' $ SmartBayWaveBuoy
+    -4.072    48.728  '6200053   ' $ 6200053
+    -4.968    48.290  '6200069   ' $ 6200069
+    -2.900    49.900  '62103     ' $ 62103
+    -2.889    49.026  '6200058   ' $ 6200058
+    -4.500    60.701  '64046     ' $ K7
+     1.000    61.100  '63112     ' $ 63112
+$ In CMEMS 202004:
+   -10.146   54.231  'AMETSBer-B' $ BelmulletB-coast-buoy
+   -10.278   54.268  'AMETSBer-A' $ BelmulletA-coast-buoy
+     1.500   59.500  '63110     ' $ 63110
+     2.227   59.574  'Heimdal   ' $ Heimdal
+     3.719   60.644  'Troll-A   ' $ Troll-A
+     3.817   55.417  'A121      ' $ A121, A122
+   -10.548   51.216  '62092     ' $ M3   
+    -6.704   51.690  '62094     ' $ M5   
+     2.269   61.204  'GullfaksC ' $ Gullfaks-C
+     1.909   58.371  'Sleipner-A' $ Sleipner-A
+$ Older, may be in inner domains:
+$     2.000   51.240  '62170     ' $ 62170
+$     2.703   51.350  'Kwintebank' $ Kwintebank
+$     2.439   51.389  'Westhinder' $ Westhinder
+$     3.370   51.393  'CadzanBoei' $ CadzandBoei
+$     2.770   51.410  'Akkaert   ' $ Akkaert
+$     3.415   51.433  'WielinNoor' $ WielingenNoord
+$     3.242   51.504  'Deurlo    ' $ Deurlo
+$     3.399   51.620  'DomburRass' $ DomburgerRassen
+$     3.312   51.748  'Schouwenba' $ Schouwenbank, Schouwenbank2
+$     3.617   51.769  'Brouwersha' $ Brouwershavensegat
+$     3.670   51.926  'LichteiGoe' $ LichteilandGoeree1, LichteilandGoeree2
+$     3.000   51.948  'EurogeuDWE' $ EurogeulDWE
+$     3.276   51.999  'Europlatf2' $ Europlatform2, Europlatform3
+$     3.737   52.010  'EurogeuE13' $ EurogeulE13
+$     4.268   52.493  'IJgeul5   ' $ IJgeul5
+$     4.058   52.550  'IJmuidMuni' $ IJmuidenMunitiestort, IJmuidenMunitiestort2
+     1.708   61.000  '63113     ' $ 63113
+    -2.295   46.833  '62067     ' $ 62067
+    -8.500   47.500  '62163     ' $ 62163
+     4.151   52.926  'Q1        ' $ Q1, Q11
+     1.700   53.000  '62130     ' $ 62130
+     2.800   53.103  '62145     ' $ 62145
+     3.220   53.218  'K13a      ' $ K13a, K13a2, K13a3
+     3.633   53.267  'K141      ' $ K141
+     4.662   53.277  'WaddenEGat' $ WaddenEierlandseGat
+     5.102   53.322  'StortemelO' $ StortemelkOost
+     4.989   53.323  'StortemelB' $ StortemelkBoei
+     1.700   53.400  '62144     ' $ 62144
+     6.203   53.470  'SchierWagg' $ SchiermonnikoogWaggen
+     5.574   53.478  'Amelander2' $ Amelander21, Amelander22
+     5.950   53.500  'AWG       ' $ AWG
+     5.482   53.509  'Amelander1' $ Amelander11, Amelander12
+     6.165   53.529  'SchierWest' $ SchiermonnikoogWestgat
+     6.167   53.596  'SchierNoor' $ SchiermonnikoogNoord
+     0.700   53.600  '62150     ' $ 62150
+     4.961   53.614  'L91       ' $ L91
+     6.518   53.617  'MeetboWEO1' $ MeetboeiWEO1
+     6.368   53.620  'MeetboWEW1' $ MeetboeiWEW1
+     1.100   53.700  '62149     ' $ 62149
+     2.800   53.800  '62146     ' $ 62146
+     2.950   53.817  'J61       ' $ J61
+     6.839   53.985  'NO1       ' $ NO1
+     0.700   54.000  '62127     ' $ 62127
+     1.100   54.000  '62165     ' $ 62165
+     8.114   54.017  'Elbe      ' $ ElbeWR
+     4.017   54.117  'F161      ' $ F161
+     7.891   54.180  'Helgoland ' $ HelgolandWR
+     7.818   54.219  'HelgolNort' $ Helgoland-NorthWR
+     4.727   54.854  'F3platform' $ F3platform
+     8.089   54.919  'Sylt      ' $ Sylt
+$    11.223   58.251  'Brofjorden' $ BrofjordenWR
+$    10.932   58.488  'Vaderoarna' $ VaderoarnaWR
+     2.826   60.491  'Oseberg-A ' $ Oseberg-A
+     1.149   61.240  'NorthCormo' $ NorthCormorant
diff --git a/fcoo/NEA9/initdir/ww3_ounf_NEA9.inp.tpl b/fcoo/NEA9/initdir/ww3_ounf_NEA9.inp.tpl
@@ -4,7 +4,7 @@ $--------------------------------------------------------------------- $
 $ First output time (yyyymmdd hhmmss), increment of output (s), 
 $ and (max) number of output times.
 $
-    __SIMSTART__  3600  721
+    __SIMSTART__   __FIELDDELTA__  4470
 $  19850101 000000  3600.  1000
 $
 $ Fields requested --------------------------------------------------- $

diff --git a/fcoo/NEA9/initdir/ww3_ounp_NEA9.inp.tpl b/fcoo/NEA9/initdir/ww3_ounp_NEA9.inp.tpl
@@ -5,7 +5,7 @@ $ First output time (yyyymmdd hhmmss), increment of output (s),
 $ and (max) number of output times.
 $
 $  19850101 000000  3600.  1000
-    __SIMSTART__  3600.  721
+    __SIMSTART__  __POINTDELTA__  4470
 $
 $ Points requested --------------------------------------------------- $
 $

diff --git a/fcoo/NEA9/initdir/ww3_shel_NEA9.nml.tpl b/fcoo/NEA9/initdir/ww3_shel_NEA9.nml.tpl
@@ -290,16 +290,16 @@
 ! -------------------------------------------------------------------- !
 &OUTPUT_DATE_NML
   DATE%FIELD%START         =  '__SIMSTART__'
-  DATE%FIELD%STRIDE        =  '3600'
+  DATE%FIELD%STRIDE        =  '__FIELDDELTA__'
   DATE%FIELD%STOP          =  '__SIMSTOP__'
   DATE%POINT%START         =  '__SIMSTART__'
-  DATE%POINT%STRIDE        =  '1200'
+  DATE%POINT%STRIDE        =  '__POINTDELTA__'
   DATE%POINT%STOP          =  '__SIMSTOP__'
   DATE%RESTART%START       =  '__SIMSTART__'
   DATE%RESTART%STRIDE      =  '__HOTDELTA__'
   DATE%RESTART%STOP        =  '__SIMSTOP__'
   DATE%BOUNDARY%START      =  '__SIMSTART__'
-  DATE%BOUNDARY%STRIDE     =  '1200'
+  DATE%BOUNDARY%STRIDE     =  '__BOUNDDELTA__'
   DATE%BOUNDARY%STOP       =  '__SIMSTOP__'
 /
 

diff --git a/fcoo/NWS9/initdir/points.NWS9 b/fcoo/NWS9/initdir/points.NWS9
@@ -0,0 +1,79 @@
+$
+$ Added 20200812:
+    -6.100    50.103  '62107     ' $ 62107
+    -9.999    55.000  '62093     ' $ 62093
+   -11.401    59.100  '64045     ' $ 64045
+    -5.430    53.484  '62091     ' $ 62091
+    -9.268    53.228  'SmartBayWa' $ SmartBayWaveBuoy
+    -4.072    48.728  '6200053   ' $ 6200053
+    -4.968    48.290  '6200069   ' $ 6200069
+    -2.900    49.900  '62103     ' $ 62103
+    -2.889    49.026  '6200058   ' $ 6200058
+    -4.500    60.701  '64046     ' $ K7
+     1.000    61.100  '63112     ' $ 63112
+$ In CMEMS 202004:
+   -10.146   54.231  'AMETSBer-B' $ BelmulletB-coast-buoy
+   -10.278   54.268  'AMETSBer-A' $ BelmulletA-coast-buoy
+     1.500   59.500  '63110     ' $ 63110
+     2.227   59.574  'Heimdal   ' $ Heimdal
+     3.719   60.644  'Troll-A   ' $ Troll-A
+     3.817   55.417  'A121      ' $ A121, A122
+   -10.548   51.216  '62092     ' $ M3   
+    -6.704   51.690  '62094     ' $ M5   
+     2.269   61.204  'GullfaksC ' $ Gullfaks-C
+     1.909   58.371  'Sleipner-A' $ Sleipner-A
+$ Older, may be in inner domains:
+$     2.000   51.240  '62170     ' $ 62170
+$     2.703   51.350  'Kwintebank' $ Kwintebank
+$     2.439   51.389  'Westhinder' $ Westhinder
+$     3.370   51.393  'CadzanBoei' $ CadzandBoei
+$     2.770   51.410  'Akkaert   ' $ Akkaert
+$     3.415   51.433  'WielinNoor' $ WielingenNoord
+$     3.242   51.504  'Deurlo    ' $ Deurlo
+$     3.399   51.620  'DomburRass' $ DomburgerRassen
+$     3.312   51.748  'Schouwenba' $ Schouwenbank, Schouwenbank2
+$     3.617   51.769  'Brouwersha' $ Brouwershavensegat
+$     3.670   51.926  'LichteiGoe' $ LichteilandGoeree1, LichteilandGoeree2
+$     3.000   51.948  'EurogeuDWE' $ EurogeulDWE
+$     3.276   51.999  'Europlatf2' $ Europlatform2, Europlatform3
+$     3.737   52.010  'EurogeuE13' $ EurogeulE13
+$     4.268   52.493  'IJgeul5   ' $ IJgeul5
+$     4.058   52.550  'IJmuidMuni' $ IJmuidenMunitiestort, IJmuidenMunitiestort2
+     1.708   61.000  '63113     ' $ 63113
+    -2.295   46.833  '62067     ' $ 62067
+    -8.500   47.500  '62163     ' $ 62163
+     4.151   52.926  'Q1        ' $ Q1, Q11
+     1.700   53.000  '62130     ' $ 62130
+     2.800   53.103  '62145     ' $ 62145
+     3.220   53.218  'K13a      ' $ K13a, K13a2, K13a3
+     3.633   53.267  'K141      ' $ K141
+     4.662   53.277  'WaddenEGat' $ WaddenEierlandseGat
+     5.102   53.322  'StortemelO' $ StortemelkOost
+     4.989   53.323  'StortemelB' $ StortemelkBoei
+     1.700   53.400  '62144     ' $ 62144
+     6.203   53.470  'SchierWagg' $ SchiermonnikoogWaggen
+     5.574   53.478  'Amelander2' $ Amelander21, Amelander22
+     5.950   53.500  'AWG       ' $ AWG
+     5.482   53.509  'Amelander1' $ Amelander11, Amelander12
+     6.165   53.529  'SchierWest' $ SchiermonnikoogWestgat
+     6.167   53.596  'SchierNoor' $ SchiermonnikoogNoord
+     0.700   53.600  '62150     ' $ 62150
+     4.961   53.614  'L91       ' $ L91
+     6.518   53.617  'MeetboWEO1' $ MeetboeiWEO1
+     6.368   53.620  'MeetboWEW1' $ MeetboeiWEW1
+     1.100   53.700  '62149     ' $ 62149
+     2.800   53.800  '62146     ' $ 62146
+     2.950   53.817  'J61       ' $ J61
+     6.839   53.985  'NO1       ' $ NO1
+     0.700   54.000  '62127     ' $ 62127
+     1.100   54.000  '62165     ' $ 62165
+     8.114   54.017  'Elbe      ' $ ElbeWR
+     4.017   54.117  'F161      ' $ F161
+     7.891   54.180  'Helgoland ' $ HelgolandWR
+     7.818   54.219  'HelgolNort' $ Helgoland-NorthWR
+     4.727   54.854  'F3platform' $ F3platform
+     8.089   54.919  'Sylt      ' $ Sylt
+$    11.223   58.251  'Brofjorden' $ BrofjordenWR
+$    10.932   58.488  'Vaderoarna' $ VaderoarnaWR
+     2.826   60.491  'Oseberg-A ' $ Oseberg-A
+     1.149   61.240  'NorthCormo' $ NorthCormorant
diff --git a/model/src/w3xsmfmd.F90 b/model/src/w3xsmfmd.F90
@@ -1099,6 +1099,9 @@ SUBROUTINE fit_expr(numD,log_KZ,Um,AM,PM,Um_min,fit_kind,P0,A0)
       ! The log of log_Um is linear in log_KZ:
       if ( present(A0) ) AM = A0 ! Best fit of PM, only, given AM=A0
       if ( present(P0) ) PM = P0 ! Best fit of log_AM, only, given PM=P0
+
+      ! Return if no depth
+      if ( numD == 0 ) return
 
       logAM0 = log(AM)
       loglog_ref(2:SPND) = logAM0 + PM * log_KZ(2:SPND)