Fb wind stress

manual/impl/switch.tex

@@ -59,6 +59,7 @@ \subsubsection{~Mandatory switches} \label{sub:man_switch}
 \sit{flx2} {Friction velocity from Tolman and Chalikov input.}
 \sit{flx3} {Idem, with cap of Eq.~(\ref{eq:Cd_cap_1}) or (\ref{eq:Cd_cap_2}).}
 \sit{flx4} {Friction velocity according to Eq.~(\ref{eq:ST607}).}
+\sit{flx5} {Friction velocity from atmospheric stress (stresses should be an input).}
 \end{slist}
 
 \noindent

model/bin/build_utils.sh

@@ -99,7 +99,7 @@ check_switches()
 #sort:stress:
       stress ) TY='one'
                ID='stress computation'
-               OK='FLX0 FLX1 FLX2 FLX3 FLX4' ;;
+               OK='FLX0 FLX1 FLX2 FLX3 FLX4 FLX5' ;;
 #sort:dstress:
       dstress) TY='upto1'
                ID='Diagnostic stress comp'
@@ -444,6 +444,7 @@ check_switches()
    FLX2) str_st1='OK' ; str_st2='OK' ; str_st3='no' ; str_st6='OK' ;;
    FLX3) str_st1='OK' ; str_st2='OK' ; str_st3='no' ; str_st6='OK' ;;
    FLX4) str_st1='OK' ; str_st2='no' ; str_st3='no' ; str_st6='OK' ;;
+   FLX5) str_st1='no' ; str_st2='no' ; str_st3='yes' ; str_st6='OK' ;;
   esac
 
   if [ -n "$thread2" ] && [ "$thread1" != 'OMPG' ]
@@ -518,7 +519,7 @@ check_switches()
   then
       echo ' '
       echo "   *** !/ST6 cannot be used in combination with !/$stress"
-      echo "       Choose from FLX1, FLX2, FLX3, or FLX4."
+      echo "       Choose from FLX1, FLX2, FLX3, or FLX4 or FLX5."
       echo ' ' ; exit 14
   fi
 
@@ -592,6 +593,8 @@ switch_files()
          flxx=$NULL ;;
    FLX4) flx='w3flx4md'
          flxx=$NULL ;;
+   FLX5) flx='w3flx5md'
+         flxx=$NULL ;;
   esac
 
   case $dstress in
@@ -767,7 +770,7 @@ create_file_list()
                core=
                data='w3wdatmd w3gdatmd w3adatmd w3idatmd w3odatmd wmmdatmd'
                prop=
-            sourcet="w3parall w3triamd w3gridmd $stx $nlx $btx $is $uostmd"
+            sourcet="w3parall w3triamd w3gridmd $stx $flx $nlx $btx $is $uostmd"
                  IO='w3iogrmd'
                 aux="constants w3servmd w3arrymd w3dispmd w3gsrumd w3timemd w3nmlgridmd $pdlibyow $memcode"
                 if [ "$scrip" = 'SCRIP' ]
@@ -784,7 +787,7 @@ create_file_list()
                core=
                data="$memcode w3gdatmd w3wdatmd w3adatmd w3idatmd w3odatmd"
                prop=
-            sourcet="$pdlibcode $pdlibyow $db $tr $trx $bt $setupcode $stx $nlx $btx $is wmmdatmd w3parall $uostmd"
+            sourcet="$pdlibcode $pdlibyow $db $tr $trx $bt $setupcode $stx $flx $nlx $btx $is wmmdatmd w3parall $uostmd"
                  IO='w3iogrmd w3iorsmd'
                 aux="constants w3triamd w3servmd w3arrymd w3dispmd w3gsrumd w3timemd" ;;
       ww3_bound)
@@ -805,14 +808,14 @@ create_file_list()
                core='w3fldsmd'
                data="$memcode w3gdatmd w3adatmd w3idatmd w3odatmd w3wdatmd wmmdatmd"
                prop=
