Commit 0b7fecad authored by kbk's avatar kbk
Browse files

added latitude dependent cloud correction factor for long wave rad. - Stips

parent d572b7b2
!$Id: fluxes.F90,v 1.11 2005-04-19 15:21:33 kbk Exp $
!$Id: fluxes.F90,v 1.12 2005-04-19 15:56:58 kbk Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -6,7 +6,7 @@
! !ROUTINE: Heat and momentum fluxes.
!
! !INTERFACE:
subroutine fluxes(u10,v10,airt,tcc,sst,hf,taux,tauy)
subroutine fluxes(lat,u10,v10,airt,tcc,sst,hf,taux,tauy)
!
! !DESCRIPTION:
! The sum of the latent and sensible heat fluxes + longwave
......@@ -35,7 +35,7 @@
IMPLICIT NONE
!
! !INPUT PARAMETERS:
REALTYPE, intent(in) :: u10,v10,airt,tcc,sst
REALTYPE, intent(in) :: lat,u10,v10,airt,tcc,sst
!
! !INPUT/OUTPUT PARAMETERS:
!
......@@ -46,7 +46,10 @@
! Original author(s): Karsten Bolding and Hans Burchard
!
! $Log: fluxes.F90,v $
! Revision 1.11 2005-04-19 15:21:33 kbk
! Revision 1.12 2005-04-19 15:56:58 kbk
! added latitude dependent cloud correction factor for long wave rad. - Stips
!
! Revision 1.11 2005/04/19 15:21:33 kbk
! cleaned long wave radiation and added Bignami and Berliand (ROMS)
!
! Revision 1.10 2005/04/19 13:02:08 kbk
......@@ -89,6 +92,29 @@
integer, parameter :: bignami=3 ! Bignami 1995 -Medsea
integer, parameter :: berliand=4 ! Berliand 1952 -ROMS
!
#ifndef OLD_WRONG_FLUXES
real, parameter, dimension(91) :: cloud_correction_factor = (/ &
0.497202, 0.501885, 0.506568, 0.511250, 0.515933, &
0.520616, 0.525299, 0.529982, 0.534665, 0.539348, &
0.544031, 0.548714, 0.553397, 0.558080, 0.562763, &
0.567446, 0.572129, 0.576812, 0.581495, 0.586178, &
0.590861, 0.595544, 0.600227, 0.604910, 0.609593, &
0.614276, 0.618959, 0.623641, 0.628324, 0.633007, &
0.637690, 0.642373, 0.647056, 0.651739, 0.656422, &
0.661105, 0.665788, 0.670471, 0.675154, 0.679837, &
0.684520, 0.689203, 0.693886, 0.698569, 0.703252, &
0.707935, 0.712618, 0.717301, 0.721984, 0.726667, &
0.731350, 0.736032, 0.740715, 0.745398, 0.750081, &
0.754764, 0.759447, 0.764130, 0.768813, 0.773496, &
0.778179, 0.782862, 0.787545, 0.792228, 0.796911, &
0.801594, 0.806277, 0.810960, 0.815643, 0.820326, &
0.825009, 0.829692, 0.834375, 0.839058, 0.843741, &
0.848423, 0.853106, 0.857789, 0.862472, 0.867155, &
0.871838, 0.876521, 0.881204, 0.885887, 0.890570, &
0.895253, 0.899936, 0.904619, 0.909302, 0.913985, &
0.918668 /)
#endif
!
! !LOCAL VARIABLES:
#ifdef EA_ZERO
REALTYPE :: ea=_ZERO_
......@@ -121,6 +147,11 @@
qh=cd_heat*cpa*rho_air*w*(tw-ta) ! sensible
qe=cd_latent*L*rho_air*w*(qs-qa) ! latent
#ifndef OLD_WRONG_FLUXES
! calculate cloud correction factor North and South are equal
ccf = cloud_correction_factor(nint(abs(lat))+1)
#endif
select case(back_radiation_method) ! back radiation
case(clark)
! AS unit of ea is Pascal, must hPa
......
!$Id: meteo.F90,v 1.11 2005-01-13 09:49:37 kbk Exp $
!$Id: meteo.F90,v 1.12 2005-04-19 15:56:58 kbk Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -78,7 +78,10 @@
! Original author(s): Karsten Bolding & Hans Burchard
!
! $Log: meteo.F90,v $
! Revision 1.11 2005-01-13 09:49:37 kbk
! Revision 1.12 2005-04-19 15:56:58 kbk
! added latitude dependent cloud correction factor for long wave rad. - Stips
!
! Revision 1.11 2005/01/13 09:49:37 kbk
! wet bulb works, es is global, cleaning - Stips
!
! Revision 1.10 2004/01/15 11:45:00 kbk
......@@ -439,8 +442,9 @@
call exchange_coefficients( &
u10(i,j),v10(i,j),t2(i,j),airp(i,j), &
sst(i,j),hum(i,j),hum_method)
call fluxes(u10(i,j),v10(i,j),t2(i,j),tcc(i,j), &
sst(i,j),shf(i,j),tausx(i,j),tausy(i,j))
call fluxes(latc(i,j),u10(i,j),v10(i,j), &
t2(i,j),tcc(i,j),sst(i,j), &
shf(i,j),tausx(i,j),tausy(i,j))
else
shf(i,j) = _ZERO_
tausx(i,j) = _ZERO_
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment