
[Science] [BAS home] [Met home] [Beowulf home]	Antarctic Meteorology

yabeu ([cp s] Atmos diags to look at QT_POS ; 4.5)

Start time: 2000/12/1, length: 10/0/3
Processors: 1 * 2, Exec name: $RUNID.exe; Timestep (mins): 30 (atmos); 60 (ocean)
Start dumps: $DUMPS/xahiya.da00c10 $DUMPS/xahiyo.da00c10
Fortran mods
Diags
Userprestash files
See-also: yabkb
Notes:
1. Died at k3310 (ts 39456) with SIGFPE 8 (underflow?) which on re-run appeared to be in freedrift, which is only there for diag purposes anyway.
2. Problem *was* in icefreed but went away when compiled with -O 0 -g. So ignore it..
History: yabeu -> s -> f -> e -> c -> a -> yabad -> a -> yaauo -> n -> yaasc -> yaaoa

Purpose of Run (long)

Yes another EVP run. This time, I've included extra atmos diags (Q-on-levels, etc) to try to see if the "QT_POS" problem actually means anything. It needn't be EVP for QT_POS to show up, BTW.

Results

Differences against HadCM3 are small:

!p.multi=[0,3,1]
!y.margin=[10,5]
pp_plot,pp_xsect(a1),la=90,/blo                                          
pp_plot,pp_xsect(pp_diff(a1,ctl3)),la=90,/blo
pp_plot,pp_xsect(ctl3),la=90,/blo            
gettwogifs,out='spec.hum.sep'

b=get_fields('/data/beowulf2/wmc_32/archive/yabeu/yabeua.pm*feb.pp',seek='f.lbuser(3) eq 10',nf=19)
ctl3b=gm('ctl3/x.01m',m=2,sc=10,nl=19)                                                             
b1=pp_avg1(b,1) 
pp_plot,pp_xsect(b1),la=90,/blo                                                                    
pp_plot,pp_xsect(pp_diff(b1,ctl3b)),la=90,/blo
pp_plot,pp_xsect(ctl3b),la=90,/blo            
gettwogifs,out='spec.hum.feb'