/[MITgcm]/MITgcm_contrib/jmc_script/grph_MON.m
ViewVC logotype

Diff of /MITgcm_contrib/jmc_script/grph_MON.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by jmc, Mon Nov 9 17:36:04 2020 UTC revision 1.4 by jmc, Wed Oct 19 15:29:58 2022 UTC
# Line 18  namLg=namA ; namLg=strrep(namLg,'_','\_' Line 18  namLg=namA ; namLg=strrep(namLg,'_','\_'
18  %-----------  %-----------
19  %- ngEn = Nb of Energy plot: = 2 (ke: Mx+Av) ; = 3 (+AM) or = 4 (add Pe ?) ;  %- ngEn = Nb of Energy plot: = 2 (ke: Mx+Av) ; = 3 (+AM) or = 4 (add Pe ?) ;
20  ngEn=2;  ngEn=2;
21    %- in case PE is missing from Monitor, use Mean-Eta and this "gdH" ratio to get PE
22  gdH=9.81/1000 ; % ratio gravity / mean_H (ocean) ; mean_Bo / mean_Pground (atmos)  gdH=9.81/1000 ; % ratio gravity / mean_H (ocean) ; mean_Bo / mean_Pground (atmos)
23  %- test if the variable krd is define :  %- test if the variable krd is define :
24  if size(who('krd'),1) > 0,  if size(who('krd'),1) > 0,
# Line 26  else Line 27  else
27   fprintf('krd undefined ; set to 1 \n'); krd=1 ;   fprintf('krd undefined ; set to 1 \n'); krd=1 ;
28  end  end
29  if krd == 1,  if krd == 1,
30    % list_on : controls which field to read-in (and plot):
31    %  1 : KE ; 2 : Eta ; 3,4 : T,S ; 5,6 : U,V ; 7 : W ; 8 : CFL ; 9 : Vort ; 10 : "sc"
32  %list_on=[1 1 1 1 0 0 1 1 0 0] ;  %list_on=[1 1 1 1 0 0 1 1 0 0] ;
33   list_on=[1 1 1 1 1 1 1 1 0 0] ;   list_on=[1 1 1 1 1 1 1 1 0 0] ;
34  %list_on=[1 1 0 0 0 0 0 0 1 1] ;  %list_on=[1 1 0 0 0 0 0 0 1 1] ;
# Line 65  if krd == 1, Line 68  if krd == 1,
68    for n=1:Nexp,    for n=1:Nexp,
69      keA(:,1,n)=etA(:,2,n)-etA(:,1,n);      keA(:,1,n)=etA(:,2,n)-etA(:,1,n);
70      ddKe=max(keA(:,4,n));      ddKe=max(keA(:,4,n));
71       if ddKe > 0, keA(:,3,n)=keA(:,4,n);       if ddKe > 0, keA(:,3,n)=keA(:,4,n); %- take PE from Monitor (since it's there)
72       else keA(:,3,n)=etA(:,4,n).*etA(:,4,n);       else keA(:,3,n)=etA(:,4,n).*etA(:,4,n);
73            keA(:,3,n)=keA(:,3,n)*gdH/2; end            keA(:,3,n)=keA(:,3,n)*gdH/2; end
74      keA(:,4,n)=keA(:,3,n)+keA(:,2,n);      keA(:,4,n)=keA(:,3,n)+keA(:,2,n);
# Line 74  if krd == 1, Line 77  if krd == 1,
77     tt4e=sprintf([tt4e,' %1.1e ;'],keA(1,4,n)); keA(:,4,n)=keA(:,4,n)-keA(1,4,n);     tt4e=sprintf([tt4e,' %1.1e ;'],keA(1,4,n)); keA(:,4,n)=keA(:,4,n)-keA(1,4,n);
78    end    end
79   else tt2e=' ' ; tt3e=' ' ; tt4e=' ' ; end;   else tt2e=' ' ; tt3e=' ' ; tt4e=' ' ; end;
80   ttA=ttA/3600; titT='hrs' ; ttA=ttA/24; titT='days';   %-- set time units:
81   ttA=ttA/30 ; titT='month'; ttA=ttA/12 ; titT='year';    titT='s'   ;  ttA=ttA/60; titT='mn' ;
82      ttA=ttA/3600; titT='hrs' ;  ttA=ttA/24; titT='days';
83      ttA=ttA/30 ; titT='month';  ttA=ttA/12 ; titT='year';
84  end  end
85  %-----------  %-----------
86    
# Line 117  for ng=1:size(list_on,2) Line 122  for ng=1:size(list_on,2)
122  %if ng == 3, flag=1*list_on(4) ; titv='Spe.Hum [g/kg]';vvA=smA; end  %if ng == 3, flag=1*list_on(4) ; titv='Spe.Hum [g/kg]';vvA=smA; end
123   if ng == 4, flag=3*list_on(8) ; titv='CFL-max'     ; vvA=cfA ; end   if ng == 4, flag=3*list_on(8) ; titv='CFL-max'     ; vvA=cfA ; end
124   if ng == 5, flag=1*list_on(2) ; titv='Eta [m]'     ; vvA=etA ; end   if ng == 5, flag=1*list_on(2) ; titv='Eta [m]'     ; vvA=etA ; end
125     if ng == 6, flag=1*list_on(7) ; titv='W-vel [m/s]' ; vvA=wmA ; end
126  %if ng == 5, flag=1*list_on(2) ; titv='Eta [mb]'; vvA=etA/100 ; end  %if ng == 5, flag=1*list_on(2) ; titv='Eta [mb]'; vvA=etA/100 ; end
127   if ng == 6, flag=1*list_on(7) ; titv='W-vel [Pa/s]'; vvA=wmA ; end  %if ng == 6, flag=1*list_on(7) ; titv='W-vel [Pa/s]'; vvA=wmA ; end
128   if ng == 7, flag=1*list_on(5) ; titv='U-vel [m/s]' ; vvA=umA ; end   if ng == 7, flag=1*list_on(5) ; titv='U-vel [m/s]' ; vvA=umA ; end
129   if ng == 8, flag=1*list_on(6) ; titv='V-vel [m/s]' ; vvA=vmA ; end   if ng == 8, flag=1*list_on(6) ; titv='V-vel [m/s]' ; vvA=vmA ; end
130   if ng == 9, flag=4*list_on(9) ; titv= 'Vort [s-1]' ; vvA=zmA ; end   if ng == 9, flag=4*list_on(9) ; titv= 'Vort [s-1]' ; vvA=zmA ; end
# Line 197  for ng=1:size(list_on,2) Line 203  for ng=1:size(list_on,2)
203      else var=squeeze(vvA(:,nv,:)); end      else var=squeeze(vvA(:,nv,:)); end
204      for n=1:Nexp,      for n=1:Nexp,
205        plot(ttA(isA(n):ieA(n),n),var(isA(n):ieA(n),n),char(linA(n)));        plot(ttA(isA(n):ieA(n),n),var(isA(n):ieA(n),n),char(linA(n)));
206          % LL(n)=plot(ttA(isA(n):ieA(n),n),var(isA(n):ieA(n),n),char(linA(n)));
207        if n == 1, hold on ; end ;        if n == 1, hold on ; end ;
208        ttmn=sprintf([ttmn,' %2.1e ;'],dd(nv,n));        ttmn=sprintf([ttmn,' %2.1e ;'],dd(nv,n));
209        ttav=sprintf([ttav,' %3.2e ;'],av(nv,n));        ttav=sprintf([ttav,' %3.2e ;'],av(nv,n));
210      end ; hold off ;      end ; hold off ;
211        % set(LL(1),'LineWidth',2);
212      if ttax1 < ttax2, AA=axis; axis([ttax1 ttax2 AA(3:4)]); end;      if ttax1 < ttax2, AA=axis; axis([ttax1 ttax2 AA(3:4)]); end;
213      grid ;      grid ;
214      if np == 1, legend(namLg,'Location','best'); end      if np == 1, legend(namLg,'Location','best'); end
# Line 221  for ng=1:size(list_on,2) Line 229  for ng=1:size(list_on,2)
229     if length(ttyax) > 0, ylabel(ttyax); end     if length(ttyax) > 0, ylabel(ttyax); end
230    end ;    end ;
231   if ngEn < 4 & list_on(10) == 1,   if ngEn < 4 & list_on(10) == 1,
232    titv='surf-Cor : ' ; vvA=scA ;    titv1='surf-Cor : ' ; vvA=scA ;
233    dd=squeeze(max(vvA)-min(vvA)); av=squeeze(mean(vvA));    dd=squeeze(max(vvA)-min(vvA)); av=squeeze(mean(vvA));
234     if Nexp == 1, av=av'; dd=dd'; end ;     if Nexp == 1, av=av'; dd=dd'; end ;
235    for nv=1:ngEx-ngEk,    for nv=1:ngEx-ngEk,
236      vvM=mean(vvA(max(2,isA(n)):ieA(n),nv,:),1);      vvM=mean(vvA(max(2,isA(n)):ieA(n),nv,:),1);
237      subplot(100*ngEx+10+nv+2); ttmn=' Mx-mn:'; ttav=' Av:';      subplot(100*ngEx+10+nv+2); ttmn=' Mx-mn:'; ttav=' Av:';
238        if nv == 1, titv2=' \theta'; else titv2=' S '; end
239        if list_log(ng) > 0,
240          var=zeros(nrec,Nexp); for n=1:Nexp, var(isA(n):ieA(n),n)=vvA(isA(n):ieA(n),nv,n); end
241          if min(var(:))*max(var(:)) >= 0, msk=var; ttyax='log';
242            var(find(msk==0))=1; var=log10(abs(var)); var(find(msk==0))=NaN;
243          end
244        else var=squeeze(vvA(:,nv,:)); end
245      for n=1:Nexp,      for n=1:Nexp,
246        plot(ttA(isA(n):ieA(n),n),var(isA(n):ieA(n),n),char(linA(n)));        plot(ttA(isA(n):ieA(n),n),var(isA(n):ieA(n),n),char(linA(n)));
247        if n == 1, hold on ; end ;        if n == 1, hold on ; end ;
# Line 235  for ng=1:size(list_on,2) Line 250  for ng=1:size(list_on,2)
250      end ; hold off ;      end ; hold off ;
251      if ttax1 < ttax2, AA=axis; axis([ttax1 ttax2 AA(3:4)]); end;      if ttax1 < ttax2, AA=axis; axis([ttax1 ttax2 AA(3:4)]); end;
252      grid ;      grid ;
253      if nv == 1, title(['mean ',titv,' \theta']); end     %title(['mean ',titv1,titv2]);
254      if nv == 2, title(['mean ',titv,' S ']); end      title(['mean ',titv1,titv2,' ; ',ttav]);
255  %--  %--
256      titX='   '; for n=1:Nexp, titX=[titX,sprintf(' %4.3e',vvM(n))]; end      titX='   '; for n=1:Nexp, titX=[titX,sprintf(' %4.3e,',vvM(n))]; end
257      AA=axis; TM=text(AA(2),AA(3)+0.1*(AA(4)-AA(3)),titX);      AA=axis; TM=text(AA(2),AA(3)+0.1*(AA(4)-AA(3)),titX(1:end-1));
258      set(TM,'HorizontalAlignment','right');      set(TM,'HorizontalAlignment','right');
259  %--  %--
260    end ;    end ;
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22