| 1 | subroutine PHY________OUT(MOTIVATION) |
|---|
| 2 | |
|---|
| 3 | !------------------------------------------------------------------------------+ |
|---|
| 4 | ! Sat 22-Jun-2013 MAR | |
|---|
| 5 | ! MAR PHY________OUT | |
|---|
| 6 | ! subroutine PHY________OUT OUTPUTs MAR PHYsical parameterizations | |
|---|
| 7 | ! | |
|---|
| 8 | ! version 3.p.4.1 created by H. Gallee, Tue 12-Mar-2013 | |
|---|
| 9 | ! Last Modification by H. Gallee, Sat 22-Jun-2013 | |
|---|
| 10 | ! | |
|---|
| 11 | !------------------------------------------------------------------------------+ |
|---|
| 12 | |
|---|
| 13 | use Mod_Real |
|---|
| 14 | use Mod_PHY____dat |
|---|
| 15 | use Mod_PHY____grd |
|---|
| 16 | use Mod_PHY_DY_kkl |
|---|
| 17 | use Mod_PHY_CM_kkl |
|---|
| 18 | use Mod_PHY_AT_kkl |
|---|
| 19 | use Mod_SISVAT_gpt |
|---|
| 20 | |
|---|
| 21 | |
|---|
| 22 | IMPLICIT NONE |
|---|
| 23 | |
|---|
| 24 | |
|---|
| 25 | character(len=50) :: MOTIVATION |
|---|
| 26 | integer :: i, j, ikl ! |
|---|
| 27 | integer :: k, kz ! |
|---|
| 28 | real(kind=real8) :: RelHum |
|---|
| 29 | |
|---|
| 30 | |
|---|
| 31 | |
|---|
| 32 | |
|---|
| 33 | ! OUTPUT |
|---|
| 34 | ! ====== |
|---|
| 35 | |
|---|
| 36 | |
|---|
| 37 | ! txt file |
|---|
| 38 | ! -------- |
|---|
| 39 | |
|---|
| 40 | i = i_x0 |
|---|
| 41 | j = j_y0 |
|---|
| 42 | ikl = ikl0 |
|---|
| 43 | |
|---|
| 44 | write(4,50) MOTIVATION |
|---|
| 45 | 50 format(3(/,1x),' OUTPUT for VERIFICATION: ',a50, & |
|---|
| 46 | & /,1x ,' ***********************',/,1x) |
|---|
| 47 | |
|---|
| 48 | write(4,55) Day_TU,Mon_TU,YearTU,HourTU,MinuTU,Sec_TU,it_EXP |
|---|
| 49 | 55 format(3x,2(i2,'-'),i4,4x,3(i2,'-'),' Simulation Iteration No ',i6) |
|---|
| 50 | |
|---|
| 51 | write(4,62) |
|---|
| 52 | write(4,60) |
|---|
| 53 | 60 format(7x,'|',' sigma ',' |',' Z ',' |',' T ',' |' &! |
|---|
| 54 | & ,' U ',' |',' V ',' |',' RH ',' |',' Qv ',' |' &! |
|---|
| 55 | & ,' Qw ',' |',' Qi ',' |' &! |
|---|
| 56 | & ,' TKE ',' |',' eps ',' |' &! |
|---|
| 57 | & ,' Kzh L',' |',' Kzh ',' |' &! |
|---|
| 58 | &,/,' |',' ',' |',' [km] ',' |',' [K] ',' |' &! |
|---|
| 59 | & ,' [m/s]',' |',' [m/s]',' |',' [-] ',' |',' [g/kg]',' |' &! |
|---|
| 60 | & ,' [g/kg]',' |',' [g/kg]',' |' &! |
|---|
| 61 | & ,'[m2/s2]',' |','[m3/s2]',' |' &! |
|---|
| 62 | & ,' [m2/s]',' |',' [m2/s]',' |' &! |
|---|
| 63 | &,/,' ------+','---------','-+','---------','-+','---------','-+' &! |
|---|
| 64 | & ,'------','-+','------','-+','------','-+','-------','-+' &! |
|---|
| 65 | & ,'-------','-+','-------','-+' &! |
|---|
| 66 | & ,'-------','-+','-------','-+' &! |
|---|
| 67 | & ,'-------','-+','-------','-+') ! |
|---|
| 68 | DO k = 1,mzp |
|---|
| 69 | kz=min(k,mzp) |
|---|
| 70 | RelHum = min(un_1,qv__DY(ikl,k) /max(qvswCM(ikl,k),eps6)) |
|---|
| 71 | write(4,61) k,sigma(k),Z___DY(ikl,k)*1.e-3 ,Ta__DY(ikl,k) &! |
|---|
| 72 | & ,Ua__DY(ikl,k) ,Va__DY(ikl,k) &! |
|---|
| 73 | & ,RelHum , 1.e3*qv__DY(ikl,k) &! |
|---|
| 74 | & ,1.e3*qw__CM(ikl,kz), 1.e3*qi__CM(ikl,kz) &! |
|---|
| 75 | & , TKE_AT(ikl,kz), eps_AT(ikl,kz) &! |
|---|
| 76 | & , Kzh_AT(ikl,kz), Kzh0AT(ikl,kz) |
|---|
| 77 | 61 format(i6,' |',f9.6,' |', f9.4,' |' , f9.3,' |' ,2(f6.1,' |') &! |
|---|
| 78 | & ,f6.2,' |', f7.3,' |' ,2(f7.3,' |') &! |
|---|
| 79 | & ,2(f7.2,' |'),2(f7.3,' |') ) |
|---|
| 80 | IF (mod(k,15).EQ.0) THEN |
|---|
| 81 | write(4,62) |
|---|
| 82 | 62 format( &! |
|---|
| 83 | & ' ------+','---------','-+','---------','-+','---------','-+' &! |
|---|
| 84 | & ,'------','-+','------','-+','------','-+','-------','-+' &! |
|---|
| 85 | & ,'-------','-+','-------','-+' &! |
|---|
| 86 | & ,'-------','-+','-------','-+' &! |
|---|
| 87 | & ,'-------','-+','-------','-+') ! |
|---|
| 88 | |
|---|
| 89 | write(4,60) |
|---|
| 90 | END IF |
|---|
| 91 | ENDDO |
|---|
| 92 | |
|---|
| 93 | write(4,62) |
|---|
| 94 | write(4,63) uts_SV_gpt(ikl),1.e3*uqs_SV_gpt(ikl),us__SV_gpt(ikl) & |
|---|
| 95 | & , HsenSV_gpt(ikl), HLatSV_gpt(ikl) |
|---|
| 96 | 63 format( &! |
|---|
| 97 | & 6x,' |', 9x ,' |', ( 9x ,' |'),' u*T* |' &! |
|---|
| 98 | & ,2( 6x ,' |'), ' u*q* |' &! |
|---|
| 99 | & , 7x ,' |' ,2( 7x ,' |'),' u* |' &! |
|---|
| 100 | & , 7x ,' |' ,2( 7x ,' |') &! |
|---|
| 101 | & ,/,6x,' |', 9x ,' |', ( 9x ,' |'), f9.6,' |' &! |
|---|
| 102 | & ,2( 6x ,' |'),f6.3,' |' &! |
|---|
| 103 | & , 7x ,' |' ,2( 7x ,' |'),f7.2,' |' &! |
|---|
| 104 | & , 7x ,' |' ,2( 7x ,' |') &! |
|---|
| 105 | & ,/,6x,' |', 9x ,' |', ( 9x ,' |'), f9.3,' |' &! |
|---|
| 106 | & ,2( 6x ,' |'),f6.0,' |' &! |
|---|
| 107 | & , 7x ,' |' ,2( 7x ,' |'), 7x ,' |' &! |
|---|
| 108 | & , 7x ,' |' ,2( 7x ,' |') &! |
|---|
| 109 | & ,/,6x,' |', 9x ,' |', ( 9x ,' |'),' [W/m2] |' &! |
|---|
| 110 | & ,2( 6x ,' |'), ' [W/m2]|' &! |
|---|
| 111 | & , 7x ,' |' ,2( 7x ,' |'),' |' &! |
|---|
| 112 | & , 7x ,' |' ,2( 7x ,' |') ) |
|---|
| 113 | |
|---|
| 114 | k = mzpp |
|---|
| 115 | write(4,62) |
|---|
| 116 | RelHum = min(un_1,qv__DY(ikl,k) /max(qvswCM(ikl,k),eps6)) |
|---|
| 117 | write(4,64) k,sigma(k),Z___DY(ikl,k)*1.e-3 ,Ta__DY(ikl,k) &! |
|---|
| 118 | & ,zer0 ,zer0 &! |
|---|
| 119 | & ,RelHum , 1.e3*qv__DY(ikl,k) |
|---|
| 120 | 64 format(i6,' |',f9.6,' |', f9.4,' |' , f9.3,' |' ,2(f6.1,' |') &! |
|---|
| 121 | & ,f6.2,' |', f7.3,' |' ,2( 7x ,' |') &! |
|---|
| 122 | & ,2( 7x ,' |'),2( 7x ,' |') ) |
|---|
| 123 | write(4,62) |
|---|
| 124 | |
|---|
| 125 | |
|---|
| 126 | |
|---|
| 127 | ! cdf file |
|---|
| 128 | ! -------- |
|---|
| 129 | |
|---|
| 130 | |
|---|
| 131 | end subroutine PHY________OUT |
|---|
