@@ -1492,79 +1492,79 @@ SUBROUTINE W3WAVE ( IMOD, ODAT, TEND, STAMP, NO_OUT &
1492
1492
1493
1493
IF ((IOBP_LOC(JSEA).eq. 1 ..or. IOBP_LOC(JSEA).eq. 3 ).and. IOBDP_LOC(JSEA).eq. 1.and .IOBPA_LOC(JSEA).eq. 0 ) THEN
1494
1494
1495
- IX = MAPSF(ISEA,1 )
1496
- IY = MAPSF(ISEA,2 )
1497
- DELA= 1 .
1498
- DELX= 1 .
1499
- DELY= 1 .
1495
+ IX = MAPSF(ISEA,1 )
1496
+ IY = MAPSF(ISEA,2 )
1497
+ DELA= 1 .
1498
+ DELX= 1 .
1499
+ DELY= 1 .
1500
1500
1501
1501
#ifdef W3_REF1
1502
- IF (GTYPE.EQ. RLGTYPE) THEN
1503
- DELX= SX* CLATS(ISEA)/ FACX
1504
- DELY= SY/ FACX
1505
- DELA= DELX* DELY
1506
- END IF
1507
- IF (GTYPE.EQ. CLGTYPE) THEN
1508
- ! Maybe what follows works also for RLGTYPE ... to be verified
1509
- DELX= HPFAC(IY,IX)/ FACX
1510
- DELY= HQFAC(IY,IX)/ FACX
1511
- DELA= DELX* DELY
1512
- END IF
1513
- REFLEC= REFLC(:,ISEA)
1514
- REFLEC(4 )= BERG(ISEA)* REFLEC(4 )
1515
- REFLED= REFLD(:,ISEA)
1502
+ IF (GTYPE.EQ. RLGTYPE) THEN
1503
+ DELX= SX* CLATS(ISEA)/ FACX
1504
+ DELY= SY/ FACX
1505
+ DELA= DELX* DELY
1506
+ END IF
1507
+ IF (GTYPE.EQ. CLGTYPE) THEN
1508
+ ! Maybe what follows works also for RLGTYPE ... to be verified
1509
+ DELX= HPFAC(IY,IX)/ FACX
1510
+ DELY= HQFAC(IY,IX)/ FACX
1511
+ DELA= DELX* DELY
1512
+ END IF
1513
+ REFLEC= REFLC(:,ISEA)
1514
+ REFLEC(4 )= BERG(ISEA)* REFLEC(4 )
1515
+ REFLED= REFLD(:,ISEA)
1516
1516
#endif
1517
1517
1518
1518
#ifdef W3_BT4
1519
- D50= SED_D50(ISEA)
1520
- PSIC= SED_PSIC(ISEA)
1519
+ D50= SED_D50(ISEA)
1520
+ PSIC= SED_PSIC(ISEA)
1521
1521
#endif
1522
- !
1522
+ !
1523
1523
#ifdef W3_DEBUGSRC
1524
- IF (IX .eq. DEBUG_NODE) THEN
1525
- WRITE (740 + IAPROC,* ) ' NODE_SRCE_IMP_PRE : IX=' , IX, ' JSEA=' , JSEA
1526
- END IF
1527
- WRITE (740 + IAPROC,* ) ' IT/IX/IY/IMOD=' , IT, IX, IY, IMOD
1528
- WRITE (740 + IAPROC,* ) ' ISEA/JSEA=' , ISEA, JSEA
1529
- WRITE (740 + IAPROC,* ) ' Before sum(VA)=' , sum (VA(:,JSEA))
1530
- FLUSH(740 + IAPROC)
1531
- #endif
1532
- CALL W3SRCE(srce_imp_pre, IT, ISEA, JSEA, IX, IY, IMOD, &
1533
- VAold(:,JSEA), VA(:,JSEA), &
1534
- VSioDummy, VDioDummy, SHAVETOT(JSEA), &
1535
- ALPHA(1 :NK,JSEA), WN(1 :NK,ISEA), &
1536
- CG(1 :NK,ISEA), CLATS(ISEA), DW(ISEA), U10(ISEA), &
1537
- U10D(ISEA), &
1524
+ IF (IX .eq. DEBUG_NODE) THEN
1525
+ WRITE (740 + IAPROC,* ) ' NODE_SRCE_IMP_PRE : IX=' , IX, ' JSEA=' , JSEA
1526
+ END IF
1527
+ WRITE (740 + IAPROC,* ) ' IT/IX/IY/IMOD=' , IT, IX, IY, IMOD
1528
+ WRITE (740 + IAPROC,* ) ' ISEA/JSEA=' , ISEA, JSEA
1529
+ WRITE (740 + IAPROC,* ) ' Before sum(VA)=' , sum (VA(:,JSEA))
1530
+ FLUSH(740 + IAPROC)
1531
+ #endif
1532
+ CALL W3SRCE(srce_imp_pre, IT, ISEA, JSEA, IX, IY, IMOD, &
1533
+ VAold(:,JSEA), VA(:,JSEA), &
1534
+ VSioDummy, VDioDummy, SHAVETOT(JSEA), &
1535
+ ALPHA(1 :NK,JSEA), WN(1 :NK,ISEA), &
1536
+ CG(1 :NK,ISEA), CLATS(ISEA), DW(ISEA), U10(ISEA), &
1537
+ U10D(ISEA), &
1538
1538
#ifdef W3_FLX5
1539
- TAUA(ISEA), TAUADIR(ISEA), &
1540
- #endif
1541
- AS(ISEA), UST(ISEA), &
1542
- USTDIR(ISEA), CX(ISEA), CY(ISEA), &
1543
- ICE(ISEA), ICEH(ISEA), ICEF(ISEA), &
1544
- ICEDMAX(ISEA), &
1545
- REFLEC, REFLED, DELX, DELY, DELA, &
1546
- TRNX(IY,IX), TRNY(IY,IX), BERG(ISEA), &
1547
- FPIS(ISEA), DTDYN(JSEA), &
1548
- FCUT(JSEA), DTGpre, TAUWX(JSEA), TAUWY(JSEA), &
1549
- TAUOX(JSEA), TAUOY(JSEA), TAUWIX(JSEA), &
1550
- TAUWIY(JSEA), TAUWNX(JSEA), &
1551
- TAUWNY(JSEA), PHIAW(JSEA), CHARN(JSEA), &
1552
- TWS(JSEA), PHIOC(JSEA), TMP1, D50, PSIC, TMP2, &
1553
- PHIBBL(JSEA), TMP3, TMP4, PHICE(JSEA), &
1554
- TAUOCX(JSEA), TAUOCY(JSEA), WNMEAN(JSEA), &
1555
- RHOAIR(ISEA), ASF(ISEA))
1556
- IF (.not. LSLOC) THEN
1557
- VSTOT(:,JSEA) = VSioDummy
1558
- VDTOT(:,JSEA) = VDioDummy
1559
- ENDIF
1539
+ TAUA(ISEA), TAUADIR(ISEA), &
1540
+ #endif
1541
+ AS(ISEA), UST(ISEA), &
1542
+ USTDIR(ISEA), CX(ISEA), CY(ISEA), &
1543
+ ICE(ISEA), ICEH(ISEA), ICEF(ISEA), &
1544
+ ICEDMAX(ISEA), &
1545
+ REFLEC, REFLED, DELX, DELY, DELA, &
1546
+ TRNX(IY,IX), TRNY(IY,IX), BERG(ISEA), &
1547
+ FPIS(ISEA), DTDYN(JSEA), &
1548
+ FCUT(JSEA), DTGpre, TAUWX(JSEA), TAUWY(JSEA), &
1549
+ TAUOX(JSEA), TAUOY(JSEA), TAUWIX(JSEA), &
1550
+ TAUWIY(JSEA), TAUWNX(JSEA), &
1551
+ TAUWNY(JSEA), PHIAW(JSEA), CHARN(JSEA), &
1552
+ TWS(JSEA), PHIOC(JSEA), TMP1, D50, PSIC, TMP2, &
1553
+ PHIBBL(JSEA), TMP3, TMP4, PHICE(JSEA), &
1554
+ TAUOCX(JSEA), TAUOCY(JSEA), WNMEAN(JSEA), &
1555
+ RHOAIR(ISEA), ASF(ISEA))
1556
+ IF (.not. LSLOC) THEN
1557
+ VSTOT(:,JSEA) = VSioDummy
1558
+ VDTOT(:,JSEA) = VDioDummy
1559
+ ENDIF
1560
1560
#ifdef W3_DEBUGSRC
1561
- WRITE (740 + IAPROC,* ) ' After sum(VA)=' , sum (VA(:,JSEA))
1562
- WRITE (740 + IAPROC,* ) ' sum(VSTOT)=' , sum (VSTOT(:,JSEA))
1563
- WRITE (740 + IAPROC,* ) ' sum(VDTOT)=' , sum (VDTOT(:,JSEA))
1564
- WRITE (740 + IAPROC,* ) ' SHAVETOT=' , SHAVETOT(JSEA)
1565
- FLUSH(740 + IAPROC)
1561
+ WRITE (740 + IAPROC,* ) ' After sum(VA)=' , sum (VA(:,JSEA))
1562
+ WRITE (740 + IAPROC,* ) ' sum(VSTOT)=' , sum (VSTOT(:,JSEA))
1563
+ WRITE (740 + IAPROC,* ) ' sum(VDTOT)=' , sum (VDTOT(:,JSEA))
1564
+ WRITE (740 + IAPROC,* ) ' SHAVETOT=' , SHAVETOT(JSEA)
1565
+ FLUSH(740 + IAPROC)
1566
1566
#endif
1567
- ENDIF
1567
+ ENDIF
1568
1568
END DO ! JSEA
1569
1569
END IF ! PDLIB
1570
1570
#endif
0 commit comments