[43] | 1 | ;+ |
---|
[44] | 2 | ; .. _fig1_no_obs_map.pro: |
---|
[43] | 3 | ; |
---|
| 4 | ; =================== |
---|
| 5 | ; fig1_no_obs_map.pro |
---|
| 6 | ; =================== |
---|
| 7 | ; |
---|
| 8 | ; DESCRIPTION |
---|
| 9 | ; =========== |
---|
| 10 | ; |
---|
[85] | 11 | ; :file:`${PROJECT_OD}/mask_oaflux_30N30S.nc` |
---|
| 12 | ; containing OAFLUX grid |
---|
[100] | 13 | ; has been produced by |
---|
[164] | 14 | ; :func:`oaflux_mask_30N30S`. |
---|
[85] | 15 | ; |
---|
[97] | 16 | ; Figure of |
---|
[85] | 17 | ; ++ |
---|
| 18 | ; is saved in |
---|
| 19 | ; :file:`${PROJECT_OD}/fig1_no_obs_map.ps`. |
---|
| 20 | ; |
---|
[44] | 21 | ; .. graphviz:: |
---|
| 22 | ; |
---|
| 23 | ; digraph fig1_no_obs_map { |
---|
| 24 | ; |
---|
[85] | 25 | ; mask [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/mask_oaflux_30N30S.nc"]; |
---|
[50] | 26 | ; lhf_no_obs [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lhf_no_obs.txt"]; |
---|
| 27 | ; sw_no_obs [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/sw_no_obs.txt"]; |
---|
| 28 | ; lw_no_obs [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lw_no_obs.txt"]; |
---|
| 29 | ; figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/fig1_no_obs_map.ps"]; |
---|
[44] | 30 | ; |
---|
[85] | 31 | ; fig1_no_obs_map [shape=box, |
---|
[44] | 32 | ; fontname=Courier, |
---|
| 33 | ; color=blue, |
---|
| 34 | ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/paper01/fig1/fig1_no_obs_map.pro", |
---|
| 35 | ; label="${TROPFLUX}/src/paper01/fig1/fig1_no_obs_map.pro"]; |
---|
| 36 | ; |
---|
| 37 | ; {mask lhf_no_obs sw_no_obs lw_no_obs} -> {fig1_no_obs_map} -> {figure} |
---|
[85] | 38 | ; |
---|
[44] | 39 | ; } |
---|
| 40 | ; |
---|
[43] | 41 | ; SEE ALSO |
---|
| 42 | ; ======== |
---|
| 43 | ; |
---|
[85] | 44 | ; :ref:`data_in_oaflux` |
---|
[44] | 45 | ; |
---|
[50] | 46 | ; :ref:`project_profile.sh` |
---|
| 47 | ; :ref:`project_init.pro` |
---|
[43] | 48 | ; |
---|
[46] | 49 | ; :func:`x_site_location` |
---|
| 50 | ; :func:`y_site_location` |
---|
| 51 | ; |
---|
[164] | 52 | ; :func:`oaflux_mask_30N30S` |
---|
[44] | 53 | ; |
---|
[43] | 54 | ; EXAMPLES |
---|
| 55 | ; ======== |
---|
| 56 | ; |
---|
| 57 | ; :: |
---|
| 58 | ; |
---|
[180] | 59 | ; fig1_no_obs_map |
---|
[43] | 60 | ; |
---|
[44] | 61 | ; TODO |
---|
| 62 | ; ==== |
---|
| 63 | ; |
---|
[85] | 64 | ; how text files have been produced ? |
---|
| 65 | ; |
---|
[44] | 66 | ; coding rules |
---|
| 67 | ; |
---|
| 68 | ; complete description |
---|
| 69 | ; |
---|
| 70 | ; handle IO error |
---|
| 71 | ; |
---|
[43] | 72 | ; EVOLUTIONS |
---|
| 73 | ; ========== |
---|
| 74 | ; |
---|
| 75 | ; $Id$ |
---|
| 76 | ; |
---|
[85] | 77 | ; - fplod 20110817T104358Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 78 | ; |
---|
| 79 | ; * get rid of spawn |
---|
| 80 | ; |
---|
[47] | 81 | ; - pinsard 2011-04-20T15:08:28Z loholt1.ipsl.polytechnique.fr (Linux) |
---|
| 82 | ; |
---|
| 83 | ; * make it work on loholt1 |
---|
| 84 | ; |
---|
[46] | 85 | ; - fplod 20110420T134821Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 86 | ; |
---|
| 87 | ; * externalize functions |
---|
| 88 | ; |
---|
[44] | 89 | ; - fplod 20110412T073818Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 90 | ; |
---|
| 91 | ; * make it work on cratos |
---|
| 92 | ; * remove hard coding path |
---|
| 93 | ; * add graphviz |
---|
| 94 | ; |
---|
[43] | 95 | ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 96 | ; |
---|
| 97 | ; * minimal header |
---|
| 98 | ; |
---|
| 99 | ;- |
---|
[41] | 100 | pro fig1_no_obs_map |
---|
[44] | 101 | ; |
---|
[85] | 102 | @cm_4cal |
---|
| 103 | @cm_4data |
---|
| 104 | @cm_4mesh |
---|
| 105 | @cm_4data |
---|
[50] | 106 | @cm_project |
---|
[44] | 107 | ; |
---|
[41] | 108 | reinitplt, /z,/invert |
---|
| 109 | key_portrait = 1 |
---|
[97] | 110 | ; |
---|
[85] | 111 | openps, FILENAME = project_od_env + 'fig1_no_obs_map.ps' |
---|
[41] | 112 | close,/all |
---|
[97] | 113 | ; |
---|
[41] | 114 | yr=365. |
---|
[97] | 115 | lh_mi=365/yr |
---|
| 116 | lh_ma=3285/yr |
---|
[94] | 117 | lh_int=0.25 |
---|
[97] | 118 | sw_mi=365/yr |
---|
| 119 | sw_ma=3285/yr |
---|
[94] | 120 | sw_int=0.25 |
---|
[97] | 121 | lw_mi=365/yr |
---|
| 122 | lw_ma=1460/yr |
---|
[94] | 123 | lw_int=0.125 |
---|
[41] | 124 | lct=65 |
---|
[97] | 125 | ; |
---|
[50] | 126 | file=project_id_env + 'mask_oaflux_30N30S.nc' |
---|
[41] | 127 | initncdf, file |
---|
| 128 | domdef, 30,390,-30,30 |
---|
| 129 | msk=ncdf_lec(file,var='msk') |
---|
[97] | 130 | ; |
---|
[50] | 131 | fi=project_id_env + 'lhf_no_obs.txt' |
---|
[41] | 132 | res=read_ascii(fi,data_start=1) |
---|
| 133 | ff=res.field1 |
---|
| 134 | lat=reform(ff(0,*)) |
---|
| 135 | lon=reform(ff(1,*)) |
---|
| 136 | no=reform(ff(2,*)) |
---|
| 137 | no=no/yr |
---|
[97] | 138 | ;ind=where(no ge lh_ma) |
---|
[94] | 139 | ; no(ind)=no-0.5 |
---|
[97] | 140 | ; |
---|
[41] | 141 | marge1=[0,0,-5,0] |
---|
| 142 | plt, .4+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,box=box,lct=lct, lh_mi, lh_ma, int=lh_int, $ |
---|
| 143 | title='a) Turbulent Flux validation sites', subtitle='', small=[1,4,1],/rempl, marge=marge1, format='(i5)' |
---|
| 144 | NN=n_elements(lat) |
---|
[97] | 145 | ; |
---|
[41] | 146 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 147 | ; |
---|
[41] | 148 | for n=0,NN-1 do begin |
---|
[44] | 149 | x=lon[n] |
---|
| 150 | y=lat[n] |
---|
| 151 | c=no[n] |
---|
[97] | 152 | cmi=lh_mi |
---|
| 153 | cma=lh_ma |
---|
[94] | 154 | dc=cma-cmi |
---|
[41] | 155 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 156 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 157 | endfor |
---|
[97] | 158 | ; |
---|
[41] | 159 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 160 | ; |
---|
[41] | 161 | for n=0,NN-1 do begin |
---|
[44] | 162 | x=lon[n] |
---|
| 163 | y=lat[n] |
---|
| 164 | c=no[n] |
---|
[97] | 165 | ; |
---|
[41] | 166 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 167 | endfor |
---|
[97] | 168 | ; |
---|
[50] | 169 | fi=project_id_env + 'sw_no_obs.txt' |
---|
[41] | 170 | res=read_ascii(fi,data_start=1) |
---|
| 171 | ff=res.field1 |
---|
| 172 | lat=reform(ff(0,*)) |
---|
| 173 | lon=reform(ff(1,*)) |
---|
| 174 | no=reform(ff(2,*)) |
---|
| 175 | no=no/yr |
---|
[97] | 176 | ;ind=where(no ge sw_ma) |
---|
[94] | 177 | ;no(ind)=sw_ma-0.5 |
---|
[41] | 178 | marge1=[0,0,-5,0] |
---|
| 179 | plt, .4+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,box=box,lct=lct, sw_mi, sw_ma, int=sw_int, $ |
---|
| 180 | title='b) SWR validation sites', subtitle='', small=[1,4,2],/noer,/rempl, marge=marge1, format='(i5)' |
---|
| 181 | NN=n_elements(lat) |
---|
[97] | 182 | ; |
---|
[41] | 183 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 184 | ; |
---|
[41] | 185 | for n=0,NN-1 do begin |
---|
[44] | 186 | x=lon[n] |
---|
| 187 | y=lat[n] |
---|
| 188 | c=no[n] |
---|
[97] | 189 | cmi=sw_mi |
---|
| 190 | cma=sw_ma |
---|
[94] | 191 | dc=cma-cmi |
---|
[41] | 192 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 193 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 194 | endfor |
---|
[97] | 195 | ; |
---|
[41] | 196 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 197 | ; |
---|
[41] | 198 | for n=0,NN-1 do begin |
---|
[44] | 199 | x=lon[n] |
---|
| 200 | y=lat[n] |
---|
| 201 | c=no[n] |
---|
[97] | 202 | ; |
---|
[41] | 203 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 204 | endfor |
---|
[97] | 205 | ; |
---|
[50] | 206 | fi=project_id_env + 'lw_no_obs.txt' |
---|
[41] | 207 | res=read_ascii(fi,data_start=1) |
---|
| 208 | ff=res.field1 |
---|
| 209 | lat=reform(ff(0,*)) |
---|
| 210 | lon=reform(ff(1,*)) |
---|
| 211 | no=reform(ff(2,*)) |
---|
| 212 | no=no/yr |
---|
[97] | 213 | ;ind=where(no ge lw_ma) |
---|
[94] | 214 | ;no(ind)=lw_ma-0.25 |
---|
[97] | 215 | ; |
---|
[41] | 216 | marge1=[0,0,-5,0] |
---|
| 217 | plt, .4+msk*1.e-5,realcont=2,/nocont,/nofill,/noer,xminor=1,yminor=1,box=box,lct=lct,lw_mi, lw_ma, int=lw_int, $ |
---|
| 218 | title='c) LWR validation sites', subtitle='', small=[1,4,3],/rempl, marge=marge1, format='(i5)' |
---|
| 219 | NN=n_elements(lat) |
---|
[97] | 220 | ; |
---|
[41] | 221 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 222 | ; |
---|
[41] | 223 | for n=0,NN-1 do begin |
---|
[44] | 224 | x=lon[n] |
---|
| 225 | y=lat[n] |
---|
| 226 | c=no[n] |
---|
[97] | 227 | cmi=lw_mi |
---|
| 228 | cma=lw_ma |
---|
[94] | 229 | dc=cma-cmi |
---|
[41] | 230 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 231 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 232 | endfor |
---|
[97] | 233 | ; |
---|
[41] | 234 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 235 | ; |
---|
[41] | 236 | for n=0,NN-1 do begin |
---|
[44] | 237 | x=lon[n] |
---|
| 238 | y=lat[n] |
---|
| 239 | c=no[n] |
---|
[97] | 240 | ; |
---|
[41] | 241 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 242 | endfor |
---|
[97] | 243 | ; |
---|
| 244 | ; independent validation sites |
---|
| 245 | ; |
---|
| 246 | plots, 61.5, 15.5,psym=2 ; as location |
---|
| 247 | plots, 331, 25.5,psym=2 ; subduction (central) |
---|
| 248 | plots, 326, 18,psym=2 ; subduction (sw) |
---|
| 249 | plots, 338, 18,psym=2 ; subduction (se) |
---|
| 250 | plots, 156, -1.7 , psym=2 ; coare |
---|
| 251 | plots, 274.2, -21, psym=2 ; stratus |
---|
| 252 | plots, 202.46, 22.45, psym=2 ; WHOTS |
---|
| 253 | ; |
---|
[41] | 254 | closeps |
---|
[97] | 255 | ; |
---|
[41] | 256 | end |
---|