-            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode w3triamd $stx $nlx $btx  $is $uostmd"
+            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode w3triamd $stx $flx $nlx $btx  $is $uostmd"
                  IO="w3iogrmd $oasismd $agcmmd $ogcmmd $igcmmd"
                 aux="constants w3servmd w3timemd $tidecode w3arrymd w3dispmd w3gsrumd w3parall" ;;
       ww3_prnc)
                core='w3fldsmd'
                data="$memcode w3gdatmd w3adatmd w3idatmd w3odatmd w3wdatmd wmmdatmd"
                prop=
-            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode w3triamd $stx $nlx $btx $is w3parall $uostmd"
+            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode w3triamd $stx $flx $nlx $btx $is w3parall $uostmd"
                  IO="w3iogrmd $oasismd $agcmmd $ogcmmd $igcmmd"
                 aux="constants w3servmd w3timemd w3arrymd w3dispmd w3gsrumd w3tidemd w3nmlprncmd" ;;
       ww3_prtide) 
@@ -882,15 +885,15 @@ create_file_list()
                core=
                data="w3parall wmmdatmd $memcode w3gdatmd w3wdatmd w3adatmd w3idatmd w3odatmd"
                prop=
-            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode $tr $trx $stx $nlx $btx $is $uostmd"
+            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode $tr $trx $stx $flx $nlx $btx $is $uostmd"
                  IO='w3iogrmd w3iogomd'
                 aux="constants w3servmd w3timemd w3arrymd w3dispmd w3gsrumd"
                 aux="$aux" ;;
       ww3_ounf)
                core='w3initmd'
                data="wmmdatmd $memcode w3gdatmd w3wdatmd w3adatmd w3idatmd w3odatmd"
                prop=
-            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode w3parall w3triamd $stx $nlx $btx  $is $uostmd"
+            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode w3parall w3triamd $stx $flx $nlx $btx  $is $uostmd"
                  IO='w3iogrmd w3iogomd w3iorsmd w3iopomd'
                 aux="constants w3servmd w3timemd w3arrymd w3dispmd w3gsrumd"
                 aux="$aux w3nmlounfmd $smco w3ounfmetamd w3metamd" ;;
@@ -927,15 +930,15 @@ create_file_list()
                core=
                data="w3parall wmmdatmd w3triamd $memcode w3gdatmd w3wdatmd w3adatmd w3idatmd w3odatmd"
                prop=
-            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode $stx $nlx $btx $is $uostmd"
+            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode $stx $flx $nlx $btx $is $uostmd"
                  IO='w3iogrmd w3iogomd'
                 aux="constants w3servmd w3timemd w3arrymd w3dispmd w3gsrumd"
                 aux="$aux" ;;
       ww3_gspl) 
                core='w3fldsmd'
                data="$memcode w3gdatmd w3wdatmd w3adatmd w3idatmd w3odatmd"
                prop=
-            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode wmmdatmd w3parall w3triamd $stx $nlx $btx $is $uostmd"
+            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode wmmdatmd w3parall w3triamd $stx $flx $nlx $btx $is $uostmd"
                  IO="w3iogrmd  $oasismd $agcmmd $ogcmmd $igcmmd"
                 aux="constants w3servmd w3timemd w3arrymd w3dispmd w3gsrumd $tidecode" ;;
       ww3_gint) 
@@ -950,7 +953,7 @@ create_file_list()
                core=
                data="$memcode w3gdatmd w3wdatmd w3adatmd w3idatmd w3odatmd"
                prop=
-            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode wmmdatmd w3parall w3triamd $stx $nlx $btx $tr $bs $is $uostmd"
+            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode wmmdatmd w3parall w3triamd $stx $flx $nlx $btx $tr $bs $is $uostmd"
                  IO='w3iogrmd w3iogomd'
                 aux="constants w3servmd w3timemd w3arrymd w3dispmd w3gsrumd"
                 aux="$aux" ;;
