diff -rupN WW3_314_13/ftn/w3updtmd.ftn WW3_314_14/ftn/w3updtmd.ftn --- WW3_314_13/ftn/w3updtmd.ftn 2012-10-02 21:25:15.000000000 -0400 +++ WW3_314_14/ftn/w3updtmd.ftn 2012-10-02 21:17:11.000000000 -0400 @@ -5,7 +5,7 @@ !/ | WAVEWATCH III NOAA/NCEP | !/ | H. L. Tolman | !/ | FORTRAN 90 | -!/ | Last update : 17-Aug-2010 | +!/ | Last update : 13-Mar-2012 | !/ +-----------------------------------+ !/ !/ 21-Jan-2000 : Origination. ( version 2.00 ) @@ -23,9 +23,11 @@ !/ 11-Jan-2007 : Clean-up W3UTRN boundary points. ( version 3.10 ) !/ 11-May-2007 : Adding NTPROC/NAPROC separation. ( version 3.11 ) !/ 29-May-2009 : Preparing distribution version. ( version 3.14 ) -!/ 17-Aug-2010 : ABPI0-N(:,0) init. bug fix. ( version 3.14.5 ) +!/ 17-Aug-2010 : ABPI0-N(:,0) init. bug fix. ( version 3.14 ) +!/ 13-Mar-2012 : Add initialization of UST on re- ( version 3.14 ) +!/ activation of grid point. !/ -!/ Copyright 2009 National Weather Service (NWS), +!/ Copyright 2009-2012 National Weather Service (NWS), !/ National Oceanic and Atmospheric Administration. All rights !/ reserved. WAVEWATCH III is a trademark of the NWS. !/ No unauthorized use without permission. @@ -872,7 +874,7 @@ !/ | WAVEWATCH III NOAA/NCEP | !/ | H. L. Tolman | !/ | FORTRAN 90 | -!/ | Last update : 11-May-2007 | +!/ | Last update : 13-Mar-2012 | !/ +-----------------------------------+ !/ !/ 19-Oct-1998 : Final FORTRAN 77 ( version 1.18 ) @@ -882,6 +884,8 @@ !/ 28-Jun-2005 : Adding MAPST2. ( version 3.07 ) !/ Taking out initilization. !/ 11-May-2007 : Adding NTPROC/NAPROC separation. ( version 3.11 ) +!/ 13-Mar-2012 : Add initialization of UST on re- ( version 3.14 ) +!/ activation of grid point. !/ ! 1. Purpose : ! @@ -938,7 +942,7 @@ !/ ------------------------------------------------------------------- / USE W3GDATMD, ONLY: NX, NY, NSEA, NSEAL, MAPSF, MAPSTA, MAPST2, & NTH, NK, NSPEC, SIG, TH, DTH, FICEN - USE W3WDATMD, ONLY: TIME, TICE, ICE + USE W3WDATMD, ONLY: TIME, TICE, ICE, UST !! USE W3ADATMD, ONLY: U10, U10D, CG USE W3ADATMD, ONLY: CG USE W3IDATMD, ONLY: TIN, ICEI @@ -1013,6 +1017,7 @@ IF ( ICEI(IX,IY).LT.FICEN .AND. MAPICE(IY,IX).EQ.1 ) THEN ! MAPICE(IY,IX) = 0 + UST(ISEA) = 0.05 ! IF ( MAPST2(IY,IX) .EQ. 0 ) THEN MAPSTA(IY,IX) = ABS(MAPSTA(IY,IX)) @@ -1069,7 +1074,7 @@ !/ | WAVEWATCH III NOAA/NCEP | !/ | H. L. Tolman | !/ | FORTRAN 90 | -!/ | Last update : 11-May-2007 | +!/ | Last update : 13-Mar-2012 | !/ +-----------------------------------+ !/ !/ 15-Jan-1998 : Final FORTRAN 77 ( version 1.18 ) @@ -1078,6 +1083,8 @@ !/ 15-Dec-2004 : Multiple grid version. ( version 3.06 ) !/ 15-Jul-2005 : Adding drying out of points. ( version 3.07 ) !/ 11-May-2007 : Adding NTPROC/NAPROC separation. ( version 3.11 ) +!/ 13-Mar-2012 : Add initialization of UST on re- ( version 3.14 ) +!/ activation of grid point. !/ ! 1. Purpose : ! @@ -1147,7 +1154,7 @@ USE W3GDATMD, ONLY: NX, NY, NSEA, NSEAL, MAPSF, MAPSTA, MAPST2, & ZB, DMIN, NK, NTH, NSPEC, SIG, DSIP, & MAPWN, MAPTH, FACHFA - USE W3WDATMD, ONLY: TIME, TLEV, WLV + USE W3WDATMD, ONLY: TIME, TLEV, WLV, UST USE W3ADATMD, ONLY: CG, WN, DW USE W3IDATMD, ONLY: TLN, WLEV USE W3SERVMD, ONLY: EXTCDE @@ -1282,6 +1289,7 @@ ! IF (WLV(ISEA)-ZB(ISEA).GT.0. .AND. MAPDRY(IY,IX).EQ.1) THEN MAPDRY(IY,IX) = 0 + UST(ISEA) = 0.05 IF ( MAPST2(IY,IX) .EQ. 0 ) THEN MAPSTA(IY,IX) = ABS(MAPSTA(IY,IX)) !/T2 WRITE (NDST,9021) ISEA, WLV(ISEA)-ZB(ISEA), &