diff --git a/model/src/w3pro2md.F90 b/model/src/w3pro2md.F90 index 4ab6a3e69..a23f893ef 100644 --- a/model/src/w3pro2md.F90 +++ b/model/src/w3pro2md.F90 @@ -967,22 +967,6 @@ SUBROUTINE W3XYP2 ( ISP, DTG, MAPSTA, MAPFS, VQ, VGX, VGY ) #endif DNN = XWIND * DCELL + (1.-XWIND) * DNND * TFAC - - CGD = 0.5 * GRAV / SIG(IK) * IOBDP_LOC(JSEA) - DSS = ( CGD * (XFR-1.))**2 * DTME / 12. - DNN = ( CGD * DTH )**2 * DTME / 12. - DCELL = CGD / 10.0 ! -> CELLP needs probably redifinition ... - KH = WN(IK,ISEA)*DW(ISEA) - XWIND = 3.3 * U10(ISEA)*WN(IK,ISEA)/SIG(IK) - 2.3 - XWIND = MAX ( 0. , MIN ( 1. , XWIND ) ) - TFAC = MIN ( 1. , (CLATS(ISEA)/CLATMN)**2 ) - DSS = SWFAC(JSEA) * (XWIND * DCELL + (1.-XWIND) * DSS * TFAC) - DNN = SWFAC(JSEA) * (XWIND * DCELL + (1.-XWIND) * DNN * TFAC) - DIFFVEC(1,JSEA) = (DSS*ECOS(ITH)**2+DNN*ESIN(ITH)**2) - DIFFVEC(2,JSEA) = (DSS*ESIN(ITH)**2+DNN*ECOS(ITH)**2) / CLATS(ISEA)**2 - DIFFVEC(3,JSEA) = ((DSS-DNN) * ESIN(ITH)*ECOS(ITH)) / CLATS(ISEA) - - VDXX(IXY) = DTLOC * (DSS*ECOS(ITH)**2+DNN*ESIN(ITH)**2) VDYY(IXY) = DTLOC * (DSS*ESIN(ITH)**2+DNN*ECOS(ITH)**2) & / CLATS(ISEA)**2