Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Thu Jul 18, 2019 8:00 pm




Post new topic Reply to topic  [ 6 posts ] 

All times are UTC

Author Message
 Post subject: about QCORRECTION
PostPosted: Wed Jul 08, 2009 3:59 pm 
Offline
User avatar

Joined: Wed Aug 02, 2006 1:21 pm
Posts: 23
Location: TAMU
Hello,

Is dqdt in set_vbc.F same as dQdSST in forcing file?

Quote:
stflx(i,j,itemp)=stflx(i,j,itemp)+ &
& dqdt(i,j)*(t(i,j,N(ng),nrhs,itemp)-sst(i,j))


If so, the unit of stflx is "degC m/s", while the unit of dQdSST is "Watts meter-2 Celsius-1".
It seems the units of both sides of the equation are not matching.

If not, what's the relationship between them?

Thanks for any reply.


Top
 Profile  
Reply with quote  
 Post subject: Re: about QCORRECTION
PostPosted: Thu Jul 09, 2009 12:58 am 
Offline

Joined: Fri Oct 12, 2007 5:59 pm
Posts: 8
Location: University of California, Irvine
I think they are the same. When you set the value of dqdt in "ana_sst.h", you need to divide it by (rho0*Cp), and then the unit of dqdt is m/s.


Top
 Profile  
Reply with quote  
 Post subject: Re: about QCORRECTION
PostPosted: Thu Jul 09, 2009 1:31 pm 
Offline

Joined: Mon Nov 10, 2008 8:44 pm
Posts: 4
Location: IOCAS, TAMU
In my ana_sst.h
Quote:
!-----------------------------------------------------------------------
! Set sea surface temperature (Celsius) and heat flux sensitivity to
! SST (Watts/m2).
!-----------------------------------------------------------------------
!
DO j=JstrR,JendR
DO i=IstrR,IendR
sst(i,j)=???
dqdt(i,j)=???
END DO
END DO
ana_sst.h: no values provided for sst and dqdt.


and
Quote:
! expressions. The forcing dQdSTT is usually computed in units of !
! (Watts/m2/degC). It needs to be scaled to (m/s) by dividing by !
! rho0*Cp.


so easy to miss it..., what's the consequence if just keep it the way as dqdt(i,j)=???.... during the compiling there is no warning


Top
 Profile  
Reply with quote  
 Post subject: Re: about QCORRECTION
PostPosted: Thu Jul 09, 2009 6:13 pm 
Offline
User avatar

Joined: Wed Aug 02, 2006 1:21 pm
Posts: 23
Location: TAMU
meiw wrote:
I think they are the same. When you set the value of dqdt in "ana_sst.h", you need to divide it by (rho0*Cp), and then the unit of dqdt is m/s.


Do I need to divide dQdSST by (rho0*Cp) when I am preparing the forcing file?

In the varinfo.dat:
Quote:
'dQdSST' ! Input
'surface net heat flux sensitivity to SST'
'watt meter-2 Celsius-1' ! Input: [Watt/m2]
'dQdSST, scalar, series' ! [Celsius m/s]
'sst_time'
'iddQdT'
'r2dvar'
1.0d0


I think it suggests that I don't need to divide it by (rho0*Cp) when I prepare dQdSST


Top
 Profile  
Reply with quote  
 Post subject: Re: about QCORRECTION
PostPosted: Thu Jul 09, 2009 7:08 pm 
Offline

Joined: Fri Oct 12, 2007 5:59 pm
Posts: 8
Location: University of California, Irvine
chimera.u wrote:
In my ana_sst.h
Quote:
!-----------------------------------------------------------------------
! Set sea surface temperature (Celsius) and heat flux sensitivity to
! SST (Watts/m2).
!-----------------------------------------------------------------------
!
DO j=JstrR,JendR
DO i=IstrR,IendR
sst(i,j)=???
dqdt(i,j)=???
END DO
END DO
ana_sst.h: no values provided for sst and dqdt.


and
Quote:
! expressions. The forcing dQdSTT is usually computed in units of !
! (Watts/m2/degC). It needs to be scaled to (m/s) by dividing by !
! rho0*Cp.


so easy to miss it..., what's the consequence if just keep it the way as dqdt(i,j)=???.... during the compiling there is no warning


I think this depends on whether "QCORRECTION" is defined in your header file or not.


Top
 Profile  
Reply with quote  
 Post subject: Re: about QCORRECTION
PostPosted: Thu Jul 09, 2009 7:16 pm 
Offline

Joined: Fri Oct 12, 2007 5:59 pm
Posts: 8
Location: University of California, Irvine
lmkli wrote:
meiw wrote:
I think they are the same. When you set the value of dqdt in "ana_sst.h", you need to divide it by (rho0*Cp), and then the unit of dqdt is m/s.


Do I need to divide dQdSST by (rho0*Cp) when I am preparing the forcing file?

In the varinfo.dat:
Quote:
'dQdSST' ! Input
'surface net heat flux sensitivity to SST'
'watt meter-2 Celsius-1' ! Input: [Watt/m2]
'dQdSST, scalar, series' ! [Celsius m/s]
'sst_time'
'iddQdT'
'r2dvar'
1.0d0


I think it suggests that I don't need to divide it by (rho0*Cp) when I prepare dQdSST


I think you still need to divide it by (rho0*Cp) in your forcing file. Actually, it is easy to check whether this is right or not since (rho0*Cp) is a very big number. You can run and compare two experiments: one with dqdt divided by (rho0*Cp) and the other without.


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 6 posts ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 2 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
cron
Powered by phpBB® Forum Software © phpBB Group