| 1 |
|
| 2 |
#include "ctrparam.h" |
| 3 |
|
| 4 |
SUBROUTINE SDRAG(WLMMAX,JWLMMAX) 7801. |
| 5 |
C**** 7802. |
| 6 |
C**** THIS SUBROUTINE PUTS A DRAG ON THE WINDS ON THE TOP LAYER OF 7803. |
| 7 |
C**** THE ATMOSPHERE 7804. |
| 8 |
C**** 7805. |
| 9 |
#include "BD2G04.COM" 7806. |
| 10 |
|
| 11 |
COMMON U,V,T,P,Q 7807. |
| 12 |
COMMON/WORK1/CONV(IM0,JM0,LM0),PK(IM0,JM0,LM0) 7808. |
| 13 |
real TAUJ(JM0+1) |
| 14 |
DO 100 J=2,JM 7809. |
| 15 |
I=IM 7810. |
| 16 |
DO 100 IP1=1,IM 7811. |
| 17 |
SP=.25*(P(I,J-1)+P(IP1,J-1)+P(I,J)+P(IP1,J)) 7812. |
| 18 |
WLM=SQRT(U(I,J,LM)*U(I,J,LM)+V(I,J,LM)*V(I,J,LM)) 7813. |
| 19 |
if(WLM.gt.WLMMAX)then |
| 20 |
WLMMAX=WLM |
| 21 |
JWLMMAX=J |
| 22 |
endif |
| 23 |
RHO=PTOP/(RGAS*T(I,J,LM)*PK(I,J,LM)) 7814. |
| 24 |
CDN=DUMMY1(1)+DUMMY1(2)*WLM 7815. |
| 25 |
TAUJ(J)=CDN*100.*RHO*WLM*U(I,J,LM) |
| 26 |
X=NDYN*DT*RHO*CDN*WLM*GRAV/(SP*DSIG(LM)) 7816. |
| 27 |
U(I,J,LM)=U(I,J,LM)*(1.-X) 7817. |
| 28 |
V(I,J,LM)=V(I,J,LM)*(1.-X) 7818. |
| 29 |
100 I=IP1 7819. |
| 30 |
C**** 5973. |
| 31 |
TAUJ(1)=TAUJ(2) |
| 32 |
TAUJ(JM+1)=TAUJ(JM) |
| 33 |
do J=1,JM |
| 34 |
PLAND=FDATA(I,J,2) 5974. |
| 35 |
PWATER=1.-PLAND 5975. |
| 36 |
PLICE=FDATA(I,J,3)*PLAND 5976. |
| 37 |
PEARTH=PLAND-PLICE 5977. |
| 38 |
POICE=ODATA(I,J,2)*PWATER 5978. |
| 39 |
POCEAN=PWATER-POICE 5979. |
| 40 |
if(POCEAN.LE.1.E-5)then |
| 41 |
POCEAN=0. |
| 42 |
POICE=PWATER |
| 43 |
endif |
| 44 |
TTOFR=PEARTH+PLICE+POICE+POCEAN |
| 45 |
if(abs(TTOFR-1.).gt.1.e-3)then |
| 46 |
print *,'From sdrag TTOFR=',TTOFR |
| 47 |
print *,' J=',J,' PLAND=',PLAND,' POCEAN=',POCEAN |
| 48 |
print *,' I=',I,' PWATER=',PWATER,' POICE=',POICE |
| 49 |
print *, 'ODATA(I,J,2)=',ODATA(I,J,2) |
| 50 |
stop |
| 51 |
end if |
| 52 |
C Fric |
| 53 |
AJ(J,38)=AJ(J,38)+(TAUJ(j)+TAUJ(j+1))*POCEAN |
| 54 |
BJ(J,38)=BJ(J,38)+(TAUJ(j)+TAUJ(j+1))*PLAND |
| 55 |
CJ(J,38)=CJ(J,38)+(TAUJ(j)+TAUJ(j+1))*POICE |
| 56 |
enddo |
| 57 |
C Fric |
| 58 |
RETURN 7820. |
| 59 |
END 7821. |
Parent Directory
| 
Revision Log
| 
Revision Graph