I know this has been posted before. While I still have a little confusion about the humidity used in bulk flux script.
IF (RH.lt.2.0_r8) THEN !RH fraction
cff=cff*RH !Vapor pres (mb)
Q(i)=0.62197_r8*(cff/(PairM-0.378_r8*cff)) !Spec hum (kg/kg)
ELSE !RH input was actually specific humidity in g/kg
Q(i)=RH/1000.0_r8 !Spec Hum (kg/kg)
END IF
This is the humidity code I found in ROMS. It looks like ROMS automatically reads in the humidity inputs as relative humidity if the value is smaller than 2 ,while read in data as specific humidity if the value is larger than 2 and divided it by 1000(duo to the different units). In my case, I only have specific humidity of around (0.013kg/kg), so I multiples it by 1000 to 13g/kg. Then the model will switch it to specific humidity by the code above and turns it to 0.013kg/kg. I suppose my understanding is correct here.
My question is does it matter if the variable name is relative humidity. As my log file still shows this:
GET_2DFLD - surface air relative humidity, t = 40117 12:00:00
(Rec=0000007, Index=1, File: Pilbara_grid_heat.nc)
(Tmin= 40116.0000 Tmax= 40176.0000)
(Min = 1.48332236E-01 Max = 1.70468596E-01)
The value is definitely not relative humidity but specific humidity(0.015~0.017). I assume ROMS identify those two humidity by value, not by name. SO even it shows relative humidity here, it actually means specific humidity.Is that right?
relative humidity
-
- Posts: 60
- Joined: Tue Sep 29, 2009 3:50 pm
- Location: School of Environment System Engineering,UWA