@@ -965,7 +968,7 @@ create_file_list()
                core='w3strkmd'
                data="$memcode w3gdatmd w3adatmd w3idatmd w3odatmd w3wdatmd"
                prop=
-            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode wmmdatmd w3dispmd w3triamd $ln $stx $nlx $btx $tr $bs $is $uostmd"
+            sourcet="$pdlibcode $pdlibyow $db $bt $setupcode wmmdatmd w3dispmd w3triamd $ln $stx $flx $nlx $btx $tr $bs $is $uostmd"
                  IO=
                 aux="constants w3servmd w3timemd w3arrymd w3gsrumd w3parall" ;;
       libww3|libww3.so) 

model/bin/make_makefile.sh

@@ -114,7 +114,7 @@
 
 
   progs="ww3_grid ww3_strt ww3_prep ww3_prnc ww3_shel ww3_multi ww3_sbs1
-         ww3_outf ww3_outp ww3_trck ww3_trnc ww3_grib gx_outf gx_outp ww3_ounf 
+         ww3_outf ww3_outp ww3_trck ww3_trnc ww3_grib gx_outf gx_outp ww3_ounf
          ww3_ounp ww3_gspl ww3_gint ww3_bound ww3_bounc ww3_systrk $tideprog"
   progs="$progs ww3_multi_esmf  ww3_uprstr"
   progs="$progs libww3"
@@ -220,7 +220,7 @@
       done
       echo "	@cd \$(aPo); ld -o $lib -shared $objs" >> makefile
       echo ' '                                   >> makefile
-      
+
     else
       echo '	@$(aPb)/link '"$filesl"          >> makefile
       echo ' '                                   >> makefile
@@ -250,7 +250,7 @@
 
     suffixes="ftn f F f90 F90 c"
     fexti=none
-    ispdlibi=no 
+    ispdlibi=no
     for s in $suffixes
     do
       if [ -f $main_dir/ftn/$file.$s ]
@@ -285,9 +285,9 @@
     string3="$NULL"
 
     if [ "$ispdlibi" = 'yes' ]
-    then 
+    then
       string1='$(aPo)/'$file'.o : PDLIB/'$file.$fexti' '
-    fi 
+    fi
 
     $main_dir/bin/ad3 $file 0 1 > ad3.out 2>&1
 
@@ -351,6 +351,7 @@
          'W3FLX2MD'     ) modtest=w3flx2md.o ;;
          'W3FLX3MD'     ) modtest=w3flx3md.o ;;
          'W3FLX4MD'     ) modtest=w3flx4md.o ;;
+         'W3FLX5MD'     ) modtest=w3flx5md.o ;;
          'W3SLN1MD'     ) modtest=w3sln1md.o ;;
          'W3SRC0MD'     ) modtest=w3src0md.o ;;
          'W3SRC1MD'     ) modtest=w3src1md.o ;;
@@ -435,7 +436,7 @@
          'W3OUNFMETAMD' ) modtest=w3ounfmetamd.o ;;
          'W3METAMD'     ) modtest=w3metamd.o ;;
          'W3GRIDMD'     ) modtest=w3gridmd.o ;;
-         *              ) modfound=no ;; 
+         *              ) modfound=no ;;
       esac
 
       if [ "$modfound" == "yes" ]
@@ -494,9 +495,9 @@
     fi
 
     if  [ "$scripnc" = 'SCRIPNC' ]
-    then 
+    then
        scrip_mk=$scrip_dir/SCRIP_NC.mk
-    else 
+    else
        scrip_mk=$scrip_dir/SCRIP.mk
     fi
     if [ ! -e $scrip_mk ]

model/bin/ww3_from_ftp.sh

@@ -10,23 +10,23 @@ curr_dir=`pwd`
 # Set WW3 code version
 ww3ver=v7.12.6
 
