Skip to content

Commit

Permalink
ww3_gse2: some leftovers ...
Browse files Browse the repository at this point in the history
  • Loading branch information
aronroland committed Dec 19, 2024
1 parent 98690ce commit 667c2dc
Showing 1 changed file with 0 additions and 16 deletions.
16 changes: 0 additions & 16 deletions model/src/w3pro2md.F90
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down

0 comments on commit 667c2dc

Please sign in to comment.