*IF DEF,W06_1A FINDPIC.2 ! FINDPIC.3 ! Description: FINDPIC.4 ! This subroutine is part of the wavetrain diagnostic output code FINDPIC.5 ! developed by Anne Guillaume at MeteoFrance and ECMWF. FINDPIC.6 ! Introduced into the unified wave moel at UM4.1 FINDPIC.7 ! FINDPIC.8 ! Method: FINDPIC.9 ! FINDPIC.10 ! FINDPIC.11 ! FINDPIC.12 ! Current Code Owner: Martin Holt FINDPIC.13 ! FINDPIC.14 ! History: FINDPIC.15 ! Version Date Comment FINDPIC.16 ! ------- ---- ------- FINDPIC.17 ! UM4.1 June 1996 Code introduced to UM. M Holt FINDPIC.18 ! FINDPIC.19 ! Code Description: FINDPIC.20 ! Language: FORTRAN 77 + common extensions. FINDPIC.21 ! FINDPIC.22 !- End of header FINDPIC.23 FINDPIC.24SUBROUTINE FINDPIC(KNRJ,KBLO,KJS,KJL,KANG,KFRE,KPICA,KPICF) 1FINDPIC.25 C FINDPIC.26 C**** *FINDPIC* - ROUTINE TO FIND PEAK OF 2D-SPECTRA FINDPIC.27 C FINDPIC.28 C A.GUILLAUME METEO-FRANCE 20/08/91 FINDPIC.29 C A.GUILLAUME ECMWF 30/06/92 FINDPIC.30 C FINDPIC.31 C PURPOSE. FINDPIC.32 C -------- FINDPIC.33 C FINDPIC.34 C *FINDPIC* FIND PEAK DIRECTION AND FREQUENCY FINDPIC.35 C FINDPIC.36 C** INTERFACE. FINDPIC.37 C ---------- FINDPIC.38 C FINDPIC.39 C *CALL* *FINDPIC(KNRJ,KBLO,KJS,KJL,KANG,KFRE,KPICA,KPICF) FINDPIC.40 C FINDPIC.41 C I/ *KNRJ* - 2D WAVE SPECTRA (*FACTOR TO BE INTEGER) FINDPIC.42 C I/ *KBLO* - DIMENSION OF ONE BLOCK. FINDPIC.43 C I/ *KJS* - INDEX OF FIRST POINT OF BLOCK TO USE. FINDPIC.44 C I/ *KJL* - INDEX OF LAST POINT OF BLOCK TO USE. FINDPIC.45 C I/ *KANG* - NUMBER OF DIRECTIONS. FINDPIC.46 C I/ *KFRE* - NUMBER OF FREQUENCIES. FINDPIC.47 C /O *KPICA* - INDEX OF PEAK DIRECTION. FINDPIC.48 C /O *KPICF* - INDEX OF PEAK FREQUENCY. FINDPIC.49 C FINDPIC.50 C METHOD. FINDPIC.51 C ------- FINDPIC.52 C FINDPIC.53 C FINDPIC.54 C EXTERNALS. FINDPIC.55 C ---------- FINDPIC.56 C NONE. FINDPIC.57 C FINDPIC.58 C REFERENCE. FINDPIC.59 C ---------- FINDPIC.60 C FINDPIC.61 C NONE. FINDPIC.62 C FINDPIC.63 DIMENSION KNRJ(KBLO,KANG,KFRE) FINDPIC.64 DIMENSION KPICA(KBLO),KPICF(KBLO) FINDPIC.65 C..WORKING ARRAYS: FINDPIC.66 DIMENSION INRJMX(KBLO) FINDPIC.67 C FINDPIC.68 C 1. CHERCHE LE MAX D'ENERGIE. FINDPIC.69 C ------- -- --- ---------- FINDPIC.70 C FINDPIC.71 100 CONTINUE FINDPIC.72 DO 101 J=KJS,KJL FINDPIC.73 INRJMX(J)=0 FINDPIC.74 101 CONTINUE FINDPIC.75 DO 102 JFRE=1,KFRE FINDPIC.76 DO 102 JANG=1,KANG FINDPIC.77 DO 102 J=KJS,KJL FINDPIC.78 INRJMX(J)=MAX0(KNRJ(J,JANG,JFRE),INRJMX(J)) FINDPIC.79 102 CONTINUE FINDPIC.80 C FINDPIC.81 C 2. ISOLE UN PIC. FINDPIC.82 C ----- -- --- FINDPIC.83 C FINDPIC.84 200 CONTINUE FINDPIC.85 DO 201 J=KJS,KJL FINDPIC.86 KPICF(J)=1 FINDPIC.87 KPICA(J)=1 FINDPIC.88 201 CONTINUE FINDPIC.89 DO 202 JFRE=1,KFRE FINDPIC.90 DO 202 JANG=1,KANG FINDPIC.91 DO 202 J=KJS,KJL FINDPIC.92 IF(KNRJ(J,JANG,JFRE).EQ.INRJMX(J)) THEN FINDPIC.93 KPICA(J)=JANG FINDPIC.94 KPICF(J)=JFRE FINDPIC.95 ENDIF FINDPIC.96 202 CONTINUE FINDPIC.97 RETURN FINDPIC.98 END FINDPIC.99 *ENDIF FINDPIC.100