-#Get top level directory of ww3 from user: 
+#Get top level directory of ww3 from user:
 echo -e "\n\n This script will download data from the ftp for WAVEWATCH III "
 echo -e "Enter the relative path to the main/top level directory, this would "
 echo -e "be '../../' if in the model/bin directory or '.' if already in the "
 echo -e "top/main directory:"
-read ww3dir 
+read ww3dir
 
-#Move to top level directory of ww3: 
-cd $ww3dir 
+#Move to top level directory of ww3:
+cd $ww3dir
 
-#Download from ftp and uptar: 
-echo -e "Downloading and untaring file from ftp:" 
+#Download from ftp and uptar:
+echo -e "Downloading and untaring file from ftp:"
 wget https://ftp.emc.ncep.noaa.gov/static_files/public/WW3/ww3_from_ftp.${ww3ver}.tar.gz
 tar -xvzf ww3_from_ftp.${ww3ver}.tar.gz
 
 #Move regtest info from data_regtests to regtests:
-echo -e "Moving data from data_regtests to regtests"  
+echo -e "Moving data from data_regtests to regtests"
 cp -r data_regtests/ww3_tp2.18/input/*.nc  regtests/ww3_tp2.18/input/
 cp -r data_regtests/ww3_tp2.15/input/wind.nc  regtests/ww3_tp2.15/input/
 cp -r data_regtests/ww3_tp2.15/input/*.nc  regtests/ww3_tp2.15/input_rho/
@@ -43,7 +43,7 @@ if [ ! -d regtests/ww3_tp2.14/input/oasis3-mct/doc ]
 then
   mkdir regtests/ww3_tp2.14/input/oasis3-mct/doc
 fi
-cp -r data_regtests/ww3_tp2.14/input/oasis3-mct/doc/* regtests/ww3_tp2.14/input/oasis3-mct/doc/ 
+cp -r data_regtests/ww3_tp2.14/input/oasis3-mct/doc/* regtests/ww3_tp2.14/input/oasis3-mct/doc/
 cp -r data_regtests/ww3_tp2.14/input/toy/*.nc.OAS*CM regtests/ww3_tp2.14/input/toy/
 cp -r data_regtests/ww3_tp2.14/input/toy/r-toy.nc.OASACM regtests/ww3_tp2.14/input/toy/r-toy.nc.OASACM2
 cp -r data_regtests/ww3_tp2.14/input/toy/toy_coupled_field.nc.OASACM regtests/ww3_tp2.14/input/toy/toy_coupled_field.nc.OASACM2
@@ -59,15 +59,15 @@ cp -r data_regtests/ww3_ufs1.1/input/*.nc  regtests/ww3_ufs1.2/input/
 cp -r data_regtests/ww3_ufs1.2/input/*     regtests/ww3_ufs1.2/input/
 cp -r data_regtests/ww3_ufs1.3/input/*nc         regtests/ww3_ufs1.3/input/
 
-#Do you want to clean up (aka delete tar file, delete the data_regtests directory) 
+#Do you want to clean up (aka delete tar file, delete the data_regtests directory)
 echo -e "\n\n Do you want to delete the tar file ww3_from_ftp.${ww3ver}.tar.gz [y|n]: "
 read wnew
 if [ "${wnew}" = "Y" ] || [ "${wnew}" = "y" ]
 then
   echo -e '\n Deleting tar file ww3_from_ftp.${ww3ver}.tar.gz'
   rm ww3_from_ftp.${ww3ver}.tar.gz
 else
-  echo -e ' Not deleting tar file.' 
+  echo -e ' Not deleting tar file.'
 fi
 
 echo -e "\n\n Files were copied from the data_regtests to the regtests folder."
@@ -78,9 +78,9 @@ then
   echo -e '\n Deleting the data_regtests folder'
   rm -rf data_regtests
 else
-  echo -e ' Not deleting the data_regtests folder.' 
+  echo -e ' Not deleting the data_regtests folder.'
 fi
 
-#move back to original directory: 
+#move back to original directory:
 cd $curr_dir
 echo -e "End of ww3_from_ftp.sh"