Fb coupling fields

manual/impl/switch.tex

@@ -186,15 +186,15 @@ \subsubsection{~Mandatory switches} \label{sub:man_switch}
 \end{slist}
 
 \noindent
-Selection of method of wind interpolation (time):
+Selection of method of wind/momentum interpolation (time):
 \begin{slist}
 \sit{wnt0}{No interpolation.}
 \sit{wnt1}{Linear interpolation.}
 \sit{wnt2}{Approximately quadratic interpolation.}
 \end{slist}
 
 \noindent
-Selection of method of wind interpolation (space):
+Selection of method of wind/momentum interpolation (space):
 \begin{slist}
 \sit{wnx0}{Vector interpolation.}
 \sit{wnx1}{Approximately linear speed interpolation.}

model/bin/ww3_from_ftp.sh

@@ -28,7 +28,8 @@ 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"  
 cp -r data_regtests/ww3_tp2.18/input/*.nc  regtests/ww3_tp2.18/input/
-cp -r data_regtests/ww3_tp2.15/input/*.nc  regtests/ww3_tp2.15/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/
 cp -r data_regtests/ww3_tp2.13/*.png       regtests/ww3_tp2.13/
 cp -r data_regtests/ww3_tic1.4/input/*.nc  regtests/ww3_tic1.4/input/
 cp -r data_regtests/ww3_tp2.8/input/*.nc   regtests/ww3_tp2.8/input/

model/ftn/gx_outf.ftn

@@ -8,7 +8,7 @@
 !/                  |           A. Chawla               |
 !/                  |           J.H.G.M. Alves          |
 !/                  |                        FORTRAN 90 |
-!/                  | Last update :         27-Aug-2015 |
+!/                  | Last update :         22-Mar-2021 |
 !/                  +-----------------------------------+
 !/
 !/    30-Jun-1999 : Final FORTRAN 77                    ( version 1.18 )
@@ -37,6 +37,7 @@
 !/    25-Feb-2013 : Adding double-index output fields   ( version 4.11 ) 
 !/    27-Aug-2015 : ICEH and ICEF added as output       ( version 5.10 )
 !/    25-Aug-2018 : Add WBT parameter                   ( version 6.06 )
+!/    22-Mar-2021 : RHOA and TAUA added as output       ( version 7.xx )
 !/
 !/
 !/    Copyright 2009 National Weather Service (NWS),
@@ -124,12 +125,13 @@
       USE W3TIMEMD, ONLY: STME21, TICK21, DSEC21
 !/
       USE W3GDATMD
-      USE W3WDATMD, ONLY: TIME, WLV, ICE, ICEH, ICEF, BERG, UST, USTDIR
+      USE W3WDATMD, ONLY: TIME, WLV, ICE, ICEH, ICEF, BERG,            &
+                          UST, USTDIR, RHOAIR
       USE W3ADATMD, ONLY: CFLXYMAX, CFLTHMAX, AS, CX, CY, UA, UD, WN,  &
                           DW, HS, WLM, T01, T0M1, T02, THM, THS, FP0,  &
                           THP0, FP1, THP1, ABA, ABD, UBA, UBD, FCUT,   &
                           SXX, SYY, SXY, PHS, PTP, PLP, PDIR, PSI, PWS,&
-                          PTM1, PT1, PT2, PEP,                         &
+                          PTM1, PT1, PT2, PEP, TAUA, TAUADIR,          &
                           PTHP0, PQP, PSW, PPE, PGW, QP,               &
                           PWST, PNR, USERO, TAUOX, TAUOY, TAUWIX,      &
                           TAUWIY, PHIAW, PHIOC, TUSX, TUSY, PRMS, TPMS,&
@@ -734,7 +736,7 @@
 !/                  | WAVEWATCH III           NOAA/NCEP |
 !/                  |           H. L. Tolman            |
 !/                  |                        FORTRAN 90 |
-!/                  | Last update :         28-Mar-2007 |
+!/                  | Last update :         22-Mar-2021 |
 !/                  +-----------------------------------+
 !/
 !/    30-Jun-1999 : Final FORTRAN 77                    ( version 1.18 )
@@ -748,6 +750,7 @@
 !/    18-Jan-2007 : Update MSOUT/MBOUT treatment.       ( version 3.10 )
 !/    28-Mar-2007 : Adding partitioned output.          ( version 3.11 )
 !/                  Adding user slots for outputs.
+!/    22-Mar-2021 : Add output fields RHOA and TAUA     ( version 7.xx )
 !/
 !  1. Purpose :
 !
@@ -971,15 +974,30 @@
                 FLONE = .TRUE.
                 CALL W3S2XY ( NSEA, NSEA, NX, NY, ICE   , MAPSF, X1 )
 !
+!     Atmospheric momentum 
+!
+              ELSE IF ( J.EQ.1 .AND. K.EQ.8 ) THEN
+                FLTWO = .TRUE.
+                CALL W3S2XY ( NSEA, NSEA, NX, NY, TAUA(1:NSEA)        &
+                                                        , MAPSF, XX )
+                CALL W3S2XY ( NSEA, NSEA, NX, NY, TAUADIR(1:NSEA)     &
+                                                        , MAPSF, XY )
+!
+!     Air density
+!
+              ELSE IF ( J.EQ.1 .AND. K.EQ.9 ) THEN
+                FLONE = .TRUE.
+                CALL W3S2XY ( NSEA, NSEA, NX, NY, RHOAIR, MAPSF, X1 )
+!
 !     Ice thickness
 !
-!/IS2        ELSE IF ( J.EQ.1 .AND. K.EQ.9 ) THEN
+!/IS2        ELSE IF ( J.EQ.1 .AND. K.EQ.11 ) THEN
 !/IS2          FLONE = .TRUE.
 !/IS2          CALL W3S2XY (NSEA, NSEA, NX, NY, ICEH  , MAPSF, X1 )
 !
 !     Average sea ice floe diameter
 !
-!/IS2        ELSE IF ( J.EQ.1 .AND. K.EQ.10) THEN
+!/IS2        ELSE IF ( J.EQ.1 .AND. K.EQ.12) THEN
 !/IS2          FLONE = .TRUE.
 !/IS2          CALL W3S2XY (NSEA, NSEA, NX, NY, ICEF  , MAPSF, X1 )
 !

model/ftn/pdlib_field_vec.ftn

@@ -6,10 +6,11 @@ MODULE PDLIB_FIELD_VEC
 !/                  | Mathieu Dutour-Sikiric (IRB)      |
 !/                  |                                   |
 !/                  |                        FORTRAN 90 |
-!/                  | Last update :         01-Jan-2010 |
+!/                  | Last update :         22-Mar-2021 |
 !/                  +-----------------------------------+
 !/
 !/    01-Jan-2010 : Origination.                        ( version 6.04 )
+!/    22-Mar-2021 : Add WNMEAN, TAUOC output            ( version 7.xx )
 !/
 !/    Copyright 2010 National Weather Service (NWS),
 !/       National Oceanic and Atmospheric Administration.  All rights
@@ -163,6 +164,9 @@ MODULE PDLIB_FIELD_VEC
       IF ( FLGRDALL( 2, 17) ) THEN
         IH = IH + 1
       END IF
+      IF ( FLGRDALL( 2, 19) ) THEN
+        IH = IH + 1
+      END IF
       IF ( FLGRDALL( 3, 1) ) THEN
         DO IK=E3DF(2,1),E3DF(3,1)
           IH = IH + 1
@@ -317,6 +321,10 @@ MODULE PDLIB_FIELD_VEC
       IF ( FLGRDALL( 6, 11) ) THEN
         IH = IH + 1
       END IF
+      IF ( FLGRDALL( 6, 13) ) THEN
+        IH = IH + 1
+        IH = IH + 1
+      END IF
       IF ( FLGRDALL( 7, 1) ) THEN
         IH = IH + 1
         IH = IH + 1
@@ -837,10 +845,11 @@ MODULE PDLIB_FIELD_VEC
 !/                  | Mathieu Dutour-Sikiric (IRB)      |
 !/                  |                                   |
 !/                  |                        FORTRAN 90 |
-!/                  | Last update :         01-Mai-2018 |
+!/                  | Last update :         22-Mar-2021 |
 !/                  +-----------------------------------+
 !/
 !/    01-Mai-2018 : Origination.                        ( version 6.04 )
+!/    22-Mar-2021 : Add WNMEAN, TAUOC output            ( version 7.xx )
 !/
 !  1. Purpose : Do communication for PDLIB output
 !  2. Method :
@@ -863,7 +872,7 @@ MODULE PDLIB_FIELD_VEC
                           DTDYN, FCUT, SPPNT, ABA, ABD, UBA, UBD,&
                           SXX, SYY, SXY, USERO, PHS, PTP, PLP,   &
                           PDIR, PSI, PWS, PWST, PNR, PHIAW,      &
-                          PHIOC,                                 &
+                          PHIOC, TAUOCX, TAUOCY, WNMEAN,         &
                           TUSX, TUSY, TAUWIX, TAUWIY, TAUOX,     &
                           TAUOY, USSX, USSY, MSSX, MSSY,         &
                           MSCX, MSCY, PRMS, TPMS, CHARN,         &
@@ -1013,6 +1022,10 @@ MODULE PDLIB_FIELD_VEC
               IH = IH + 1
               Arrexch(IH,JSEA)=WBT(JSEA)
             END IF
+            IF ( FLGRDALL( 2, 19) ) THEN
+              IH = IH + 1
+              Arrexch(IH,JSEA)=WNMEAN(JSEA)
+            END IF
             IF ( FLGRDALL( 3, 1) ) THEN 
               DO IK=E3DF(2,1),E3DF(3,1)
                 IH = IH + 1
@@ -1251,6 +1264,12 @@ MODULE PDLIB_FIELD_VEC
               IH = IH + 1
               Arrexch(IH,JSEA)=PHICE(JSEA)
             END IF
+            IF ( FLGRDALL( 6, 13) ) THEN
+              IH = IH + 1
+              Arrexch(IH,JSEA)=TAUOCX(JSEA)
+              IH = IH + 1
+              Arrexch(IH,JSEA)=TAUOCY(JSEA)
+            END IF
             IF ( FLGRDALL( 7, 1) ) THEN
               IH = IH + 1
               Arrexch(IH,JSEA)=ABA(JSEA)
@@ -1474,6 +1493,10 @@ MODULE PDLIB_FIELD_VEC
                 IH = IH + 1
                 WBT(1:NSEA) = ARRtotal(IH,:)
               END IF
+              IF ( FLGRDALL( 2, 19) ) THEN
+                IH = IH + 1
+                WNMEAN(1:NSEA) = ARRtotal(IH,:)
+              END IF
               IF ( FLGRDALL( 3, 1) ) THEN 
                 DO IK=E3DF(2,1),E3DF(3,1)
                   IH = IH + 1
@@ -1712,6 +1735,12 @@ MODULE PDLIB_FIELD_VEC
                 IH = IH + 1
                 PHICE(1:NSEA) = ARRtotal(IH,:)
               END IF
+              IF ( FLGRDALL( 6, 13) ) THEN
+                IH = IH + 1
+                TAUOCX(1:NSEA) = ARRtotal(IH,:)
+                IH = IH + 1
+                TAUOCY(1:NSEA) = ARRtotal(IH,:)
+              END IF
               IF ( FLGRDALL( 7, 1) ) THEN
                 IH = IH + 1
                 ABA(1:NSEA) = ARRtotal(IH,:)