/[MITgcm]/MITgcm/pkg/ptracers/ptracers_mnc_init.F
ViewVC logotype

Annotation of /MITgcm/pkg/ptracers/ptracers_mnc_init.F

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


Revision 1.2 - (hide annotations) (download)
Sat Sep 4 03:08:03 2004 UTC (19 years, 9 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint55b_post, checkpoint55, checkpoint54f_post, checkpoint55a_post
Changes since 1.1: +10 -3 lines
 o add ability to output ptracers timeave fields to NetCDF
   - small cleanups to verification/dic_example

1 edhill 1.2 C $Header: /u/gcmpack/MITgcm/pkg/ptracers/ptracers_mnc_init.F,v 1.1 2004/09/03 20:10:47 edhill Exp $
2 edhill 1.1 C $Name: $
3    
4     #include "PTRACERS_OPTIONS.h"
5    
6     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7     CBOP
8     C !ROUTINE:
9    
10     C !INTERFACE:
11     SUBROUTINE PTRACERS_MNC_INIT( myThid )
12    
13     C !DESCRIPTION:
14     C Initialize PTRACERS parameters, read in data.ptracers
15    
16     C !USES:
17     IMPLICIT NONE
18     #include "SIZE.h"
19     #include "EEPARAMS.h"
20     #include "PTRACERS_SIZE.h"
21     #include "PTRACERS.h"
22     #include "PARAMS.h"
23    
24     C !INPUT PARAMETERS:
25     INTEGER myThid
26     CEOP
27    
28     #ifdef ALLOW_PTRACERS
29     #ifdef ALLOW_MNC
30    
31     C !LOCAL VARIABLES:
32 edhill 1.2 INTEGER ii, jj, ilnb
33     character*(80) name
34 edhill 1.1 character*(MAX_LEN_MBUF) msgbuf
35    
36     C Functions
37     integer ILNBLNK
38    
39 edhill 1.2 DO ii = 1,PTRACERS_numInUse
40    
41     DO jj = 1,80
42     name(jj:jj) = ' '
43     ENDDO
44 edhill 1.1
45     C Create the tracer
46     ilnb = ILNBLNK(PTRACERS_names(ii))
47     IF (ilnb .LT. 1) THEN
48     write(msgbuf,'(2a,i3,2a)') 'PTRACERS_MNC_INIT ERROR: ',
49     & 'PTRACERS_names(',ii,') is blank--please name it in',
50     & ' ''data.ptracers'''
51     CALL print_error(msgbuf, mythid)
52     stop 'ABNORMAL END: S/R PTRACERS_MNC_INIT'
53     ENDIF
54     CALL MNC_CW_ADD_VNAME(
55     & PTRACERS_names(ii),'Cen_xy_Hn__C__t',4,5,myThid)
56 edhill 1.2 WRITE(name,'(a,a)') 'surf_',PTRACERS_names(ii)(1:ilnb)
57     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__-__t',3,4,myThid)
58 edhill 1.1
59     C Add the long_name attribute if its defined
60     ilnb = ILNBLNK(PTRACERS_long_names(ii))
61     IF (ilnb .GT. 0) THEN
62     CALL MNC_CW_ADD_VATTR_TEXT(PTRACERS_names(ii),1,'long_name',
63     & PTRACERS_long_names(ii),myThid)
64     ENDIF
65    
66     C Add the units attribute if its defined
67     ilnb = ILNBLNK(PTRACERS_units(ii))
68     IF (ilnb .GT. 0) THEN
69     CALL MNC_CW_ADD_VATTR_TEXT(PTRACERS_names(ii),1,'units',
70     & PTRACERS_units(ii),myThid)
71     ENDIF
72    
73     ENDDO
74    
75     #endif /* ALLOW_MNC */
76     #endif /* ALLOW_PTRACERS */
77    
78     RETURN
79     END
80    

  ViewVC Help
Powered by ViewVC 1.1.22